I have been trying to use gdal_polygonize to convert a TIF to an ESRI
Shapefile, but it seems to be losing the projection information (i.e. no
*.prj file is created).
The following output illustrates the loss of projection information.
I'm currently using gdal trunk revision 15512.
Any suggestions on how to keep the projection info?
-------------------------------------------------------
$ gdalinfo testimage.tif
Driver: GTiff/GeoTIFF
Files: testimage.tif
Size is 5233, 1966
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235629972,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-130.538185630190014,39.997456711953717)
Pixel Size = (0.005086568914294,-0.005086568914294)
Metadata:
AREA_OR_POINT=Point
TIFFTAG_SOFTWARE=MODIS Reprojection Tool v4.0 February 2008
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-130.5381856, 39.9974567) (130d32'17.47"W, 39d59'50.84"N)
Lower Left (-130.5381856, 29.9972622) (130d32'17.47"W, 29d59'50.14"N)
Upper Right (-103.9201705, 39.9974567) (103d55'12.61"W, 39d59'50.84"N)
Lower Right (-103.9201705, 29.9972622) (103d55'12.61"W, 29d59'50.14"N)
Center (-117.2291781, 34.9973595) (117d13'45.04"W, 34d59'50.49"N)
Band 1 Block=5233x1 Type=Int16, ColorInterp=Gray
$ /usr/local/gdal/trunk/bin/gdal_polygonize.py testimage.tif -f 'ESRI
Shapefile' test.shp
Creating output test.shp of format ESRI Shapefile.
0...10...20...30...40...50...60...70...80...90...100 - done.
$ ogrinfo test.shp test -so
INFO: Open of `test.shp'
using driver `ESRI Shapefile' successful.
Layer name: test
Geometry: Polygon
Feature Count: 5592
Extent: (-130.538186, 29.997262) - (-103.920171, 39.997457)
Layer SRS WKT:
(unknown)
DN: Real (11.0)
--
Tyler A. Erickson, Ph.D.
Research Scientist, Michigan Tech Research Institute, and
Adjunct Assistant Professor, Dept. of Civil & Environmental Engineering
Michigan Technological University
3600 Green Court, Suite 100
Ann Arbor, MI 48105
W 83.6889°, N 42.3021° (WGS84)
[EMAIL PROTECTED]
(734) 913-6846
www.mtri.org
www.michiganview.org
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
