Re: [gdal-dev] extract vector/raster data from GeoPDF

2009-09-03 Thread Brent Fraser

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

2009-09-03 Thread Joaquim Luis

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

2009-09-03 Thread Peter J Halls
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

2009-09-03 Thread Joaquim Luis

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

2009-09-03 Thread Jonas Malmqvist

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

2009-09-03 Thread Brent Fraser

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

2009-09-03 Thread jor sion
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

2009-09-03 Thread Frank Warmerdam

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

2009-09-03 Thread Joaquim Luis

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

2009-09-03 Thread Hermann Peifer

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

2009-09-03 Thread 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

2009-09-03 Thread Matt Bartolome
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

2009-09-03 Thread Frank Warmerdam

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

2009-09-03 Thread Daniel Morissette

+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

2009-09-03 Thread Even Rouault
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

2009-09-03 Thread Tamas Szekeres
+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

2009-09-03 Thread Rydzak, Carol-P28503
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

2009-09-03 Thread Martin Hvidberg

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

2009-09-03 Thread Even Rouault
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 ?

2009-09-03 Thread Martin Hvidberg

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

2009-09-03 Thread Joaquim Luis

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

2009-09-03 Thread Martin Chapman
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

2009-09-03 Thread Martin Chapman
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

2009-09-03 Thread Martin Chapman
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...

2009-09-03 Thread ryabis

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