<div dir="ltr">Hi Folks,<div>I have a directory of 783 .img format Lidar - based  DEMs in the same projection, North Carolina State Plane  Feet (2011) , NAD 83 , NVD88.  I was going to use gdalbuildvrt to create a single virtual image for the area, but gdalbuildvrt complained, informing me that they were not in the same projection.</div><div><br></div><div>A couple of quick bash scripts/commands</div><div><br></div><div>for x in *.img; do gdalinfo $x >>img_info.txt;done<br></div><div><br></div><div>and </div><div>grep PROJCS img_info.txt|sort|uniq -c<br></div><div><br></div><div>gives me the following:</div><div><br></div><div><div>    437 PROJCS["NAD_1983_2011_StatePlane_North_Carolina_FIPS_3200_Ft_US",</div><div>    346 PROJCS["NAD_1983_StatePlane_North_Carolina_FIPS_3200_Feet_2011",</div><div><br></div><div>gdalinfo gives the following for each type of file ( origin section clipped out) :</div><div><br></div><div><div>Driver: HFA/Erdas Imagine Images (.img)</div><div>Files: D05_37_20878102_20141231.img</div><div>       D05_37_20878102_20141231.img.aux.xml</div><div>Size is 1000, 1000</div><div>Coordinate System is:</div><div>PROJCS["NAD_1983_2011_StatePlane_North_Carolina_FIPS_3200_Ft_US",</div><div>    GEOGCS["GCS_NAD_1983_2011",</div><div>        DATUM["NAD_1983_2011",</div><div>            SPHEROID["GRS_1980",6378137.0,298.257222101]],</div><div>        PRIMEM["Greenwich",0.0],</div><div>        UNIT["Degree",0.017453292519943295]],</div><div>    PROJECTION["Lambert_Conformal_Conic_2SP"],</div><div>    PARAMETER["False_Easting",2000000.0],</div><div>    PARAMETER["False_Northing",0.0],</div><div>    PARAMETER["Central_Meridian",-79.0],</div><div>    PARAMETER["Standard_Parallel_1",34.3333333333333],</div><div>    PARAMETER["Standard_Parallel_2",36.1666666666667],</div><div>    PARAMETER["Latitude_Of_Origin",33.75],</div><div>    UNIT["Foot_US",0.30480060960121924],</div><div>    VERTCS["NAVD_1988_Foot_US",</div><div>        VDATUM["North_American_Vertical_Datum_1988"],</div><div>        PARAMETER["Vertical_Shift",0.0],</div><div>        PARAMETER["Direction",1.0],</div><div>        UNIT["Foot_US",0.3048006096012192]]]</div></div><div><br></div><div><br></div><div><div>Driver: HFA/Erdas Imagine Images (.img)</div><div>Files: D05_37_20957301_20141231.img</div><div>       D05_37_20957301_20141231.img.aux.xml</div><div>Size is 1000, 1000</div><div>Coordinate System is:</div><div>PROJCS["NAD_1983_StatePlane_North_Carolina_FIPS_3200_Feet_2011",</div><div>    GEOGCS["GCS_NAD_1983_2011",</div><div>        DATUM["NAD_1983_2011",</div><div>            SPHEROID["GRS_1980",6378137.0,298.257222101]],</div><div>        PRIMEM["Greenwich",0.0],</div><div>        UNIT["Degree",0.0174532925199433]],</div><div>    PROJECTION["Lambert_Conformal_Conic_2SP"],</div><div>    PARAMETER["False_Easting",2000000.002616666],</div><div>    PARAMETER["False_Northing",0.0],</div><div>    PARAMETER["Central_Meridian",-79.0],</div><div>    PARAMETER["Standard_Parallel_1",34.33333333333334],</div><div>    PARAMETER["Standard_Parallel_2",36.16666666666666],</div><div>    PARAMETER["Latitude_Of_Origin",33.75],</div><div>    UNIT["Foot_US",0.30480060960121924],</div><div>    VERTCS["NAVD_1988_Foot_US",</div><div>        VDATUM["North_American_Vertical_Datum_1988"],</div><div>        PARAMETER["Vertical_Shift",0.0],</div><div>        PARAMETER["Direction",1.0],</div><div>        UNIT["Foot_US",0.3048006096012192]]]</div></div><div><br></div><div><br></div><div><br></div><div>In theory, these should all be EPSG:6543, so just assigning the correct  horizontal projection/datum with the epsg code should make things usable. (i.e,</div><div>gdal_translate -a_"srs epsg:6543" --config GDAL_CACHEMAX 512 -of GTiff -co COMPRESS=DEFLATE -co PREDICTOR=3 in.img out.tif )  ( Thank you for the EPSG:6543 projection support in GDAL 2.0!)  <br></div><div><br></div><div>However, this only assigns the horizontal infomation, how does one assign a vertical datum with a horizontal EPSG code? </div><div><br></div><div>I see the NVD88 code in vertcs.csv , but how would I implement it in the above command?</div><div><br></div><div>Using gdal 2.0 Beta2</div><div><br></div><div><br></div><div>Doug</div>-- <br><div class="gmail_signature"><div>Doug Newcomb</div><div>USFWS</div><div>Raleigh, NC</div><div>919-856-4520 ext. 14 <a href="mailto:doug_newcomb@fws.gov" target="_blank">doug_newcomb@fws.gov</a></div><div>---------------------------------------------------------------------------------------------------------</div><div>The opinions I express are my own and are not representative of the official policy of the U.S.Fish and Wildlife Service or Dept. of the Interior.   Life is too short for undocumented, proprietary data formats.</div></div>
</div></div>