Re: [gdal-dev] extract vector/raster data from GeoPDF
Joaquim, Doing the format conversion from .ps to pdf is one thing (and there are several ways to do it), but embedding the georeferencing in the PDF to make it a GeoPDF is the interesting bit. I haven't found any open source project capable of do that. Or even capable of writing the PDF elements (e.g. Frames) necessary for writing the GeoPDF information. But I remain hopeful... Brent Joaquim Luis wrote: Brent Fraser wrote: I see the example calculation of the transformation matrix, and a statement I created a GeoPDF by running the Postscript file through Ghostscript to create a Postscript file that looks like this But I don't see where the georeferencing is written into PDF to make it a GeoPDF. Is it a Ghostscript command line? Brent, One only has to convert it to pdf using ghostscript. I did it with GMT's ps2raster, like that ps2raster quad.ps -Tf -A Joaquim Luis Brent Klokan Petr Přidal wrote: Hi, There is a great blog post (and the linked worked example post with details): http://geopdf.blogspot.com/2009/02/geopdf-map-for-worked-example.html It shows you how to create geopdf via GhostScript - so there is already a practical open-source example how to encode the georeference into the PDF/PS according the OGC standard - for use in Acrobat Reader. To add support for such tag in MapServer, which generates pdf dynamicaly via pdflib, should not be totally problematic. Decoding is not as hard either, there are nice libraries like poppler (http://poppler.freedesktop.org/), which allows you to parse vectors (and convert them to SVG for example) or rasterize the PDF files (into TIFF,...) via Cairo. The work is in assigning correct geographic coordinates to the coordinate system internally used in PDF files and especially write the bridge to the outside world (with GDAL/OGR). I am afraid that authors of the GeoPDF standard would not like this, as it seems that the idea of GeoPDF is see it in the Acrobat, print it, but that's all. At least I think so, because they discontinued their Geopdf2geotiff product and all the conversion tools are just one way - into GeoPDF. Please correct me... Anyway, in this moment you can quite easily use utility like pdfimages to extract full quality image tiles from any GeoPDF (like those from USGS) and merge it based on their location in PDF into one GDAL file via VRT (gdalbuildvrt) with a bit of hacking. This is what I did for my favorite USGS DRG of Grand Canyon ;-). Look at: http://klokan.mzk.cz/~klokan/geopdf/ - soon I will update the MapTiler.org overlay examples... Unfortunately all PDF parsing libraries I know are GPL, and that means we can't use them for the gdal driver - because of the license issues. But to create a GPL utility for converting GeoPDF to anything what GDAL/OGR supports should be OK. Poppler can be the best base of such GDAL-based utility for reading/rasterizing of the GeoPDF files. Now just find a sponsor and time to make it ;-). Best, Klokan Petr Pridal ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] extract vector/raster data from GeoPDF
Brent, I was also digging a bit on this matter. It's true that the blogger only opens our appetite as there is no more details on how to do it generally. I mean, for UTM maps it has become fairly obvious how one can do it. The relevant part in the PS file is [ {ThisPage} /LGIDict /Description Title /CTM [(35.28267) (0) (0) (35.28267) (205188.64) (3207094.8)] /Projection /Description (WGS 84 UTM 16N) /ProjectionType (TC) /Datum (WGE) /CentralMeridian (-87.0) /OriginLatitude (0.0) /FalseEasting (50.0) /FalseNorthing (0.0) /ScaleFactor (0.999600) /Type /Projection /Type /LGIDict /Version (2.1) /PUT pdfmark However, there is no more info on how to program in this TerraGo way A next post shads a bit more lite on the matter. http://geopdf.blogspot.com/2009/02/geopdf-and-geops-with-adobe-style.html As you can see in the quada.ps file there is an alternative Adobe way, which encodes the referencing as % embed georegistation info [ {ThisPage} /VP [ /Type /Viewport /BBox[0 0 dh dv] /Name Title /Measure /Type /Measure /Subtype /GEO /Bounds[0 0 0 1 1 1 1 0] /GPTS[29.0 -90.0 30.0 -90.0 30.0 -89.0 29.0 -89.0] /LPTS[0.024324 0.039461 0.051689 0.978774 0.975675 0.961136 0.957432 0.022122] /GCS /Type /PROJCS /WKT (PROJCS[WGS_1984_UTM_Zone_16N,GEOGCS[GCS_WGS_1984,DATUM[D_WGS_1984,SPHEROID[WGS_1984,6378137,298.257223563]],PRIMEM[Greenwich,0],UNIT[Degree,0.017453292519943295]],PROJECTION[Transverse_Mercator],PARAMETER[latitude_of_origin,0],PARAMETER[central_meridian,-87],PARAMETER[scale_factor,0.9996],PARAMETER[false_easting,50],PARAMETER[false_northing,0],UNIT[Meter,1]]) ] /PUT pdfmark Now this is more familiar, except for the /LPTS[0.024324 0.039461 0.051689 0.978774 0.975675 0.961136 0.957432 0.022122] part that, contrary to the promise has not been explained where it comes from. I think I could implement the first encoding type in GMT's ps2raster (if I knew all the specs) Second type would be even easier ... if we had a way to translate the GMT projection syntax into WKT, which we don't. Joaquim Joaquim, Doing the format conversion from .ps to pdf is one thing (and there are several ways to do it), but embedding the georeferencing in the PDF to make it a GeoPDF is the interesting bit. I haven't found any open source project capable of do that. Or even capable of writing the PDF elements (e.g. Frames) necessary for writing the GeoPDF information. But I remain hopeful... Brent Joaquim Luis wrote: Brent Fraser wrote: I see the example calculation of the transformation matrix, and a statement I created a GeoPDF by running the Postscript file through Ghostscript to create a Postscript file that looks like this But I don't see where the georeferencing is written into PDF to make it a GeoPDF. Is it a Ghostscript command line? Brent, One only has to convert it to pdf using ghostscript. I did it with GMT's ps2raster, like that ps2raster quad.ps -Tf -A Joaquim Luis Brent Klokan Petr Přidal wrote: Hi, There is a great blog post (and the linked worked example post with details): http://geopdf.blogspot.com/2009/02/geopdf-map-for-worked-example.html It shows you how to create geopdf via GhostScript - so there is already a practical open-source example how to encode the georeference into the PDF/PS according the OGC standard - for use in Acrobat Reader. To add support for such tag in MapServer, which generates pdf dynamicaly via pdflib, should not be totally problematic. Decoding is not as hard either, there are nice libraries like poppler (http://poppler.freedesktop.org/), which allows you to parse vectors (and convert them to SVG for example) or rasterize the PDF files (into TIFF,...) via Cairo. The work is in assigning correct geographic coordinates to the coordinate system internally used in PDF files and especially write the bridge to the outside world (with GDAL/OGR). I am afraid that authors of the GeoPDF standard would not like this, as it seems that the idea of GeoPDF is see it in the Acrobat, print it, but that's all. At least I think so, because they discontinued their Geopdf2geotiff product and all the conversion tools are just one way - into GeoPDF. Please correct me... Anyway, in this moment you can quite easily use utility like pdfimages to extract full quality image tiles from any GeoPDF (like those from USGS) and merge it based on their location in PDF into one GDAL file via VRT (gdalbuildvrt) with a bit of hacking. This is what I did for my favorite USGS DRG of Grand Canyon ;-). Look at: http://klokan.mzk.cz/~klokan/geopdf/ - soon I will update the MapTiler.org overlay examples... Unfortunately all PDF parsing libraries I know are GPL, and that means we can't use them for the gdal driver - because of the license issues. But to create a GPL utility for converting GeoPDF to anything what GDAL/OGR supports should be
Re: [gdal-dev] extract vector/raster data from GeoPDF
A bit of googling determines that the 'Adobe way' structure, together with its elements, is defined in the ISO PDF definition - ISO 32000-1 (2008). Looking at this with a PostScript Developer hat, the /LPTS has to be associated with scaling / transforming the GPTS coordinates to the page coordinate system. There is a need to refer to the standard! Best wishes, Peter Joaquim Luis wrote: Brent, I was also digging a bit on this matter. It's true that the blogger only opens our appetite as there is no more details on how to do it generally. I mean, for UTM maps it has become fairly obvious how one can do it. The relevant part in the PS file is [ {ThisPage} /LGIDict /Description Title /CTM [(35.28267) (0) (0) (35.28267) (205188.64) (3207094.8)] /Projection /Description (WGS 84 UTM 16N) /ProjectionType (TC) /Datum (WGE) /CentralMeridian (-87.0) /OriginLatitude (0.0) /FalseEasting (50.0) /FalseNorthing (0.0) /ScaleFactor (0.999600) /Type /Projection /Type /LGIDict /Version (2.1) /PUT pdfmark However, there is no more info on how to program in this TerraGo way A next post shads a bit more lite on the matter. http://geopdf.blogspot.com/2009/02/geopdf-and-geops-with-adobe-style.html As you can see in the quada.ps file there is an alternative Adobe way, which encodes the referencing as % embed georegistation info [ {ThisPage} /VP [ /Type /Viewport /BBox[0 0 dh dv] /Name Title /Measure /Type /Measure /Subtype /GEO /Bounds[0 0 0 1 1 1 1 0] /GPTS[29.0 -90.0 30.0 -90.0 30.0 -89.0 29.0 -89.0] /LPTS[0.024324 0.039461 0.051689 0.978774 0.975675 0.961136 0.957432 0.022122] /GCS /Type /PROJCS /WKT (PROJCS[WGS_1984_UTM_Zone_16N,GEOGCS[GCS_WGS_1984,DATUM[D_WGS_1984,SPHEROID[WGS_1984,6378137,298.257223563]],PRIMEM[Greenwich,0],UNIT[Degree,0.017453292519943295]],PROJECTION[Transverse_Mercator],PARAMETER[latitude_of_origin,0],PARAMETER[central_meridian,-87],PARAMETER[scale_factor,0.9996],PARAMETER[false_easting,50],PARAMETER[false_northing,0],UNIT[Meter,1]]) ] /PUT pdfmark Now this is more familiar, except for the /LPTS[0.024324 0.039461 0.051689 0.978774 0.975675 0.961136 0.957432 0.022122] part that, contrary to the promise has not been explained where it comes from. I think I could implement the first encoding type in GMT's ps2raster (if I knew all the specs) Second type would be even easier ... if we had a way to translate the GMT projection syntax into WKT, which we don't. Joaquim Joaquim, Doing the format conversion from .ps to pdf is one thing (and there are several ways to do it), but embedding the georeferencing in the PDF to make it a GeoPDF is the interesting bit. I haven't found any open source project capable of do that. Or even capable of writing the PDF elements (e.g. Frames) necessary for writing the GeoPDF information. But I remain hopeful... Brent Joaquim Luis wrote: Brent Fraser wrote: I see the example calculation of the transformation matrix, and a statement I created a GeoPDF by running the Postscript file through Ghostscript to create a Postscript file that looks like this But I don't see where the georeferencing is written into PDF to make it a GeoPDF. Is it a Ghostscript command line? Brent, One only has to convert it to pdf using ghostscript. I did it with GMT's ps2raster, like that ps2raster quad.ps -Tf -A Joaquim Luis Brent Klokan Petr Přidal wrote: Hi, There is a great blog post (and the linked worked example post with details): http://geopdf.blogspot.com/2009/02/geopdf-map-for-worked-example.html It shows you how to create geopdf via GhostScript - so there is already a practical open-source example how to encode the georeference into the PDF/PS according the OGC standard - for use in Acrobat Reader. To add support for such tag in MapServer, which generates pdf dynamicaly via pdflib, should not be totally problematic. Decoding is not as hard either, there are nice libraries like poppler (http://poppler.freedesktop.org/), which allows you to parse vectors (and convert them to SVG for example) or rasterize the PDF files (into TIFF,...) via Cairo. The work is in assigning correct geographic coordinates to the coordinate system internally used in PDF files and especially write the bridge to the outside world (with GDAL/OGR). I am afraid that authors of the GeoPDF standard would not like this, as it seems that the idea of GeoPDF is see it in the Acrobat, print it, but that's all. At least I think so, because they discontinued their Geopdf2geotiff product and all the conversion tools are just one way - into GeoPDF. Please correct me... Anyway, in this moment you can quite easily use utility like pdfimages to extract full quality image tiles from any GeoPDF (like those from USGS) and merge it based on their location in PDF into one GDAL file via VRT (gdalbuildvrt) with a bit of
Re: [gdal-dev] extract vector/raster data from GeoPDF
Peter J Halls wrote: A bit of googling determines that the 'Adobe way' structure, together with its elements, is defined in the ISO PDF definition - ISO 32000-1 (2008). Looking at this with a PostScript Developer hat, the /LPTS has to be associated with scaling / transforming the GPTS coordinates to the page coordinate system. There is a need to refer to the standard! Yes, And the TerraGo way comes possibly from this document GeoPDF Encoding Best Practice Version 2.2 http://portal.opengeospatial.org/files/?artifact_id=2 ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Shapefile to Arc/info ascii
Hi all, I've been looking around in your mailing list searching for a way to convert a vector shapefile to a raster Arc/Info ascii format but haven't found any applicable example. It seems like rasterize could be a solution in the GDAL library but I don't get how. The steps I'm aiming to perform are the ones below, can GDAL do that? *Put together a mosaic of several shape files into one single shape file. *Crop a rectangular area of the new shape file. *Rasterize the new shapefile file (all shapes should be used) with a specific resolution in x and y and specify which attribute column that should be rasterized. Thank you, Jonas Malmqvist -- ¤¤ Jonas Malmqvist FS Dynamics AB Jungmansgatan 31 SE-413 11 Göteborg Sweden Phone: +46 (0)31 720 71 45 Fax: +46 (0)31 761 99 49 E-mail: jonas.malmqv...@fsdynamics.se ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] extract vector/raster data from GeoPDF
Joaquim, Interesting! As Peter pointed out, the LPTS tag is the transformation matrix to go from projected coords to page coords, and is created by forming the UTM-page coordinates tranformation (from the information in the CTM tag below) and inverting it (there's a brief mention near the end of http://geopdf.blogspot.com/2009/02/georegistration-worked-example.html) Since the TerraGo way has been superceded by the Adobe/ISO way (thank goodness!), it makes it easier to identify the reference system definition in GDAL, so there may be some hope for GDAL and GeoPDF. Now I need to think about my requirements for reading and writing GeoPDFs. Write a Postscript driver for GDAL? Figure out how to write the equivalent info into a PDF using Mapserver/Cairo? Yikes! Thanks, Brent Joaquim Luis wrote: Brent, I was also digging a bit on this matter. It's true that the blogger only opens our appetite as there is no more details on how to do it generally. I mean, for UTM maps it has become fairly obvious how one can do it. The relevant part in the PS file is [ {ThisPage} /LGIDict /Description Title /CTM [(35.28267) (0) (0) (35.28267) (205188.64) (3207094.8)] /Projection /Description (WGS 84 UTM 16N) /ProjectionType (TC) /Datum (WGE) /CentralMeridian (-87.0) /OriginLatitude (0.0) /FalseEasting (50.0) /FalseNorthing (0.0) /ScaleFactor (0.999600) /Type /Projection /Type /LGIDict /Version (2.1) /PUT pdfmark However, there is no more info on how to program in this TerraGo way A next post shads a bit more lite on the matter. http://geopdf.blogspot.com/2009/02/geopdf-and-geops-with-adobe-style.html As you can see in the quada.ps file there is an alternative Adobe way, which encodes the referencing as % embed georegistation info [ {ThisPage} /VP [ /Type /Viewport /BBox[0 0 dh dv] /Name Title /Measure /Type /Measure /Subtype /GEO /Bounds[0 0 0 1 1 1 1 0] /GPTS[29.0 -90.0 30.0 -90.0 30.0 -89.0 29.0 -89.0] /LPTS[0.024324 0.039461 0.051689 0.978774 0.975675 0.961136 0.957432 0.022122] /GCS /Type /PROJCS /WKT (PROJCS[WGS_1984_UTM_Zone_16N,GEOGCS[GCS_WGS_1984,DATUM[D_WGS_1984,SPHEROID[WGS_1984,6378137,298.257223563]],PRIMEM[Greenwich,0],UNIT[Degree,0.017453292519943295]],PROJECTION[Transverse_Mercator],PARAMETER[latitude_of_origin,0],PARAMETER[central_meridian,-87],PARAMETER[scale_factor,0.9996],PARAMETER[false_easting,50],PARAMETER[false_northing,0],UNIT[Meter,1]]) ] /PUT pdfmark Now this is more familiar, except for the /LPTS[0.024324 0.039461 0.051689 0.978774 0.975675 0.961136 0.957432 0.022122] part that, contrary to the promise has not been explained where it comes from. I think I could implement the first encoding type in GMT's ps2raster (if I knew all the specs) Second type would be even easier ... if we had a way to translate the GMT projection syntax into WKT, which we don't. Joaquim Joaquim, Doing the format conversion from .ps to pdf is one thing (and there are several ways to do it), but embedding the georeferencing in the PDF to make it a GeoPDF is the interesting bit. I haven't found any open source project capable of do that. Or even capable of writing the PDF elements (e.g. Frames) necessary for writing the GeoPDF information. But I remain hopeful... Brent Joaquim Luis wrote: Brent Fraser wrote: I see the example calculation of the transformation matrix, and a statement I created a GeoPDF by running the Postscript file through Ghostscript to create a Postscript file that looks like this But I don't see where the georeferencing is written into PDF to make it a GeoPDF. Is it a Ghostscript command line? Brent, One only has to convert it to pdf using ghostscript. I did it with GMT's ps2raster, like that ps2raster quad.ps -Tf -A Joaquim Luis Brent Klokan Petr Přidal wrote: Hi, There is a great blog post (and the linked worked example post with details): http://geopdf.blogspot.com/2009/02/geopdf-map-for-worked-example.html It shows you how to create geopdf via GhostScript - so there is already a practical open-source example how to encode the georeference into the PDF/PS according the OGC standard - for use in Acrobat Reader. To add support for such tag in MapServer, which generates pdf dynamicaly via pdflib, should not be totally problematic. Decoding is not as hard either, there are nice libraries like poppler (http://poppler.freedesktop.org/), which allows you to parse vectors (and convert them to SVG for example) or rasterize the PDF files (into TIFF,...) via Cairo. The work is in assigning correct geographic coordinates to the coordinate system internally used in PDF files and especially write the bridge to the outside world (with GDAL/OGR). I am afraid that authors of the GeoPDF standard would not like this, as it seems that the idea of GeoPDF is see it in the Acrobat, print
[gdal-dev] about one-dimensional vector dataset
Hi, How to get values of one-dimensional vector dataset like this [image: 1.png] the X_Axis dataset is showed in HDF Explorer like this: [image: 2.png] I could not access these dataset via GetMetadata(SUBDATASETS) method, Anybody could tell me how to access these datasets? Best wishes, -- JoSn 2.png1.png___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Shapefile to Arc/info ascii
Jonas Malmqvist wrote: Hi all, I've been looking around in your mailing list searching for a way to convert a vector shapefile to a raster Arc/Info ascii format but haven't found any applicable example. It seems like rasterize could be a solution in the GDAL library but I don't get how. The steps I'm aiming to perform are the ones below, can GDAL do that? *Put together a mosaic of several shape files into one single shape file. Jonas, ogr2ogr with the -update and -append switches can be used to concatenate several compatible shapefiles into one. The source files should have essentially the same set of attributes and a compatible geometry type. *Crop a rectangular area of the new shape file. ogr2ogr can extract only the features that intersect a particular rectangle using the -spat switch. However, the features that are partly inside, and partly outside are not clipped to the rectangle - you get the whole feature. *Rasterize the new shapefile file (all shapes should be used) with a specific resolution in x and y and specify which attribute column that should be rasterized. The gdal_rasterize command can handle the rasterization (use -a to select the attribute field to burn in), but you will have to precreate the file with the desired georeferenced region and resolution. This can often be accomplished through suitable use of gdal_translate or it can be easily done in a gdal script (python, etc). Best regards, -- ---+-- I set the clouds in motion - turn up | Frank Warmerdam, warmer...@pobox.com light and sound - activate the windows | http://pobox.com/~warmerdam and watch the world go round - Rush| Geospatial Programmer for Rent ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] extract vector/raster data from GeoPDF
Brent Fraser wrote: Joaquim, Interesting! As Peter pointed out, the LPTS tag is the transformation matrix to go from projected coords to page coords, and is created by forming the UTM-page coordinates tranformation (from the information in the CTM tag below) and inverting it (there's a brief mention near the end of http://geopdf.blogspot.com/2009/02/georegistration-worked-example.html) Since the TerraGo way has been superceded by the Adobe/ISO way (thank goodness!), it makes it easier to identify the reference system definition in GDAL, so there may be some hope for GDAL and GeoPDF. Now I need to think about my requirements for reading and writing GeoPDFs. Write a Postscript driver for GDAL? Figure out how to write the equivalent info into a PDF using Mapserver/Cairo? Yikes! Brent, I am on this list more to keep learning (as the present case) than as a core GDAL user (but I use GDAL a lot in Mirone). In fact I am much more a GMT user where I also do some development. I mentioned it before here but probably it slipped away, the GMT program ps2raster can write world files and create GeoTIFFs (from system calls to GDAL) from postscript files created by GMT. I think that it wouldn't be very hard to extend ps2raster's ability to write GeoPDFs if, as I said before, there was a way to get a WKT from a -J... GMT syntax. Anyway, it wouldn't be a GDAL solution and a write-only one so I don't know what else to tell you. Joaquim ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Re: Shapefile to Arc/info ascii
Jonas Malmqvist wrote: *Put together a mosaic of several shape files into one single shape file. http://www.gdal.org/ogr/drv_shapefile.html (see the example at the bottom of the page) *Crop a rectangular area of the new shape file. *Rasterize the new shapefile file (all shapes should be used) with a specific resolution in x and y and specify which attribute column that should be rasterized. Create a blank raster file in the desired extents and resolution, e.g. via gdal_translate, see http://lists.osgeo.org/pipermail/gdal-dev/2008-February/016061.html Finally: gdal_rasterize -a ... see http://www.gdal.org/gdal_rasterize.html Hermann ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
RE: [gdal-dev] Motion: Commit access for Harsh Govind
Good Afternoon All, I have read and understood the RFC 3 GDAL committer guidelines (http://trac.osgeo.org/gdal/wiki/rfc3_commiters) and I accept the agreement. Regards, Harsh Govind SPADAC Inc. -Original Message- From: gdal-dev-boun...@lists.osgeo.org [mailto:gdal-dev-boun...@lists.osgeo.org] On Behalf Of Frank Warmerdam Sent: Wednesday, September 02, 2009 3:56 PM To: gdal-dev Subject: [gdal-dev] Motion: Commit access for Harsh Govind Motion: Extent commit access to Harsh Govind. --- PSC members, Harsh reports: I am a GIS Developer at SPADAC Inc. and responsibilities include compiling and upgrading GDAL for our products. In the past 4+ years at SPADAC I have setup automated build process for compiling and building GDAL along with its SDE (multiple versions of ArcSDE) and OCI drivers as plug-ins and C# wrapper for .NET version of our product: Signature Analyst(r) (SA) [Details: http://www.spadac.com/SignatureAnalyst/]. The build process keeps debug and release versions separate so programmers may step into the code if required during development. My upcoming tasks include creating separate builds for 32bit and 64bit machines. I have implemented multiple export formats (vector) for results of SA utilizing GDAL/ OGR. I have been involved in development of tools for reading/writing vector data using GDAL/OGR and also using C# wrapper giving me experience in both C++ and C# (SDE, SHP and PGDB formats in particular). E.g. Export SA assessment as Shapefile, Read SDE data source as factor for SA, Read PGDB data source as factor in SA, etc... I recently fixed a bug in ogrsdedatasource.cpp which would cause a failure to return a particular layer if SDE version is SDE.DEFAULT. Also, I am involved in multiple tasks relating GDAL/OGR with other Team members that include KML Super Overlay driver and supporting date query in SHP files. Harsh and SPADAC have been contributing fixes through me and Trac and now have several substantial improvements including a new KML Super Overlay driver they would like to contribute. I have suggested that Harsh become a committer so that SPADAC can more directly contribute improvements. I'll let Harsh confirm to the mailing list that he is prepared to adhere to the RFC 3 GDAL Committer Guidelines. Best regards, -- ---+-- I set the clouds in motion - turn up | Frank Warmerdam, warmer...@pobox.com light and sound - activate the windows | http://pobox.com/~warmerdam and watch the world go round - Rush| Geospatial Programmer for Rent ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Seek error accessing TIFF directory
Hello list, I just compiled gdal 1.6.2 on my fedora 10 box from source. I'm in the process of evaluating some 2009 imagery in geotiff format. I've tried compiling with external libs and internal options with the same result. # ./configure --with-libtiff=internal --with-geotiff=internal When processing the tiff with mapnik rendering library I get: UserWarning: TIFFReadDirectory [myfile].tif Seek error accessing TIFF directory (encountered during parsing of layer 'imagery') I'm not sure what function of gdal is being called by mapnik but things like gdalinfo, gdalwarp seem to work without any issues its just when rendering a tile with mapnik. When I compile against gdal version 1.5.3 (package install with fedora 10) everything seems to work as expected. Any suggestions or is this a bug in the newer version? -Matt ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] about one-dimensional vector dataset
jor sion wrote: Hi, How to get values of one-dimensional vector dataset like this 1.png the X_Axis dataset is showed in HDF Explorer like this: 2.png I could not access these dataset via GetMetadata(SUBDATASETS) method, Anybody could tell me how to access these datasets? JoSn, The graphic isn't much for me to go on. Perhaps you could provide a pointer to a dataset that we could examine? Best regards, -- ---+-- I set the clouds in motion - turn up | Frank Warmerdam, warmer...@pobox.com light and sound - activate the windows | http://pobox.com/~warmerdam and watch the world go round - Rush| Geospatial Programmer for Rent ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Motion: Commit access for Harsh Govind
+1. Welcome :) Daniel Harsh Govind wrote: Good Afternoon All, I have read and understood the RFC 3 GDAL committer guidelines (http://trac.osgeo.org/gdal/wiki/rfc3_commiters) and I accept the agreement. Regards, Harsh Govind SPADAC Inc. -Original Message- From: gdal-dev-boun...@lists.osgeo.org [mailto:gdal-dev-boun...@lists.osgeo.org] On Behalf Of Frank Warmerdam Sent: Wednesday, September 02, 2009 3:56 PM To: gdal-dev Subject: [gdal-dev] Motion: Commit access for Harsh Govind Motion: Extent commit access to Harsh Govind. --- PSC members, Harsh reports: I am a GIS Developer at SPADAC Inc. and responsibilities include compiling and upgrading GDAL for our products. In the past 4+ years at SPADAC I have setup automated build process for compiling and building GDAL along with its SDE (multiple versions of ArcSDE) and OCI drivers as plug-ins and C# wrapper for .NET version of our product: Signature Analyst(r) (SA) [Details: http://www.spadac.com/SignatureAnalyst/]. The build process keeps debug and release versions separate so programmers may step into the code if required during development. My upcoming tasks include creating separate builds for 32bit and 64bit machines. I have implemented multiple export formats (vector) for results of SA utilizing GDAL/ OGR. I have been involved in development of tools for reading/writing vector data using GDAL/OGR and also using C# wrapper giving me experience in both C++ and C# (SDE, SHP and PGDB formats in particular). E.g. Export SA assessment as Shapefile, Read SDE data source as factor for SA, Read PGDB data source as factor in SA, etc... I recently fixed a bug in ogrsdedatasource.cpp which would cause a failure to return a particular layer if SDE version is SDE.DEFAULT. Also, I am involved in multiple tasks relating GDAL/OGR with other Team members that include KML Super Overlay driver and supporting date query in SHP files. Harsh and SPADAC have been contributing fixes through me and Trac and now have several substantial improvements including a new KML Super Overlay driver they would like to contribute. I have suggested that Harsh become a committer so that SPADAC can more directly contribute improvements. I'll let Harsh confirm to the mailing list that he is prepared to adhere to the RFC 3 GDAL Committer Guidelines. Best regards, -- Daniel Morissette http://www.mapgears.com/ ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Motion: Commit access for Harsh Govind
Selon Frank Warmerdam warmer...@pobox.com: +1 Motion: Extent commit access to Harsh Govind. --- PSC members, Harsh reports: I am a GIS Developer at SPADAC Inc. and responsibilities include compiling and upgrading GDAL for our products. In the past 4+ years at SPADAC I have setup automated build process for compiling and building GDAL along with its SDE (multiple versions of ArcSDE) and OCI drivers as plug-ins and C# wrapper for .NET version of our product: Signature Analyst® (SA) [Details: http://www.spadac.com/SignatureAnalyst/]. The build process keeps debug and release versions separate so programmers may step into the code if required during development. My upcoming tasks include creating separate builds for 32bit and 64bit machines. I have implemented multiple export formats (vector) for results of SA utilizing GDAL/ OGR. I have been involved in development of tools for reading/writing vector data using GDAL/OGR and also using C# wrapper giving me experience in both C++ and C# (SDE, SHP and PGDB formats in particular). E.g. Export SA assessment as Shapefile, Read SDE data source as factor for SA, Read PGDB data source as factor in SA, etc... I recently fixed a bug in ogrsdedatasource.cpp which would cause a failure to return a particular layer if SDE version is SDE.DEFAULT. Also, I am involved in multiple tasks relating GDAL/OGR with other Team members that include KML Super Overlay driver and supporting date query in SHP files. Harsh and SPADAC have been contributing fixes through me and Trac and now have several substantial improvements including a new KML Super Overlay driver they would like to contribute. I have suggested that Harsh become a committer so that SPADAC can more directly contribute improvements. I'll let Harsh confirm to the mailing list that he is prepared to adhere to the RFC 3 GDAL Committer Guidelines. Best regards, -- ---+-- I set the clouds in motion - turn up | Frank Warmerdam, warmer...@pobox.com light and sound - activate the windows | http://pobox.com/~warmerdam and watch the world go round - Rush| Geospatial Programmer for Rent ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Motion: Commit access for Harsh Govind
+1 I would also be interested in having SDE for Win64 to be compiled somehow ;-) Best regards, Tamas 2009/9/2 Frank Warmerdam warmer...@pobox.com Motion: Extent commit access to Harsh Govind. --- PSC members, Harsh reports: I am a GIS Developer at SPADAC Inc. and responsibilities include compiling and upgrading GDAL for our products. In the past 4+ years at SPADAC I have setup automated build process for compiling and building GDAL along with its SDE (multiple versions of ArcSDE) and OCI drivers as plug-ins and C# wrapper for .NET version of our product: Signature Analyst® (SA) [Details: http://www.spadac.com/SignatureAnalyst/]. The build process keeps debug and release versions separate so programmers may step into the code if required during development. My upcoming tasks include creating separate builds for 32bit and 64bit machines. I have implemented multiple export formats (vector) for results of SA utilizing GDAL/ OGR. I have been involved in development of tools for reading/writing vector data using GDAL/OGR and also using C# wrapper giving me experience in both C++ and C# (SDE, SHP and PGDB formats in particular). E.g. Export SA assessment as Shapefile, Read SDE data source as factor for SA, Read PGDB data source as factor in SA, etc... I recently fixed a bug in ogrsdedatasource.cpp which would cause a failure to return a particular layer if SDE version is SDE.DEFAULT. Also, I am involved in multiple tasks relating GDAL/OGR with other Team members that include KML Super Overlay driver and supporting date query in SHP files. Harsh and SPADAC have been contributing fixes through me and Trac and now have several substantial improvements including a new KML Super Overlay driver they would like to contribute. I have suggested that Harsh become a committer so that SPADAC can more directly contribute improvements. I'll let Harsh confirm to the mailing list that he is prepared to adhere to the RFC 3 GDAL Committer Guidelines. Best regards, -- ---+-- I set the clouds in motion - turn up | Frank Warmerdam, warmer...@pobox.com light and sound - activate the windows | http://pobox.com/~warmerdamhttp://pobox.com/%7Ewarmerdam and watch the world go round - Rush| Geospatial Programmer for Rent ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] GDAL: How to create one NITF image from CIB CD
I am trying to use GDAL utilities (on a Solaris 10 system) against CIB files (from a CD) to convert them to a single NITF image. When I run gdal_translate, it doesn't seem to accept the a.toc, but it does create a one-to-one nitf image for each CIB file the command is run against. Is there a GDAL command that takes the table of contents file and processes all the CIB frames into one NITF image? Or takes all the .ntf files that get created and stitches them together into one .ntf image? Thank you, Carol Rydzak General Dynamics C4 Systems Phone: (480) 441-7405 Email: carol.ryd...@gdc4s.com This email message is for the sole use of the intended recipient(s) and may contain GDC4S confidential or privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not an intended recipient, please contact the sender by reply email and destroy all copies of the original message. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Re: Can't get 1*N dataset in HDFfile
jor sion wrote: Hi, I can't get 1*N dataset in HDF file by GDAL, The dataset looks like this: dataset.png It shows in HDF Explorer show.png Any suggestions? Well it's not a GDAL solution... But have a look at PyTables (depends NumPy) for Python. PyTables reads HDF files (4 and 5) and can do what you want. Martin Hvidberg ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] GDAL: How to create one NITF image from CIB CD
Selon Even Rouault even.roua...@mines-paris.org: Actually, I read too fast. You want to convert all the NITF tiles into a single NITF file. That is possible too of course. Just add -of NITF in the gdal_translate command line. Note that the resulting NITF will be larger in size than the sum of source NITF tiles, as CIB uses the Vector Quantization compression scheme that is only supported as read-only by GDAL. Selon Rydzak, Carol-P28503 carol.ryd...@gdc4s.com: I assume you are using a recent enough GDAL version (= 1.5.0) that has the RPFTOC driver. If so, gdalinfo on the a.toc should return a list of subdatasets (at least one). You can then select one and translate it to a GeoTIFF. It will mosaic all the underlying NITF tiles. See http://gdal.org/frmt_various.html#RPFTOC for more info. Someone recently reported problems with CIB but I've not heard any news since then. The driver was tested against CADRG datasets but should theoritically work with CIB too. I am trying to use GDAL utilities (on a Solaris 10 system) against CIB files (from a CD) to convert them to a single NITF image. When I run gdal_translate, it doesn't seem to accept the a.toc, but it does create a one-to-one nitf image for each CIB file the command is run against. Is there a GDAL command that takes the table of contents file and processes all the CIB frames into one NITF image? Or takes all the .ntf files that get created and stitches them together into one .ntf image? Thank you, Carol Rydzak General Dynamics C4 Systems Phone: (480) 441-7405 Email: carol.ryd...@gdc4s.com This email message is for the sole use of the intended recipient(s) and may contain GDC4S confidential or privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not an intended recipient, please contact the sender by reply email and destroy all copies of the original message. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Missing msvcr80.dll ?
Dear Group Thanks for GDAL I am running a script in Python, using GDAL. Would like to run it on another (faster) computer. Installed Python and GDAL When I run the program it asks for a MSVCR80.dll file Cant find it anywhere on the computer Tried to copy it from the old computer, but it won't accept to load that. What the h... is MSCVR80 and how do I get it onto my computer ??? Best Regards Martin Hvidberg ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] extract vector/raster data from GeoPDF
Brent Fraser wrote: Since the TerraGo way has been superceded by the Adobe/ISO way (thank goodness!), it makes it easier to identify the reference system definition in GDAL, so there may be some hope for GDAL and GeoPDF. Still on this matter, the GeoPDF http://www.terragotech.com/blog/quada.pdf written in the Adobe way is not recognized as a GeoPDF by my Acrobat with the TerraGo Toolbar and it crashes Chrome when trying to download it!!! ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
RE: [gdal-dev] GDAL: How to create one NITF image from CIB CD
Carol, http://www.gdal.org/gdal_translate.html says the command switch: -sds: Copy all subdatasets of this file to individual output files. Use withformats like HDF or OGDI that have subdatasets. OGDI is the old RP driver so I imagine it applies. Maybe that's what you need. Also, you can find the source code for gdal_translate.cpp in the apps directory of your gdal download. Take a look at the code and you can see exactly what's going on. If you have a debugger you can step through the code. You may have to morph the code to get exactly what you want though. Martin -Original Message- From: gdal-dev-boun...@lists.osgeo.org [mailto:gdal-dev-boun...@lists.osgeo.org] On Behalf Of Rydzak, Carol-P28503 Sent: Thursday, September 03, 2009 8:26 PM To: Even Rouault Cc: gdal-dev@lists.osgeo.org Subject: RE: [gdal-dev] GDAL: How to create one NITF image from CIB CD Hi. Thank you for getting back to me on this! Yes, I am using the 1.6.0 version of GDAL. I did what you suggested and it gives me a bus error when I run the gdalinfo command on the a.toc file. Do you recall what the problems were with the CIB data that you heard about? I have a 1 meter, 5 meter and 10 meter data sets that I tried and I got a bus error on all of them. Any other ideas? Thanks, Carol -Original Message- From: Even Rouault [mailto:even.roua...@mines-paris.org] Sent: Thursday, September 03, 2009 12:50 PM To: Rydzak, Carol-P28503 Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] GDAL: How to create one NITF image from CIB CD Selon Rydzak, Carol-P28503 carol.ryd...@gdc4s.com: I assume you are using a recent enough GDAL version (= 1.5.0) that has the RPFTOC driver. If so, gdalinfo on the a.toc should return a list of subdatasets (at least one). You can then select one and translate it to a GeoTIFF. It will mosaic all the underlying NITF tiles. See http://gdal.org/frmt_various.html#RPFTOC for more info. Someone recently reported problems with CIB but I've not heard any news since then. The driver was tested against CADRG datasets but should theoritically work with CIB too. I am trying to use GDAL utilities (on a Solaris 10 system) against CIB files (from a CD) to convert them to a single NITF image. When I run gdal_translate, it doesn't seem to accept the a.toc, but it does create a one-to-one nitf image for each CIB file the command is run against. Is there a GDAL command that takes the table of contents file and processes all the CIB frames into one NITF image? Or takes all the .ntf files that get created and stitches them together into one .ntf image? Thank you, Carol Rydzak General Dynamics C4 Systems Phone: (480) 441-7405 Email: carol.ryd...@gdc4s.com This email message is for the sole use of the intended recipient(s) and may contain GDC4S confidential or privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not an intended recipient, please contact the sender by reply email and destroy all copies of the original message. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
RE: [gdal-dev] GDAL: How to create one NITF image from CIB CD
Carol, Also see the switch: src_dataset: The source dataset name. It can be either file name, URL of data source or subdataset name for multi-dataset files. Martin -Original Message- From: gdal-dev-boun...@lists.osgeo.org [mailto:gdal-dev-boun...@lists.osgeo.org] On Behalf Of Rydzak, Carol-P28503 Sent: Thursday, September 03, 2009 8:26 PM To: Even Rouault Cc: gdal-dev@lists.osgeo.org Subject: RE: [gdal-dev] GDAL: How to create one NITF image from CIB CD Hi. Thank you for getting back to me on this! Yes, I am using the 1.6.0 version of GDAL. I did what you suggested and it gives me a bus error when I run the gdalinfo command on the a.toc file. Do you recall what the problems were with the CIB data that you heard about? I have a 1 meter, 5 meter and 10 meter data sets that I tried and I got a bus error on all of them. Any other ideas? Thanks, Carol -Original Message- From: Even Rouault [mailto:even.roua...@mines-paris.org] Sent: Thursday, September 03, 2009 12:50 PM To: Rydzak, Carol-P28503 Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] GDAL: How to create one NITF image from CIB CD Selon Rydzak, Carol-P28503 carol.ryd...@gdc4s.com: I assume you are using a recent enough GDAL version (= 1.5.0) that has the RPFTOC driver. If so, gdalinfo on the a.toc should return a list of subdatasets (at least one). You can then select one and translate it to a GeoTIFF. It will mosaic all the underlying NITF tiles. See http://gdal.org/frmt_various.html#RPFTOC for more info. Someone recently reported problems with CIB but I've not heard any news since then. The driver was tested against CADRG datasets but should theoritically work with CIB too. I am trying to use GDAL utilities (on a Solaris 10 system) against CIB files (from a CD) to convert them to a single NITF image. When I run gdal_translate, it doesn't seem to accept the a.toc, but it does create a one-to-one nitf image for each CIB file the command is run against. Is there a GDAL command that takes the table of contents file and processes all the CIB frames into one NITF image? Or takes all the .ntf files that get created and stitches them together into one .ntf image? Thank you, Carol Rydzak General Dynamics C4 Systems Phone: (480) 441-7405 Email: carol.ryd...@gdc4s.com This email message is for the sole use of the intended recipient(s) and may contain GDC4S confidential or privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not an intended recipient, please contact the sender by reply email and destroy all copies of the original message. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
RE: [gdal-dev] GDAL: How to create one NITF image from CIB CD
Carol, If you write code the way to go about mosaicing all rpf frame files into one image would be to create a GDAL MEM dataset that is the width and height of an entire RPF Boundary. Then loop through each Frame file in the Boundary, open up each Frame file with the GDAL NITF driver and copy the bytes from the NITF file into the correct tile location in the MEM grid. Then save the GDAL MEM file out to the file format of your choice. Also, I think www.ossim.org will have the utilities you need to do this already built. Look at OrthoIgen.exe. Martin -Original Message- From: gdal-dev-boun...@lists.osgeo.org [mailto:gdal-dev-boun...@lists.osgeo.org] On Behalf Of Rydzak, Carol-P28503 Sent: Thursday, September 03, 2009 8:26 PM To: Even Rouault Cc: gdal-dev@lists.osgeo.org Subject: RE: [gdal-dev] GDAL: How to create one NITF image from CIB CD Hi. Thank you for getting back to me on this! Yes, I am using the 1.6.0 version of GDAL. I did what you suggested and it gives me a bus error when I run the gdalinfo command on the a.toc file. Do you recall what the problems were with the CIB data that you heard about? I have a 1 meter, 5 meter and 10 meter data sets that I tried and I got a bus error on all of them. Any other ideas? Thanks, Carol -Original Message- From: Even Rouault [mailto:even.roua...@mines-paris.org] Sent: Thursday, September 03, 2009 12:50 PM To: Rydzak, Carol-P28503 Cc: gdal-dev@lists.osgeo.org Subject: Re: [gdal-dev] GDAL: How to create one NITF image from CIB CD Selon Rydzak, Carol-P28503 carol.ryd...@gdc4s.com: I assume you are using a recent enough GDAL version (= 1.5.0) that has the RPFTOC driver. If so, gdalinfo on the a.toc should return a list of subdatasets (at least one). You can then select one and translate it to a GeoTIFF. It will mosaic all the underlying NITF tiles. See http://gdal.org/frmt_various.html#RPFTOC for more info. Someone recently reported problems with CIB but I've not heard any news since then. The driver was tested against CADRG datasets but should theoritically work with CIB too. I am trying to use GDAL utilities (on a Solaris 10 system) against CIB files (from a CD) to convert them to a single NITF image. When I run gdal_translate, it doesn't seem to accept the a.toc, but it does create a one-to-one nitf image for each CIB file the command is run against. Is there a GDAL command that takes the table of contents file and processes all the CIB frames into one NITF image? Or takes all the .ntf files that get created and stitches them together into one .ntf image? Thank you, Carol Rydzak General Dynamics C4 Systems Phone: (480) 441-7405 Email: carol.ryd...@gdc4s.com This email message is for the sole use of the intended recipient(s) and may contain GDC4S confidential or privileged information. Any unauthorized review, use, disclosure or distribution is prohibited. If you are not an intended recipient, please contact the sender by reply email and destroy all copies of the original message. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[Gdal-dev] Rotating image...
I have bmp file, and I would like to rotate it on some angle... (For example 38... 45,7... 2,3 degrees). How can I make it by using GDAL? -- View this message in context: http://n2.nabble.com/Rotating-image-tp3577921p3577921.html Sent from the GDAL - Dev mailing list archive at Nabble.com. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev