[gdal-dev] RPC

2016-01-18 Thread Surya Prakash 
Hi ALL,
I have some imagery associated with RPC files. The locational accuracy of the 
images defined by the RPCs is quite poor and I would like to update the RPCs by 
using some additional GCPs. How to do this, Please guide me.

Thanks & regards
SPR___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] RPC

2016-01-18 Thread Even Rouault 
Le lundi 18 janvier 2016 08:39:06, Surya Prakash a écrit :
> Hi ALL,
> I have some imagery associated with RPC files. The locational accuracy of
> the images defined by the RPCs is quite poor and I would like to update
> the RPCs by using some additional GCPs. How to do this, Please guide me.

You might find this interesting :
http://fwarmerdam.blogspot.fr/2013_09_01_archive.html

Are you sure the RPC you got is poor ? Several points need to be taken into 
account to get precise wrapping when using RPCs: if the angle of the image is 
far from nadir, you will need to use a DEM. And one point to take into account 
when doing so is that the elevations of the DEM must be relative to the WGS84 
ellipsoid. For that either you need to correct the DEM before using GDAL, or 
you must use GDAL trunk and have the DEM SRS be a COMPD_CS with a vertical 
datum grid shift. There also were issues with half-pixel shift in RPC 
application that have been fixed very recently in GDAL.

> 
> Thanks & regards
> SPR

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] gdal_proximity: possible bug

2016-01-18 Thread Filipe Dias 
Hi,

I used gdal_proximity on a raster and did not pass the "Values" parameter
because in the help file it said " all non-zero pixels will be considered
target pixels.". The result was an incorrect Distances raster. Afterwards I
set the "Values" parameter and it worked as expected.

I reported this initially in the QGIS bug tracker, see more here:

http://hub.qgis.org/issues/14115#change-67271

I'm using GDAL 1.11.2.

Thanks

Best regards
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] map info and projections from GeoTIFF file to create ENVI header

2016-01-18 Thread Dmitry Baryshnikov 

Hi Nevzat,

you need to get spatial reference from GDALDataset (via GetProjectionRef 
http://gdal.org/classGDALDataset.html#aa42537e1062ce254d124b29ff3ebe857) 
and convert it to OGRSpatialReference.


Than use exportToMICoordSys to get MI CRS string.

http://gdal.org/classOGRSpatialReference.html#a1f2908cd5ca33609844ef0c0ff2186ea

Best regards,
Dmitry

18.01.2016 20:27, ngu...@jlab.org пишет:

It is solved. Thank you very much Even and Dmitry for your helps. I have
missed the GDAL_DATA variable and it caused me all that trouble.

Now I need to get the map info string:
map info = {UTM, 1, 1, 145185, 2189715, 30, 30, 19, North,WGS-84}

Using hDS->GetGeoTransform(adfGeoTransform) does not give me all the
information in the map info string. It only gives me elements [4-7], how to
get the others?

Best regards, - Nevzat


Yes, Thank you Even, I was missing the GDAL_DATA environmental variable. Now I
get the same gdalinfo output. Do I need to set this variable while using the
C++ classes you suggested before: GetGeoTransform and GetProjectionRef()..?
Because, there also I had the same problem.

Best regards, - Nevzat



Le lundi 18 janvier 2016 17:13:14, ngu...@jlab.org a écrit :

Dmitry and All,
I just downloaded the set of tiff files from earthexplorer and run gdalinfo
on it. Please see below. I am missing some fields in the output
information. It is different from what you got. I don't know what can
cause this? My gdalinfo.exe is not good? Any idea?

Have a look at the FAQ entries about GDAL_DATA :
https://trac.osgeo.org/gdal/wiki/FAQInstallationAndBuilding#WhatisGDAL_DATAenvironmentvariable


C:\Cubes\LC80080472013174LGN00.tar>gdalinfo LC80080472013174LGN00_B1.TIF
Driver: GTiff/GeoTIFF
Files: LC80080472013174LGN00_B1.TIF
.\LC80080472013174LGN00_MTL.txt
Size is 7571, 7411
Coordinate System is:
PROJCS["unnamed",
 PROJECTION["Transverse_Mercator"],
 PARAMETER["latitude_of_origin",0],
 PARAMETER["central_meridian",-69],
 PARAMETER["scale_factor",0.9996],
 PARAMETER["false_easting",50],
 PARAMETER["false_northing",0],
 UNIT["metre",1,
 AUTHORITY["EPSG","9001"]]]
Origin = (145185.000,2189715.000)
Pixel Size = (30.000,-30.000)
Metadata:
   AREA_OR_POINT=Point
   METADATATYPE=ODL
Image Structure Metadata:
   INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  145185.000, 2189715.000)
Lower Left  (  145185.000, 1967385.000)
Upper Right (  372315.000, 2189715.000)
Lower Right (  372315.000, 1967385.000)
Center  (  258750.000, 2078550.000)
Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray

C:\Cubes\LC80080472013174LGN00.tar>gdalinfo --version
GDAL 2.0.1, released 2015/09/15

C:\Cubes\LC80080472013174LGN00.tar>


I have the 2.0.1 version:
C:\warmerda\bld\bin>gdalinfo --version
GDAL 2.0.1, released 2015/09/15

I will check the history of the file I have. Maybe it is modified. I will
also download the same file from earthexplorer and try my routines on
it. Thanks for the suggestions. I appreciate.
- Nevzat

-Original Message-
From: Dmitry Baryshnikov [mailto:bishop@gmail.com]
Sent: Sunday, January 17, 2016 12:57 PM
To: Nevzat Guler ; gdal-dev@lists.osgeo.org
Subject: Re: [gdal-dev] map info and projections from GeoTIFF file to
create ENVI header

I downloaded this scene from earthexplorer, and everything is fine:

$ gdalinfo LC80080472013174LGN00_B1.TIF
Driver: GTiff/GeoTIFF
Files: LC80080472013174LGN00_B1.TIF

 ./LC80080472013174LGN00_MTL.txt

Size is 7571, 7411
Coordinate System is:
PROJCS["WGS 84 / UTM zone 19N",

  GEOGCS["WGS 84",

  DATUM["WGS_1984",

  SPHEROID["WGS 84",6378137,298.257223563,

  AUTHORITY["EPSG","7030"]],

  AUTHORITY["EPSG","6326"]],

  PRIMEM["Greenwich",0,

  AUTHORITY["EPSG","8901"]],

  UNIT["degree",0.0174532925199433,

  AUTHORITY["EPSG","9122"]],

  AUTHORITY["EPSG","4326"]],

  PROJECTION["Transverse_Mercator"],
  PARAMETER["latitude_of_origin",0],
  PARAMETER["central_meridian",-69],
  PARAMETER["scale_factor",0.9996],
  PARAMETER["false_easting",50],
  PARAMETER["false_northing",0],
  UNIT["metre",1,

  AUTHORITY["EPSG","9001"]],

  AXIS["Easting",EAST],
  AXIS["Northing",NORTH],
  AUTHORITY["EPSG","32619"]]

Origin = (145185.000,2189715.000)
Pixel Size = (30.000,-30.000)

Metadata:
AREA_OR_POINT=Point
METADATATYPE=ODL

Image Structure Metadata:
INTERLEAVE=BAND

Corner Coordinates:
Upper Left  (  145185.000, 2189715.000) ( 72d23' 7.97"W, 19d46'16.43"N)
Lower Left  (  145185.000, 1967385.000) ( 72d20'44.51"W, 17d45'55.15"N)
Upper Right (  372315.000, 2189715.000) ( 70d13' 8.50"W, 19d47'56.92"N)
Lower Right ( 372315.000, 1967385.000) ( 70d12'16.74"W, 17d47'24.74"N)
Center  (  258750.000, 2078550.000) ( 71d17'19.50"W, 18d47' 4.68"N)
Band

Re: [gdal-dev] map info and projections from GeoTIFF file to create ENVI header

2016-01-18 Thread nguler 
It is solved. Thank you very much Even and Dmitry for your helps. I have
missed the GDAL_DATA variable and it caused me all that trouble.

Now I need to get the map info string:
map info = {UTM, 1, 1, 145185, 2189715, 30, 30, 19, North,WGS-84}

Using hDS->GetGeoTransform(adfGeoTransform) does not give me all the
information in the map info string. It only gives me elements [4-7], how to
get the others?

Best regards, - Nevzat

> Yes, Thank you Even, I was missing the GDAL_DATA environmental variable. Now I
> get the same gdalinfo output. Do I need to set this variable while using the
> C++ classes you suggested before: GetGeoTransform and GetProjectionRef()..?
> Because, there also I had the same problem.
>
> Best regards, - Nevzat
>
>
>> Le lundi 18 janvier 2016 17:13:14, ngu...@jlab.org a écrit :
>>> Dmitry and All,
>>> I just downloaded the set of tiff files from earthexplorer and run gdalinfo
>>> on it. Please see below. I am missing some fields in the output
>>> information. It is different from what you got. I don't know what can
>>> cause this? My gdalinfo.exe is not good? Any idea?
>>
>> Have a look at the FAQ entries about GDAL_DATA :
>> https://trac.osgeo.org/gdal/wiki/FAQInstallationAndBuilding#WhatisGDAL_DATAenvironmentvariable
>>
>>>
>>> C:\Cubes\LC80080472013174LGN00.tar>gdalinfo LC80080472013174LGN00_B1.TIF
>>> Driver: GTiff/GeoTIFF
>>> Files: LC80080472013174LGN00_B1.TIF
>>>.\LC80080472013174LGN00_MTL.txt
>>> Size is 7571, 7411
>>> Coordinate System is:
>>> PROJCS["unnamed",
>>> PROJECTION["Transverse_Mercator"],
>>> PARAMETER["latitude_of_origin",0],
>>> PARAMETER["central_meridian",-69],
>>> PARAMETER["scale_factor",0.9996],
>>> PARAMETER["false_easting",50],
>>> PARAMETER["false_northing",0],
>>> UNIT["metre",1,
>>> AUTHORITY["EPSG","9001"]]]
>>> Origin = (145185.000,2189715.000)
>>> Pixel Size = (30.000,-30.000)
>>> Metadata:
>>>   AREA_OR_POINT=Point
>>>   METADATATYPE=ODL
>>> Image Structure Metadata:
>>>   INTERLEAVE=BAND
>>> Corner Coordinates:
>>> Upper Left  (  145185.000, 2189715.000)
>>> Lower Left  (  145185.000, 1967385.000)
>>> Upper Right (  372315.000, 2189715.000)
>>> Lower Right (  372315.000, 1967385.000)
>>> Center  (  258750.000, 2078550.000)
>>> Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray
>>>
>>> C:\Cubes\LC80080472013174LGN00.tar>gdalinfo --version
>>> GDAL 2.0.1, released 2015/09/15
>>>
>>> C:\Cubes\LC80080472013174LGN00.tar>
>>>
>>> > I have the 2.0.1 version:
>>> > C:\warmerda\bld\bin>gdalinfo --version
>>> > GDAL 2.0.1, released 2015/09/15
>>> >
>>> > I will check the history of the file I have. Maybe it is modified. I will
>>> > also download the same file from earthexplorer and try my routines on
>>> > it. Thanks for the suggestions. I appreciate.
>>> > - Nevzat
>>> >
>>> > -Original Message-
>>> > From: Dmitry Baryshnikov [mailto:bishop@gmail.com]
>>> > Sent: Sunday, January 17, 2016 12:57 PM
>>> > To: Nevzat Guler ; gdal-dev@lists.osgeo.org
>>> > Subject: Re: [gdal-dev] map info and projections from GeoTIFF file to
>>> > create ENVI header
>>> >
>>> > I downloaded this scene from earthexplorer, and everything is fine:
>>> >
>>> > $ gdalinfo LC80080472013174LGN00_B1.TIF
>>> > Driver: GTiff/GeoTIFF
>>> > Files: LC80080472013174LGN00_B1.TIF
>>> >
>>> > ./LC80080472013174LGN00_MTL.txt
>>> >
>>> > Size is 7571, 7411
>>> > Coordinate System is:
>>> > PROJCS["WGS 84 / UTM zone 19N",
>>> >
>>> >  GEOGCS["WGS 84",
>>> >
>>> >  DATUM["WGS_1984",
>>> >
>>> >  SPHEROID["WGS 84",6378137,298.257223563,
>>> >
>>> >  AUTHORITY["EPSG","7030"]],
>>> >
>>> >  AUTHORITY["EPSG","6326"]],
>>> >
>>> >  PRIMEM["Greenwich",0,
>>> >
>>> >  AUTHORITY["EPSG","8901"]],
>>> >
>>> >  UNIT["degree",0.0174532925199433,
>>> >
>>> >  AUTHORITY["EPSG","9122"]],
>>> >
>>> >  AUTHORITY["EPSG","4326"]],
>>> >
>>> >  PROJECTION["Transverse_Mercator"],
>>> >  PARAMETER["latitude_of_origin",0],
>>> >  PARAMETER["central_meridian",-69],
>>> >  PARAMETER["scale_factor",0.9996],
>>> >  PARAMETER["false_easting",50],
>>> >  PARAMETER["false_northing",0],
>>> >  UNIT["metre",1,
>>> >
>>> >  AUTHORITY["EPSG","9001"]],
>>> >
>>> >  AXIS["Easting",EAST],
>>> >  AXIS["Northing",NORTH],
>>> >  AUTHORITY["EPSG","32619"]]
>>> >
>>> > Origin = (145185.000,2189715.000)
>>> > Pixel Size = (30.000,-30.000)
>>> >
>>> > Metadata:
>>> >AREA_OR_POINT=Point
>>> >METADATATYPE=ODL
>>> >
>>> > Image Structure Metadata:
>>> >INTERLEAVE=BAND
>>> >
>>> > Corner Coordinates:
>>> > Upper Left  (  145185.000, 2189715.000) ( 72d23' 7.97"W, 19d46'16.43"N)
>>> > Lower Left  (  145185.000, 1967385.000) ( 72d20'44.51"W, 17d45'55.15"N)
>>> > Upper Right (  372315.000, 2189715.000) (

[gdal-dev] GDAL license and copyright

2016-01-18 Thread Aaron Boxer 
Hello,

I am trying to learn more about how open source projects manage copyright
of contributors.
I see that various GDAL files have different copyright assignees, although
they all use the MIT
license.

Is this difficult to manage? For example, if GDAL is packaged for a
distribution such as Fedora,
how do you manage the license for the package?  Also, if a new contributor
were to add new functionality in a new file, is it ok for them to assign
copyright to themselves when submitting to the trunk ?

Thanks very much,
Aaron
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] Intersection of two Layers

2016-01-18 Thread Kai Muehlbauer 

Hi,

at the moment I calculate the intersection of two OGR Layers (A, B) by 
iterating over features/geometries (GA) of layer A and then calling 
GA.Intersection(GB) for every Feature/Geometry (GB) of B. I set a 
spatial filter for layer B with every GA. The intersections 
features/geometries are written to a result layer (C).


In my use case the features/geometries of layer B are much smaller than 
the features of layer A. Means that quite a big number of B's geometries 
are fully contained in A's geometries.


To speed things up, I check for containment GA.contains(GB). If GB is 
contained in GA, I just copy GB to C, otherwise I calculate 
GA.Intersection(GB) and write this to C.


I found that I also could use the layer intersection function 
(OGRLayer::Intersection). Unfortunately the intersection is calculated 
for every feature/geometry.


Would it be possible to extend the layer intersection function to allow 
for containment check (eg. as an option) and just copy the contained 
geometries without calculation of expensive intersection? Suggestions, 
opinions?


In my use case I get improvements of around 20% of processing time.

Cheers,
Kai
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] map info and projections from GeoTIFF file to create ENVI header

2016-01-18 Thread nguler 
Dmitry and All,
I just downloaded the set of tiff files from earthexplorer and run gdalinfo on
it. Please see below. I am missing some fields in the output information. It
is different from what you got. I don't know what can cause this? My
gdalinfo.exe is not good? Any idea?

C:\Cubes\LC80080472013174LGN00.tar>gdalinfo LC80080472013174LGN00_B1.TIF
Driver: GTiff/GeoTIFF
Files: LC80080472013174LGN00_B1.TIF
   .\LC80080472013174LGN00_MTL.txt
Size is 7571, 7411
Coordinate System is:
PROJCS["unnamed",
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-69],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",50],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (145185.000,2189715.000)
Pixel Size = (30.000,-30.000)
Metadata:
  AREA_OR_POINT=Point
  METADATATYPE=ODL
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  145185.000, 2189715.000)
Lower Left  (  145185.000, 1967385.000)
Upper Right (  372315.000, 2189715.000)
Lower Right (  372315.000, 1967385.000)
Center  (  258750.000, 2078550.000)
Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray

C:\Cubes\LC80080472013174LGN00.tar>gdalinfo --version
GDAL 2.0.1, released 2015/09/15

C:\Cubes\LC80080472013174LGN00.tar>


> I have the 2.0.1 version:
> C:\warmerda\bld\bin>gdalinfo --version
> GDAL 2.0.1, released 2015/09/15
>
> I will check the history of the file I have. Maybe it is modified. I will also
> download the same file from earthexplorer and try my routines on it.
> Thanks for the suggestions. I appreciate.
> - Nevzat
>
> -Original Message-
> From: Dmitry Baryshnikov [mailto:bishop@gmail.com]
> Sent: Sunday, January 17, 2016 12:57 PM
> To: Nevzat Guler ; gdal-dev@lists.osgeo.org
> Subject: Re: [gdal-dev] map info and projections from GeoTIFF file to create
> ENVI header
>
> I downloaded this scene from earthexplorer, and everything is fine:
>
> $ gdalinfo LC80080472013174LGN00_B1.TIF
> Driver: GTiff/GeoTIFF
> Files: LC80080472013174LGN00_B1.TIF
> ./LC80080472013174LGN00_MTL.txt
> Size is 7571, 7411
> Coordinate System is:
> PROJCS["WGS 84 / UTM zone 19N",
>  GEOGCS["WGS 84",
>  DATUM["WGS_1984",
>  SPHEROID["WGS 84",6378137,298.257223563,
>  AUTHORITY["EPSG","7030"]],
>  AUTHORITY["EPSG","6326"]],
>  PRIMEM["Greenwich",0,
>  AUTHORITY["EPSG","8901"]],
>  UNIT["degree",0.0174532925199433,
>  AUTHORITY["EPSG","9122"]],
>  AUTHORITY["EPSG","4326"]],
>  PROJECTION["Transverse_Mercator"],
>  PARAMETER["latitude_of_origin",0],
>  PARAMETER["central_meridian",-69],
>  PARAMETER["scale_factor",0.9996],
>  PARAMETER["false_easting",50],
>  PARAMETER["false_northing",0],
>  UNIT["metre",1,
>  AUTHORITY["EPSG","9001"]],
>  AXIS["Easting",EAST],
>  AXIS["Northing",NORTH],
>  AUTHORITY["EPSG","32619"]]
> Origin = (145185.000,2189715.000)
> Pixel Size = (30.000,-30.000)
> Metadata:
>AREA_OR_POINT=Point
>METADATATYPE=ODL
> Image Structure Metadata:
>INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (  145185.000, 2189715.000) ( 72d23' 7.97"W, 19d46'16.43"N) Lower
> Left  (  145185.000, 1967385.000) ( 72d20'44.51"W, 17d45'55.15"N) Upper Right
> (  372315.000, 2189715.000) ( 70d13' 8.50"W, 19d47'56.92"N) Lower Right (
> 372315.000, 1967385.000) ( 70d12'16.74"W, 17d47'24.74"N)
> Center  (  258750.000, 2078550.000) ( 71d17'19.50"W, 18d47' 4.68"N)
> Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray
>
> $ gdalinfo --version
> GDAL 2.1.0dev, released 2015/99/99
>
> What is your gdal version? Maybe the Landsat geotiff was modified from some
> software?
>
> Best regards,
>  Dmitry
>
> 17.01.2016 20:15, Nevzat Guler пишет:
>> Thank you Dmitry,
>> What is the best way to extract the information we got using gdalinfo.exe by
>> using GDAL classes, within a C++ code? As I mentioned, I am using
>> GetGeoTransform and GetProjectionRef() but the results miss many of the
>> information (specifically the map info elements as I explained in the
>> previous email).
>>
>> On another note, I actually get a different output from gdalinfo.exe using
>> my file (please see the printed result below). I got the geotiff file from
>> my colleague. When I open the file with ENVI, I get all the information,
>> including GEOGCS information.
>>
>> Can it be that my file is just corrupted and that's why I don't get all the
>> information from GetGeoTransform and GetProjectionRef()? But then how would
>> ENVI is able to extract this information if it were not in the file?
>>
>> Thanks again, - Nevzat
>>
>> C:\Cubes\Landsat\LC80080472013174LGN00>gdalinfo
>> LC80080472013174LGN00_B1.TIF
>> Driver: GTiff/GeoTIFF
>> Files:

Re: [gdal-dev] map info and projections from GeoTIFF file to create ENVI header

2016-01-18 Thread Even Rouault 
Le lundi 18 janvier 2016 17:13:14, ngu...@jlab.org a écrit :
> Dmitry and All,
> I just downloaded the set of tiff files from earthexplorer and run gdalinfo
> on it. Please see below. I am missing some fields in the output
> information. It is different from what you got. I don't know what can
> cause this? My gdalinfo.exe is not good? Any idea?

Have a look at the FAQ entries about GDAL_DATA :
https://trac.osgeo.org/gdal/wiki/FAQInstallationAndBuilding#WhatisGDAL_DATAenvironmentvariable

> 
> C:\Cubes\LC80080472013174LGN00.tar>gdalinfo LC80080472013174LGN00_B1.TIF
> Driver: GTiff/GeoTIFF
> Files: LC80080472013174LGN00_B1.TIF
>.\LC80080472013174LGN00_MTL.txt
> Size is 7571, 7411
> Coordinate System is:
> PROJCS["unnamed",
> PROJECTION["Transverse_Mercator"],
> PARAMETER["latitude_of_origin",0],
> PARAMETER["central_meridian",-69],
> PARAMETER["scale_factor",0.9996],
> PARAMETER["false_easting",50],
> PARAMETER["false_northing",0],
> UNIT["metre",1,
> AUTHORITY["EPSG","9001"]]]
> Origin = (145185.000,2189715.000)
> Pixel Size = (30.000,-30.000)
> Metadata:
>   AREA_OR_POINT=Point
>   METADATATYPE=ODL
> Image Structure Metadata:
>   INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (  145185.000, 2189715.000)
> Lower Left  (  145185.000, 1967385.000)
> Upper Right (  372315.000, 2189715.000)
> Lower Right (  372315.000, 1967385.000)
> Center  (  258750.000, 2078550.000)
> Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray
> 
> C:\Cubes\LC80080472013174LGN00.tar>gdalinfo --version
> GDAL 2.0.1, released 2015/09/15
> 
> C:\Cubes\LC80080472013174LGN00.tar>
> 
> > I have the 2.0.1 version:
> > C:\warmerda\bld\bin>gdalinfo --version
> > GDAL 2.0.1, released 2015/09/15
> > 
> > I will check the history of the file I have. Maybe it is modified. I will
> > also download the same file from earthexplorer and try my routines on
> > it. Thanks for the suggestions. I appreciate.
> > - Nevzat
> > 
> > -Original Message-
> > From: Dmitry Baryshnikov [mailto:bishop@gmail.com]
> > Sent: Sunday, January 17, 2016 12:57 PM
> > To: Nevzat Guler ; gdal-dev@lists.osgeo.org
> > Subject: Re: [gdal-dev] map info and projections from GeoTIFF file to
> > create ENVI header
> > 
> > I downloaded this scene from earthexplorer, and everything is fine:
> > 
> > $ gdalinfo LC80080472013174LGN00_B1.TIF
> > Driver: GTiff/GeoTIFF
> > Files: LC80080472013174LGN00_B1.TIF
> > 
> > ./LC80080472013174LGN00_MTL.txt
> > 
> > Size is 7571, 7411
> > Coordinate System is:
> > PROJCS["WGS 84 / UTM zone 19N",
> > 
> >  GEOGCS["WGS 84",
> >  
> >  DATUM["WGS_1984",
> >  
> >  SPHEROID["WGS 84",6378137,298.257223563,
> >  
> >  AUTHORITY["EPSG","7030"]],
> >  
> >  AUTHORITY["EPSG","6326"]],
> >  
> >  PRIMEM["Greenwich",0,
> >  
> >  AUTHORITY["EPSG","8901"]],
> >  
> >  UNIT["degree",0.0174532925199433,
> >  
> >  AUTHORITY["EPSG","9122"]],
> >  
> >  AUTHORITY["EPSG","4326"]],
> >  
> >  PROJECTION["Transverse_Mercator"],
> >  PARAMETER["latitude_of_origin",0],
> >  PARAMETER["central_meridian",-69],
> >  PARAMETER["scale_factor",0.9996],
> >  PARAMETER["false_easting",50],
> >  PARAMETER["false_northing",0],
> >  UNIT["metre",1,
> >  
> >  AUTHORITY["EPSG","9001"]],
> >  
> >  AXIS["Easting",EAST],
> >  AXIS["Northing",NORTH],
> >  AUTHORITY["EPSG","32619"]]
> > 
> > Origin = (145185.000,2189715.000)
> > Pixel Size = (30.000,-30.000)
> > 
> > Metadata:
> >AREA_OR_POINT=Point
> >METADATATYPE=ODL
> > 
> > Image Structure Metadata:
> >INTERLEAVE=BAND
> > 
> > Corner Coordinates:
> > Upper Left  (  145185.000, 2189715.000) ( 72d23' 7.97"W, 19d46'16.43"N)
> > Lower Left  (  145185.000, 1967385.000) ( 72d20'44.51"W, 17d45'55.15"N)
> > Upper Right (  372315.000, 2189715.000) ( 70d13' 8.50"W, 19d47'56.92"N)
> > Lower Right ( 372315.000, 1967385.000) ( 70d12'16.74"W, 17d47'24.74"N)
> > Center  (  258750.000, 2078550.000) ( 71d17'19.50"W, 18d47' 4.68"N)
> > Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray
> > 
> > $ gdalinfo --version
> > GDAL 2.1.0dev, released 2015/99/99
> > 
> > What is your gdal version? Maybe the Landsat geotiff was modified from
> > some software?
> > 
> > Best regards,
> > 
> >  Dmitry
> > 
> > 17.01.2016 20:15, Nevzat Guler пишет:
> >> Thank you Dmitry,
> >> What is the best way to extract the information we got using
> >> gdalinfo.exe by using GDAL classes, within a C++ code? As I mentioned,
> >> I am using GetGeoTransform and GetProjectionRef() but the results miss
> >> many of the information (specifically the map info elements as I
> >> explained in the previous email).
> >> 
> >> On

Re: [gdal-dev] map info and projections from GeoTIFF file to create ENVI header

2016-01-18 Thread nguler 
Yes, Thank you Even, I was missing the GDAL_DATA environmental variable. Now I
get the same gdalinfo output. Do I need to set this variable while using the
C++ classes you suggested before: GetGeoTransform and GetProjectionRef()..?
Because, there also I had the same problem.

Best regards, - Nevzat


> Le lundi 18 janvier 2016 17:13:14, ngu...@jlab.org a écrit :
>> Dmitry and All,
>> I just downloaded the set of tiff files from earthexplorer and run gdalinfo
>> on it. Please see below. I am missing some fields in the output
>> information. It is different from what you got. I don't know what can
>> cause this? My gdalinfo.exe is not good? Any idea?
>
> Have a look at the FAQ entries about GDAL_DATA :
> https://trac.osgeo.org/gdal/wiki/FAQInstallationAndBuilding#WhatisGDAL_DATAenvironmentvariable
>
>>
>> C:\Cubes\LC80080472013174LGN00.tar>gdalinfo LC80080472013174LGN00_B1.TIF
>> Driver: GTiff/GeoTIFF
>> Files: LC80080472013174LGN00_B1.TIF
>>.\LC80080472013174LGN00_MTL.txt
>> Size is 7571, 7411
>> Coordinate System is:
>> PROJCS["unnamed",
>> PROJECTION["Transverse_Mercator"],
>> PARAMETER["latitude_of_origin",0],
>> PARAMETER["central_meridian",-69],
>> PARAMETER["scale_factor",0.9996],
>> PARAMETER["false_easting",50],
>> PARAMETER["false_northing",0],
>> UNIT["metre",1,
>> AUTHORITY["EPSG","9001"]]]
>> Origin = (145185.000,2189715.000)
>> Pixel Size = (30.000,-30.000)
>> Metadata:
>>   AREA_OR_POINT=Point
>>   METADATATYPE=ODL
>> Image Structure Metadata:
>>   INTERLEAVE=BAND
>> Corner Coordinates:
>> Upper Left  (  145185.000, 2189715.000)
>> Lower Left  (  145185.000, 1967385.000)
>> Upper Right (  372315.000, 2189715.000)
>> Lower Right (  372315.000, 1967385.000)
>> Center  (  258750.000, 2078550.000)
>> Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray
>>
>> C:\Cubes\LC80080472013174LGN00.tar>gdalinfo --version
>> GDAL 2.0.1, released 2015/09/15
>>
>> C:\Cubes\LC80080472013174LGN00.tar>
>>
>> > I have the 2.0.1 version:
>> > C:\warmerda\bld\bin>gdalinfo --version
>> > GDAL 2.0.1, released 2015/09/15
>> >
>> > I will check the history of the file I have. Maybe it is modified. I will
>> > also download the same file from earthexplorer and try my routines on
>> > it. Thanks for the suggestions. I appreciate.
>> > - Nevzat
>> >
>> > -Original Message-
>> > From: Dmitry Baryshnikov [mailto:bishop@gmail.com]
>> > Sent: Sunday, January 17, 2016 12:57 PM
>> > To: Nevzat Guler ; gdal-dev@lists.osgeo.org
>> > Subject: Re: [gdal-dev] map info and projections from GeoTIFF file to
>> > create ENVI header
>> >
>> > I downloaded this scene from earthexplorer, and everything is fine:
>> >
>> > $ gdalinfo LC80080472013174LGN00_B1.TIF
>> > Driver: GTiff/GeoTIFF
>> > Files: LC80080472013174LGN00_B1.TIF
>> >
>> > ./LC80080472013174LGN00_MTL.txt
>> >
>> > Size is 7571, 7411
>> > Coordinate System is:
>> > PROJCS["WGS 84 / UTM zone 19N",
>> >
>> >  GEOGCS["WGS 84",
>> >
>> >  DATUM["WGS_1984",
>> >
>> >  SPHEROID["WGS 84",6378137,298.257223563,
>> >
>> >  AUTHORITY["EPSG","7030"]],
>> >
>> >  AUTHORITY["EPSG","6326"]],
>> >
>> >  PRIMEM["Greenwich",0,
>> >
>> >  AUTHORITY["EPSG","8901"]],
>> >
>> >  UNIT["degree",0.0174532925199433,
>> >
>> >  AUTHORITY["EPSG","9122"]],
>> >
>> >  AUTHORITY["EPSG","4326"]],
>> >
>> >  PROJECTION["Transverse_Mercator"],
>> >  PARAMETER["latitude_of_origin",0],
>> >  PARAMETER["central_meridian",-69],
>> >  PARAMETER["scale_factor",0.9996],
>> >  PARAMETER["false_easting",50],
>> >  PARAMETER["false_northing",0],
>> >  UNIT["metre",1,
>> >
>> >  AUTHORITY["EPSG","9001"]],
>> >
>> >  AXIS["Easting",EAST],
>> >  AXIS["Northing",NORTH],
>> >  AUTHORITY["EPSG","32619"]]
>> >
>> > Origin = (145185.000,2189715.000)
>> > Pixel Size = (30.000,-30.000)
>> >
>> > Metadata:
>> >AREA_OR_POINT=Point
>> >METADATATYPE=ODL
>> >
>> > Image Structure Metadata:
>> >INTERLEAVE=BAND
>> >
>> > Corner Coordinates:
>> > Upper Left  (  145185.000, 2189715.000) ( 72d23' 7.97"W, 19d46'16.43"N)
>> > Lower Left  (  145185.000, 1967385.000) ( 72d20'44.51"W, 17d45'55.15"N)
>> > Upper Right (  372315.000, 2189715.000) ( 70d13' 8.50"W, 19d47'56.92"N)
>> > Lower Right ( 372315.000, 1967385.000) ( 70d12'16.74"W, 17d47'24.74"N)
>> > Center  (  258750.000, 2078550.000) ( 71d17'19.50"W, 18d47' 4.68"N)
>> > Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray
>> >
>> > $ gdalinfo --version
>> > GDAL 2.1.0dev, released 2015/99/99
>> >
>> > What is your gdal version? Maybe the Landsat geotiff was modified from
>> > some software?
>> >
>> > Best regards,
>> >
>> >  Dmitry
>> >
>> > 17.01.2016 20:15, Nevzat Guler пишет:
>> >> Thank you Dmitry,
>> >> What is the best

Re: [gdal-dev] GDAL license and copyright

2016-01-18 Thread Even Rouault 
Le lundi 18 janvier 2016 19:59:01, Aaron Boxer a écrit :
> Hello,
> 
> I am trying to learn more about how open source projects manage copyright
> of contributors.
> I see that various GDAL files have different copyright assignees, although
> they all use the MIT
> license.
> 
> Is this difficult to manage? 

Not until now. Copyright assignment to a single entity can be useful when 
dealing with potential future license changes, but given that X/MIT is one of 
the less restrictive license, it is unlikely to be needed to be changed.

> For example, if GDAL is packaged for a
> distribution such as Fedora,
> how do you manage the license for the package? 

I think distributions refer to the LICENSE.TXT file at the root of GDAL source 
tree. The global licensing implications can however by more complex if you 
link GDAL with third-party libraries that can be GPL, LGPL, proprietary, 
etc...

> Also, if a new contributor
> were to add new functionality in a new file, is it ok for them to assign
> copyright to themselves when submitting to the trunk ?

Yes.

> 
> Thanks very much,
> Aaron

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] GDAL license and copyright

2016-01-18 Thread Aaron Boxer 
On Mon, Jan 18, 2016 at 2:22 PM, Even Rouault 
wrote:

> Le lundi 18 janvier 2016 19:59:01, Aaron Boxer a écrit :
> > Hello,
> >
> > I am trying to learn more about how open source projects manage copyright
> > of contributors.
> > I see that various GDAL files have different copyright assignees,
> although
> > they all use the MIT
> > license.
> >
> > Is this difficult to manage?
>
> Not until now. Copyright assignment to a single entity can be useful when
> dealing with potential future license changes, but given that X/MIT is one
> of
> the less restrictive license, it is unlikely to be needed to be changed.
>
> > For example, if GDAL is packaged for a
> > distribution such as Fedora,
> > how do you manage the license for the package?
>
> I think distributions refer to the LICENSE.TXT file at the root of GDAL
> source
> tree. The global licensing implications can however by more complex if you
> link GDAL with third-party libraries that can be GPL, LGPL, proprietary,
> etc...
>
> > Also, if a new contributor
> > were to add new functionality in a new file, is it ok for them to assign
> > copyright to themselves when submitting to the trunk ?
>
> Yes.
>
>
Awesome, thanks!
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] GDAL license and copyright

2016-01-18 Thread Greg Troxel 

Aaron Boxer  writes:

> I am trying to learn more about how open source projects manage
> copyright of contributors.  I see that various GDAL files have
> different copyright assignees, although they all use the MIT license.

(I'm not trying to counter anything Evan said; he is not only an
authoritative speaker but IMHO said everythign totally right.  I am the
chief licensing geek in pkgsrc, a multi-OS multi-arch packaging system.
This is more than what you asked, but hopefully useful background.)

As always, IANAL, TINLA.  And this is a US-centric view, although I
think not so far off for Europe and many other Berne Convention
countries.

You said "different assignees".  I think you meant "different copyright
holders".  When a person creates a work (writes code, in nerd-speak),
the basic consequence is that they hold copyright, and there is no
assignment involved.  If the code is written within the scope of
employment, then the copyright will be held by the employer, under the
"work for hire" doctrine.

Beyond that ownership at creation, copyright can be assigned by a
contract from the original holder to a new holder. (As an exception,
Europe has the concept of "moral rights" as an element of copyright
which cannot be assigned, but this has so far not been a significant
issue in software.)

Some open source projects ask that copyright be assigned to some entity
(e.g., emacs asks for an assignment to FSF, and will not merge code
until that is done).  Others have a very large number of copyright
holders.

As a packager, I don't care how many copyright holders there are.  What
I look for is a clear license, since without a license, there is no
permission to copy or to create derivative works.  The common practice
of "inbound = outbound", where contributions are only accepted if there
is a license grant equal to the project's outgoing license, both works
well and fits very well with most people's senses of fairness and
reciprocity.

As a contributor, I am happy when I get to keep copyright, not because I
really want to keep it, but because it is easy and does not require
signing something, which usually has additional terms I may not like
(e.g., indemnification).  I do not object conceptually to assigning to a
501(c)3 organization that has promotion of Free Software as a
fundamental goal (e.g. FSF, SFC, ASF).  I am not comfortable assiging to
a company, and object unless the assignment guarantees that the code
will *only* be distributed by them under Free licenses.

Greg, lurker



signature.asc
Description: PGP signature
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] OpenJPEG tiled encoding

2016-01-18 Thread Aaron Boxer 
On Mon, Jan 18, 2016 at 5:14 PM, Even Rouault 
wrote:

> Le lundi 18 janvier 2016 22:41:38, Aaron Boxer a écrit :
> > Hello,
> >
> > It looks like OpenJPEG tiled encoding is broken for non power of 2
> images.
> > This seems to be quite an old bug.
>
> I don't reproduce that. For example, the following works :
>
> gdal_translate ../autotest/gdrivers/data/utm.tif out.jp2 -of jp2openjpeg
> -co
> blockxsize=33 -co blockysize=35 -outsize 100 100
>



You're right. This is actually a problem with kdu_show, not OpenJPEG.  I
was assuming
that Kakadu was the gold standard, but it has bugs too :)






>
> >
> > Does this type of encoding get used in GDAL ? Or are all images power of
> 2
> > dimensions.
> >
> > Thanks,
> > Aaron
>
> --
> Spatialys - Geospatial professional services
> http://www.spatialys.com
>
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Intersection of two Layers

2016-01-18 Thread Ari Jolma 

18.01.2016, 11:56, Kai Muehlbauer kirjoitti:

Hi,

at the moment I calculate the intersection of two OGR Layers (A, B) by 
iterating over features/geometries (GA) of layer A and then calling 
GA.Intersection(GB) for every Feature/Geometry (GB) of B. I set a 
spatial filter for layer B with every GA. The intersections 
features/geometries are written to a result layer (C).


In my use case the features/geometries of layer B are much smaller 
than the features of layer A. Means that quite a big number of B's 
geometries are fully contained in A's geometries.


To speed things up, I check for containment GA.contains(GB). If GB is 
contained in GA, I just copy GB to C, otherwise I calculate 
GA.Intersection(GB) and write this to C.


I found that I also could use the layer intersection function 
(OGRLayer::Intersection). Unfortunately the intersection is calculated 
for every feature/geometry.


Would it be possible to extend the layer intersection function to 
allow for containment check (eg. as an option) and just copy the 
contained geometries without calculation of expensive intersection? 
Suggestions, opinions?


There's the Options argument, which could be used. Something like 
"PRETEST=CONTAINMENT" could indicate a wish to test containment before 
intersection.


A test dataset and a patch would be useful.

Ari



In my use case I get improvements of around 20% of processing time.

Cheers,
Kai
___
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] OpenJPEG tiled encoding

2016-01-18 Thread Even Rouault 
Le lundi 18 janvier 2016 22:41:38, Aaron Boxer a écrit :
> Hello,
> 
> It looks like OpenJPEG tiled encoding is broken for non power of 2 images.
> This seems to be quite an old bug.

I don't reproduce that. For example, the following works :

gdal_translate ../autotest/gdrivers/data/utm.tif out.jp2 -of jp2openjpeg -co 
blockxsize=33 -co blockysize=35 -outsize 100 100

> 
> Does this type of encoding get used in GDAL ? Or are all images power of 2
> dimensions.
> 
> Thanks,
> Aaron

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] map info and projections from GeoTIFF file to create ENVI header

2016-01-18 Thread nguler 
Dmitry,

After doing this, I following your suggested steps, I get the string:

Earth Projection 8, 104, "m", -69, 0, 0.9996, 50, 0

Bit, the map info string I need to use looks like this:

map info = {UTM, 1, 1, 145185, 2189715, 30, 30, 19, North}

I tried other ExportTo functions but none gives me the specific map info
string I want. This is ENVI map info string.

I guess I will just have to stitch it together once I understand what the
second and third elements are used for..

I really appreciated your help.
Best regards, - Nevzat


> Hi Nevzat,
>
> you need to get spatial reference from GDALDataset (via GetProjectionRef
> http://gdal.org/classGDALDataset.html#aa42537e1062ce254d124b29ff3ebe857)
> and convert it to OGRSpatialReference.
>
> Than use exportToMICoordSys to get MI CRS string.
>
> http://gdal.org/classOGRSpatialReference.html#a1f2908cd5ca33609844ef0c0ff2186ea
>
> Best regards,
>  Dmitry
>
> 18.01.2016 20:27, ngu...@jlab.org пишет:
>> It is solved. Thank you very much Even and Dmitry for your helps. I have
>> missed the GDAL_DATA variable and it caused me all that trouble.
>>
>> Now I need to get the map info string:
>> map info = {UTM, 1, 1, 145185, 2189715, 30, 30, 19, North,WGS-84}
>>
>> Using hDS->GetGeoTransform(adfGeoTransform) does not give me all the
>> information in the map info string. It only gives me elements [4-7], how to
>> get the others?
>>
>> Best regards, - Nevzat
>>
>>> Yes, Thank you Even, I was missing the GDAL_DATA environmental variable.
>>> Now I
>>> get the same gdalinfo output. Do I need to set this variable while using
>>> the
>>> C++ classes you suggested before: GetGeoTransform and GetProjectionRef()..?
>>> Because, there also I had the same problem.
>>>
>>> Best regards, - Nevzat
>>>
>>>
 Le lundi 18 janvier 2016 17:13:14, ngu...@jlab.org a écrit :
> Dmitry and All,
> I just downloaded the set of tiff files from earthexplorer and run
> gdalinfo
> on it. Please see below. I am missing some fields in the output
> information. It is different from what you got. I don't know what can
> cause this? My gdalinfo.exe is not good? Any idea?
 Have a look at the FAQ entries about GDAL_DATA :
 https://trac.osgeo.org/gdal/wiki/FAQInstallationAndBuilding#WhatisGDAL_DATAenvironmentvariable

> C:\Cubes\LC80080472013174LGN00.tar>gdalinfo LC80080472013174LGN00_B1.TIF
> Driver: GTiff/GeoTIFF
> Files: LC80080472013174LGN00_B1.TIF
> .\LC80080472013174LGN00_MTL.txt
> Size is 7571, 7411
> Coordinate System is:
> PROJCS["unnamed",
>  PROJECTION["Transverse_Mercator"],
>  PARAMETER["latitude_of_origin",0],
>  PARAMETER["central_meridian",-69],
>  PARAMETER["scale_factor",0.9996],
>  PARAMETER["false_easting",50],
>  PARAMETER["false_northing",0],
>  UNIT["metre",1,
>  AUTHORITY["EPSG","9001"]]]
> Origin = (145185.000,2189715.000)
> Pixel Size = (30.000,-30.000)
> Metadata:
>AREA_OR_POINT=Point
>METADATATYPE=ODL
> Image Structure Metadata:
>INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (  145185.000, 2189715.000)
> Lower Left  (  145185.000, 1967385.000)
> Upper Right (  372315.000, 2189715.000)
> Lower Right (  372315.000, 1967385.000)
> Center  (  258750.000, 2078550.000)
> Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray
>
> C:\Cubes\LC80080472013174LGN00.tar>gdalinfo --version
> GDAL 2.0.1, released 2015/09/15
>
> C:\Cubes\LC80080472013174LGN00.tar>
>
>> I have the 2.0.1 version:
>> C:\warmerda\bld\bin>gdalinfo --version
>> GDAL 2.0.1, released 2015/09/15
>>
>> I will check the history of the file I have. Maybe it is modified. I
>> will
>> also download the same file from earthexplorer and try my routines on
>> it. Thanks for the suggestions. I appreciate.
>> - Nevzat
>>
>> -Original Message-
>> From: Dmitry Baryshnikov [mailto:bishop@gmail.com]
>> Sent: Sunday, January 17, 2016 12:57 PM
>> To: Nevzat Guler ; gdal-dev@lists.osgeo.org
>> Subject: Re: [gdal-dev] map info and projections from GeoTIFF file to
>> create ENVI header
>>
>> I downloaded this scene from earthexplorer, and everything is fine:
>>
>> $ gdalinfo LC80080472013174LGN00_B1.TIF
>> Driver: GTiff/GeoTIFF
>> Files: LC80080472013174LGN00_B1.TIF
>>
>>  ./LC80080472013174LGN00_MTL.txt
>>
>> Size is 7571, 7411
>> Coordinate System is:
>> PROJCS["WGS 84 / UTM zone 19N",
>>
>>   GEOGCS["WGS 84",
>>
>>   DATUM["WGS_1984",
>>
>>   SPHEROID["WGS 84",6378137,298.257223563,
>>
>>   AUTHORITY["EPSG","7030"]],
>>
>>   AUTHORITY["EPSG","6326"]],
>>
>>

Re: [gdal-dev] map info and projections from GeoTIFF file to create ENVI header

2016-01-18 Thread Dmitry Baryshnikov 

Oh, I'm confused MapInfo and ENVI map info.

In your case you need to form this string by yourself using values from 
OGRSpatialReference methods.
Also take a look to 
https://github.com/OSGeo/gdal/blob/trunk/gdal/frmts/raw/envidataset.cpp


Best regards,
Dmitry

19.01.2016 00:33, ngu...@jlab.org пишет:

Dmitry,

After doing this, I following your suggested steps, I get the string:

Earth Projection 8, 104, "m", -69, 0, 0.9996, 50, 0

Bit, the map info string I need to use looks like this:

map info = {UTM, 1, 1, 145185, 2189715, 30, 30, 19, North}

I tried other ExportTo functions but none gives me the specific map info
string I want. This is ENVI map info string.

I guess I will just have to stitch it together once I understand what the
second and third elements are used for..

I really appreciated your help.
Best regards, - Nevzat



Hi Nevzat,

you need to get spatial reference from GDALDataset (via GetProjectionRef
http://gdal.org/classGDALDataset.html#aa42537e1062ce254d124b29ff3ebe857)
and convert it to OGRSpatialReference.

Than use exportToMICoordSys to get MI CRS string.

http://gdal.org/classOGRSpatialReference.html#a1f2908cd5ca33609844ef0c0ff2186ea

Best regards,
  Dmitry

18.01.2016 20:27, ngu...@jlab.org пишет:

It is solved. Thank you very much Even and Dmitry for your helps. I have
missed the GDAL_DATA variable and it caused me all that trouble.

Now I need to get the map info string:
map info = {UTM, 1, 1, 145185, 2189715, 30, 30, 19, North,WGS-84}

Using hDS->GetGeoTransform(adfGeoTransform) does not give me all the
information in the map info string. It only gives me elements [4-7], how to
get the others?

Best regards, - Nevzat


Yes, Thank you Even, I was missing the GDAL_DATA environmental variable.
Now I
get the same gdalinfo output. Do I need to set this variable while using
the
C++ classes you suggested before: GetGeoTransform and GetProjectionRef()..?
Because, there also I had the same problem.

Best regards, - Nevzat



Le lundi 18 janvier 2016 17:13:14, ngu...@jlab.org a écrit :

Dmitry and All,
I just downloaded the set of tiff files from earthexplorer and run
gdalinfo
on it. Please see below. I am missing some fields in the output
information. It is different from what you got. I don't know what can
cause this? My gdalinfo.exe is not good? Any idea?

Have a look at the FAQ entries about GDAL_DATA :
https://trac.osgeo.org/gdal/wiki/FAQInstallationAndBuilding#WhatisGDAL_DATAenvironmentvariable


C:\Cubes\LC80080472013174LGN00.tar>gdalinfo LC80080472013174LGN00_B1.TIF
Driver: GTiff/GeoTIFF
Files: LC80080472013174LGN00_B1.TIF
 .\LC80080472013174LGN00_MTL.txt
Size is 7571, 7411
Coordinate System is:
PROJCS["unnamed",
  PROJECTION["Transverse_Mercator"],
  PARAMETER["latitude_of_origin",0],
  PARAMETER["central_meridian",-69],
  PARAMETER["scale_factor",0.9996],
  PARAMETER["false_easting",50],
  PARAMETER["false_northing",0],
  UNIT["metre",1,
  AUTHORITY["EPSG","9001"]]]
Origin = (145185.000,2189715.000)
Pixel Size = (30.000,-30.000)
Metadata:
AREA_OR_POINT=Point
METADATATYPE=ODL
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  145185.000, 2189715.000)
Lower Left  (  145185.000, 1967385.000)
Upper Right (  372315.000, 2189715.000)
Lower Right (  372315.000, 1967385.000)
Center  (  258750.000, 2078550.000)
Band 1 Block=7571x1 Type=UInt16, ColorInterp=Gray

C:\Cubes\LC80080472013174LGN00.tar>gdalinfo --version
GDAL 2.0.1, released 2015/09/15

C:\Cubes\LC80080472013174LGN00.tar>


I have the 2.0.1 version:
C:\warmerda\bld\bin>gdalinfo --version
GDAL 2.0.1, released 2015/09/15

I will check the history of the file I have. Maybe it is modified. I
will
also download the same file from earthexplorer and try my routines on
it. Thanks for the suggestions. I appreciate.
- Nevzat

-Original Message-
From: Dmitry Baryshnikov [mailto:bishop@gmail.com]
Sent: Sunday, January 17, 2016 12:57 PM
To: Nevzat Guler ; gdal-dev@lists.osgeo.org
Subject: Re: [gdal-dev] map info and projections from GeoTIFF file to
create ENVI header

I downloaded this scene from earthexplorer, and everything is fine:

$ gdalinfo LC80080472013174LGN00_B1.TIF
Driver: GTiff/GeoTIFF
Files: LC80080472013174LGN00_B1.TIF

  ./LC80080472013174LGN00_MTL.txt

Size is 7571, 7411
Coordinate System is:
PROJCS["WGS 84 / UTM zone 19N",

   GEOGCS["WGS 84",

   DATUM["WGS_1984",

   SPHEROID["WGS 84",6378137,298.257223563,

   AUTHORITY["EPSG","7030"]],

   AUTHORITY["EPSG","6326"]],

   PRIMEM["Greenwich",0,

   AUTHORITY["EPSG","8901"]],

   UNIT["degree",0.0174532925199433,

   AUTHORITY["EPSG","9122"]],

   AUTHORITY["EPSG","4326"]],

   PROJECTION["Transverse_Mercator"],
   PARAMETER["latitude_of_origin",0],

[gdal-dev] OpenJPEG tiled encoding

2016-01-18 Thread Aaron Boxer 
Hello,

It looks like OpenJPEG tiled encoding is broken for non power of 2 images.
This seems to be quite an old bug.

Does this type of encoding get used in GDAL ? Or are all images power of 2
dimensions.

Thanks,
Aaron
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev