Re: [gdal-dev] (no subject) force PAM

[Rahkonen Jukka via gdal-dev]

Hi,

I only managed to force the creation of the .aux.gml file with a workaround.

gdal_translate p4433h.tif test.tif -co profile=baseline

In this case it seems that not even --config GDAL_PAM_ENABLED YES is needed but 
test.tif.aux.xml gets created. Maybe the same PAM file works also for the 
non-baseline versions of that tiff file.

With GeoTIFF, COG, or VRT as output format  --config GDAL_PAM_ENABLED YES does 
not trigger the creation of the PAM file.
I think that storing the same metadata into two distinct places (GeoTIFF tags 
and .aux.xml) is generally speaking a bad idea because sooner or later the two 
versions do not match any more. But I guess that it could be convenient to read 
the plain text aux.xml file from the cloud before touching the main TIFF.

-Jukka Rahkonen-

Lähettäjä: gdal-dev  Puolesta Shaima Alsharif 
via gdal-dev
Lähetetty: perjantai 8. maaliskuuta 2024 11.58
Vastaanottaja: gdal-dev@lists.osgeo.org
Aihe: [gdal-dev] (no subject)

Dear Developers,


I am reaching out to seek advice or potential solutions regarding the creation 
of .aux.xml files for GeoTIFF Or GOG using GDAL. I have attempted using the 
following commands, like so:



os.environ['GDAL_PAM_ENABLED'] = 'YES'

os.environ['GDAL_GEOREF_SOURCES'] = 'WORLDFILE,PAM,INTERNAL'

os.environ['ESRI_XML_PAM'] = 'TRUE'



gdalwarp -overwrite -t_srs EPSG:4326 -r cubic JED240207122504.MAXBERJ  
output.tiff



gdaldem color-relief -alpha output.tiff titan.txt output.colored.tiff



gdal_translate -of COG -co TILING_SCHEME=GoogleMapsCompatible -co COMPRESS=LZW  
output.colored.tiff  output.cog.tiff



gdal_translate -of GTiff -co TILED=YES -co COMPRESS=DEFLATE input.tif 
output.tif --config GDAL_PAM_ENABLED YES



gdalinfo --build
PAM_ENABLED=YES
OGR_ENABLED=YES
GEOS_ENABLED=YES
GEOS_VERSION=3.11.1-CAPI-1.17.1
PROJ_BUILD_VERSION=9.1.1
PROJ_RUNTIME_VERSION=9.1.1
COMPILER=GCC 12.2.0



 gdalinfo --version

GDAL 3.8.4, released 2024/02/08

also tried  gdal-3.6.2



OS is Ubuntu Aarch64





However, this does not result in the creation of an .aux.xml file. Could you 
provide any insights into why the .aux.xml file may not be generated in this 
instance?



Any assistance or suggestions you could offer would be greatly appreciated.



Best regards,

Shaima A.

___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Shaima Alsharif via gdal-dev]

Dear Developers,

I am reaching out to seek advice or potential solutions regarding the
creation of .aux.xml files for GeoTIFF Or GOG using GDAL. I have attempted
using the following commands, like so:


os.environ['GDAL_PAM_ENABLED'] = 'YES'

os.environ['GDAL_GEOREF_SOURCES'] = 'WORLDFILE,PAM,INTERNAL'

os.environ['ESRI_XML_PAM'] = 'TRUE'


gdalwarp -overwrite -t_srs EPSG:4326 -r cubic
JED240207122504.MAXBERJ  output.tiff


gdaldem color-relief -alpha output.tiff titan.txt output.colored.tiff


gdal_translate -of COG -co TILING_SCHEME=GoogleMapsCompatible -co
COMPRESS=LZW  output.colored.tiff  output.cog.tiff


gdal_translate -of GTiff -co TILED=YES -co COMPRESS=DEFLATE input.tif
output.tif --config GDAL_PAM_ENABLED YES


gdalinfo --build
PAM_ENABLED=YES
OGR_ENABLED=YES
GEOS_ENABLED=YES
GEOS_VERSION=3.11.1-CAPI-1.17.1
PROJ_BUILD_VERSION=9.1.1
PROJ_RUNTIME_VERSION=9.1.1
COMPILER=GCC 12.2.0


 gdalinfo --version

GDAL 3.8.4, released 2024/02/08

also tried  gdal-3.6.2


OS is Ubuntu Aarch64



However, this does not result in the creation of an .aux.xml file. Could
you provide any insights into why the .aux.xml file may not be generated in
this instance?


Any assistance or suggestions you could offer would be greatly appreciated.


Best regards,

Shaima A.
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Even Rouault]

Hi,

addressed per https://github.com/OSGeo/gdal/pull/6407

Even

Le 19/09/2022 à 13:07, Moises Calzado via gdal-dev a écrit :


Hello everyone,

We're performing some tests with ogrinfo trying to read GPKG files, 
and we're facing some issues executing the command with remote GPKG 
files in WAL mode. As can be checked in the following command output, 
ogr fails while reading the file content:


ERROR: Error: Command failed: ogrinfo /vsicurl/MY_FILE_URL has
GPKG application_id, but non conformant file extension ERROR 1:
unable to open database file: this file is a WAL-enabled database.
It cannot be opened because it is presumably read-only or in a
read-only directory. ERROR 1: Only read-only mode is supported for
/vsicurl


What is strange here is that the same file works correctly if we 
execute the command using the downloaded file.


Is there any way we can execute the command using a remote file?

Thanks so much in advance,
Kind regards!
--
*Moises Calzado*

Support Engineer

(US) +1 917 463 3232 | (ES) +34 911 165 823 | mcalz...@carto.com



___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev


--
http://www.spatialys.com
My software is free, but my time generally not.
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Rahkonen Jukka]

I have, it does not help if the journal mode is WAL. 

ogrinfo 
/vsicurl/"http://localhost:8060/SuomenHallinnollisetKuntajakopohjaisetAluejaot_2022_10k.gpkg;
 -ro
ERROR 1: unable to open database file: this file is a WAL-enabled database. It 
cannot be opened because it is presumably read-only or in a read-only directory.

-Jukka Rahkonen-

-Alkuperäinen viesti-
Lähettäjä: gdal-dev  Puolesta Andrew C 
Aitchison
Lähetetty: maanantai 19. syyskuuta 2022 17.18
Vastaanottaja: Moises Calzado 
Kopio: gdal-dev@lists.osgeo.org
Aihe: Re: [gdal-dev] (no subject)

On Mon, 19 Sep 2022, Moises Calzado via gdal-dev wrote:

> Hello everyone,
>
> We're performing some tests with ogrinfo trying to read GPKG files, 
> and we're facing some issues executing the command with remote GPKG 
> files in WAL mode. As can be checked in the following command output, 
> ogr fails while reading the file content:
>
> ERROR: Error: Command failed: ogrinfo /vsicurl/MY_FILE_URL has GPKG
>> application_id, but non conformant file extension ERROR 1: unable to 
>> open database file: this file is a WAL-enabled database. It cannot be 
>> opened because it is presumably read-only or in a read-only directory. ERROR 
>> 1:
>> Only read-only mode is supported for /vsicurl
>
>
> What is strange here is that the same file works correctly if we 
> execute the command using the downloaded file.
>
> Is there any way we can execute the command using a remote file?

'man ogrinfo' includes the option:
-roOpen the data source in read-only mode.

So, have you tried
ogrinfo -ro /vsicurl/MY_FILE_URL
?


-- 
Andrew C. Aitchison  Kendal, UK
and...@aitchison.me.uk 
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.osgeo.org%2Fmailman%2Flistinfo%2Fgdal-devdata=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C19520984e6db45c9348708da9a4b1be2%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C0%7C637991944628587871%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7Csdata=vUEpaIT8yFEr%2Bc77371qb2NVR1K1YQNGbBcTBHU70P8%3Dreserved=0
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Andrew C Aitchison]

On Mon, 19 Sep 2022, Moises Calzado via gdal-dev wrote:


Hello everyone,

We're performing some tests with ogrinfo trying to read GPKG files, and
we're facing some issues executing the command with remote GPKG files in
WAL mode. As can be checked in the following command output, ogr fails
while reading the file content:

ERROR: Error: Command failed: ogrinfo /vsicurl/MY_FILE_URL has GPKG

application_id, but non conformant file extension ERROR 1: unable to open
database file: this file is a WAL-enabled database. It cannot be opened
because it is presumably read-only or in a read-only directory. ERROR 1:
Only read-only mode is supported for /vsicurl



What is strange here is that the same file works correctly if we execute
the command using the downloaded file.

Is there any way we can execute the command using a remote file?


'man ogrinfo' includes the option:
   -roOpen the data source in read-only mode.

So, have you tried
ogrinfo -ro /vsicurl/MY_FILE_URL
?


--
Andrew C. Aitchison  Kendal, UK
   and...@aitchison.me.uk
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Rahkonen Jukka]

Could you write a server side script that runs the ogrinfo command with r/w 
rights to the upload directory?

-Jukka Rahkonen-

Lähettäjä: Moises Calzado 
Lähetetty: maanantai 19. syyskuuta 2022 16.00
Vastaanottaja: Rahkonen Jukka 
Aihe: Re: [gdal-dev] (no subject)

Oh no! We use ogrinfo to check the file type, as these files are uploaded by 
other users into a GCS bucket.
We don't apply any change to these files, as we just check them with ogrinfo to 
perform other operations afterwards, so that's why we can't change the journal 
mode before executing the ogrinfo command.
If there's no way that we can complete this without modifying it, we'll think 
about an alternative solution!

Thanks so much for your help!

El lun, 19 sept 2022 a las 14:07, Rahkonen Jukka 
(mailto:jukka.rahko...@maanmittauslaitos.fi>>)
 escribió:
Hi,

What do you mean by processing? Http users cannot edit anything in any case. Do 
you edit the GPKG file locally at the same time? "Immutable" could be a 
solution if GDAL can support it somehow, but if the database is edited at the 
same it does not feel right. Taken from the SQLite documentation:
"immutable: The immutable parameter is a boolean query parameter that indicates 
that the database file is stored on read-only media. When immutable is set, 
SQLite assumes that the database file cannot be changed, even by a process with 
higher privilege, and so the database is opened read-only and all locking and 
change detection is disabled. Caution: Setting the immutable property on a 
database file that does in fact change can result in incorrect query results 
and/or SQLITE_CORRUPT errors. See also: SQLITE_IOCAP_IMMUTABLE."

Do you have some good reason not to change the pragma?

-Jukka Rahkonen-

Lähettäjä: Moises Calzado mailto:mcalz...@carto.com>>
Lähetetty: maanantai 19. syyskuuta 2022 14.46
Vastaanottaja: Rahkonen Jukka 
mailto:jukka.rahko...@maanmittauslaitos.fi>>
Aihe: Re: [gdal-dev] (no subject)

Hi Jukka,

Thanks so much for taking a look at this! We already tried to change the 
journal mode, and it works as expected, but we would like to know if there is 
any way of processing this files without changing anything on the file.

El lun, 19 sept 2022 a las 13:25, Rahkonen Jukka 
(mailto:jukka.rahko...@maanmittauslaitos.fi>>)
 escribió:
Hi,

Read 
https://www.sqlite.org/wal.html<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.sqlite.org%2Fwal.html=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C92f8891853344ee1f57008da9a3ee90a%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C0%7C637991892247543363%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=pa6FZF5AI0XqoyrfkqiqhWhsbfpzBMopKuk%2FSWB9kE4%3D=0>.
 I suggest to change the journal mode into the default mode before inserting 
the GPKG file into server. How to: Execute "PRAGMA journal_mode=DELETE". I 
tested that it will now work:

INFO: Open of 
`/vsicurl/http://localhost:8060/SuomenHallinnollisetKuntajakopohjaisetAluejaot_2022_10k.gpkg'
  using driver `GPKG' successful.
1: Aluehallintovirasto
2: Valtakunta
3: Maakunta
4: Kunta

If the GPKG file in WAL mode is already placed on the server it might be 
possible to open it somehow as "immutable". I do not know how to do that with 
GDAL.

-Jukka Rahkonen-


Lähettäjä: gdal-dev 
mailto:gdal-dev-boun...@lists.osgeo.org>> 
Puolesta Moises Calzado via gdal-dev
Lähetetty: maanantai 19. syyskuuta 2022 14.08
Vastaanottaja: gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org>
Aihe: [gdal-dev] (no subject)


Hello everyone,

We're performing some tests with ogrinfo trying to read GPKG files, and we're 
facing some issues executing the command with remote GPKG files in WAL mode. As 
can be checked in the following command output, ogr fails while reading the 
file content:

ERROR: Error: Command failed: ogrinfo /vsicurl/MY_FILE_URL has GPKG 
application_id, but non conformant file extension ERROR 1: unable to open 
database file: this file is a WAL-enabled database. It cannot be opened because 
it is presumably read-only or in a read-only directory. ERROR 1: Only read-only 
mode is supported for /vsicurl

What is strange here is that the same file works correctly if we execute the 
command using the downloaded file.

Is there any way we can execute the command using a remote file?

Thanks so much in advance,
Kind regards!
--
Moises Calzado

Support Engineer

(US) +1 917 463 3232 | (ES) +34 911 165 823 | 
mcalz...@carto.com<mailto:mcalz...@carto.com>
[https://storage.googleapis.com/carto-it-files/signature/SDSC-2022-NY_Signature-mini.jpg]<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fspatial-data-science-conference.com%2F2022%2Fnewyork%2F=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C92f8891853344ee1f57008da9a3ee90a%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C0%7C637991892247543363%7CUnknown%7CTWF

Re: [gdal-dev] (no subject)

[Rahkonen Jukka]

Hi,

What do you mean by processing? Http users cannot edit anything in any case. Do 
you edit the GPKG file locally at the same time? "Immutable" could be a 
solution if GDAL can support it somehow, but if the database is edited at the 
same it does not feel right. Taken from the SQLite documentation:
"immutable: The immutable parameter is a boolean query parameter that indicates 
that the database file is stored on read-only media. When immutable is set, 
SQLite assumes that the database file cannot be changed, even by a process with 
higher privilege, and so the database is opened read-only and all locking and 
change detection is disabled. Caution: Setting the immutable property on a 
database file that does in fact change can result in incorrect query results 
and/or SQLITE_CORRUPT errors. See also: SQLITE_IOCAP_IMMUTABLE."

Do you have some good reason not to change the pragma?

-Jukka Rahkonen-

Lähettäjä: Moises Calzado 
Lähetetty: maanantai 19. syyskuuta 2022 14.46
Vastaanottaja: Rahkonen Jukka 
Aihe: Re: [gdal-dev] (no subject)

Hi Jukka,

Thanks so much for taking a look at this! We already tried to change the 
journal mode, and it works as expected, but we would like to know if there is 
any way of processing this files without changing anything on the file.

El lun, 19 sept 2022 a las 13:25, Rahkonen Jukka 
(mailto:jukka.rahko...@maanmittauslaitos.fi>>)
 escribió:
Hi,

Read 
https://www.sqlite.org/wal.html<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.sqlite.org%2Fwal.html=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C778702221d7043994bb408da9a3484f6%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C0%7C637991847599144072%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=dj6xkZhVJSYck6IqCVx85KcTxIuJHWR34O0BZ%2BNU3io%3D=0>.
 I suggest to change the journal mode into the default mode before inserting 
the GPKG file into server. How to: Execute "PRAGMA journal_mode=DELETE". I 
tested that it will now work:

INFO: Open of 
`/vsicurl/http://localhost:8060/SuomenHallinnollisetKuntajakopohjaisetAluejaot_2022_10k.gpkg'
  using driver `GPKG' successful.
1: Aluehallintovirasto
2: Valtakunta
3: Maakunta
4: Kunta

If the GPKG file in WAL mode is already placed on the server it might be 
possible to open it somehow as "immutable". I do not know how to do that with 
GDAL.

-Jukka Rahkonen-


Lähettäjä: gdal-dev 
mailto:gdal-dev-boun...@lists.osgeo.org>> 
Puolesta Moises Calzado via gdal-dev
Lähetetty: maanantai 19. syyskuuta 2022 14.08
Vastaanottaja: gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org>
Aihe: [gdal-dev] (no subject)


Hello everyone,

We're performing some tests with ogrinfo trying to read GPKG files, and we're 
facing some issues executing the command with remote GPKG files in WAL mode. As 
can be checked in the following command output, ogr fails while reading the 
file content:

ERROR: Error: Command failed: ogrinfo /vsicurl/MY_FILE_URL has GPKG 
application_id, but non conformant file extension ERROR 1: unable to open 
database file: this file is a WAL-enabled database. It cannot be opened because 
it is presumably read-only or in a read-only directory. ERROR 1: Only read-only 
mode is supported for /vsicurl

What is strange here is that the same file works correctly if we execute the 
command using the downloaded file.

Is there any way we can execute the command using a remote file?

Thanks so much in advance,
Kind regards!
--
Moises Calzado

Support Engineer

(US) +1 917 463 3232 | (ES) +34 911 165 823 | 
mcalz...@carto.com<mailto:mcalz...@carto.com>
[https://storage.googleapis.com/carto-it-files/signature/SDSC-2022-NY_Signature-mini.jpg]<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fspatial-data-science-conference.com%2F2022%2Fnewyork%2F=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C778702221d7043994bb408da9a3484f6%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C0%7C637991847599144072%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=m0xGittslja0%2BX5e8yNAmv6XZSk9FQJU0yNEZdZwO9I%3D=0>


--
Moises Calzado

Support Engineer

(US) +1 917 463 3232 | (ES) +34 911 165 823 | 
mcalz...@carto.com<mailto:mcalz...@carto.com>
[https://storage.googleapis.com/carto-it-files/signature/SDSC-2022-NY_Signature-mini.jpg]<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fspatial-data-science-conference.com%2F2022%2Fnewyork%2F=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C778702221d7043994bb408da9a3484f6%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C0%7C637991847599144072%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=m0xGittslja0%2BX5e8yNAmv6XZSk9FQJU0yNEZdZwO9I%3D=0>
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Rahkonen Jukka]

Hi,

Read https://www.sqlite.org/wal.html. I suggest to change the journal mode into 
the default mode before inserting the GPKG file into server. How to: Execute 
"PRAGMA journal_mode=DELETE". I tested that it will now work:

INFO: Open of 
`/vsicurl/http://localhost:8060/SuomenHallinnollisetKuntajakopohjaisetAluejaot_2022_10k.gpkg'
  using driver `GPKG' successful.
1: Aluehallintovirasto
2: Valtakunta
3: Maakunta
4: Kunta

If the GPKG file in WAL mode is already placed on the server it might be 
possible to open it somehow as "immutable". I do not know how to do that with 
GDAL.

-Jukka Rahkonen-


Lähettäjä: gdal-dev  Puolesta Moises Calzado 
via gdal-dev
Lähetetty: maanantai 19. syyskuuta 2022 14.08
Vastaanottaja: gdal-dev@lists.osgeo.org
Aihe: [gdal-dev] (no subject)


Hello everyone,

We're performing some tests with ogrinfo trying to read GPKG files, and we're 
facing some issues executing the command with remote GPKG files in WAL mode. As 
can be checked in the following command output, ogr fails while reading the 
file content:

ERROR: Error: Command failed: ogrinfo /vsicurl/MY_FILE_URL has GPKG 
application_id, but non conformant file extension ERROR 1: unable to open 
database file: this file is a WAL-enabled database. It cannot be opened because 
it is presumably read-only or in a read-only directory. ERROR 1: Only read-only 
mode is supported for /vsicurl

What is strange here is that the same file works correctly if we execute the 
command using the downloaded file.

Is there any way we can execute the command using a remote file?

Thanks so much in advance,
Kind regards!
--
Moises Calzado

Support Engineer

(US) +1 917 463 3232 | (ES) +34 911 165 823 | 
mcalz...@carto.com<mailto:mcalz...@carto.com>
[https://storage.googleapis.com/carto-it-files/signature/SDSC-2022-NY_Signature-mini.jpg]<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fspatial-data-science-conference.com%2F2022%2Fnewyork%2F=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C59c27b74927549e7612f08da9a2f47d9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C0%7C637991825108785824%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=hAhLefsE2N2%2ByrpGgP%2F%2FH1Hwqi9WeTUK6Onk10RkbsM%3D=0>
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Moises Calzado via gdal-dev]

Hello everyone,

We're performing some tests with ogrinfo trying to read GPKG files, and
we're facing some issues executing the command with remote GPKG files in
WAL mode. As can be checked in the following command output, ogr fails
while reading the file content:

ERROR: Error: Command failed: ogrinfo /vsicurl/MY_FILE_URL has GPKG
> application_id, but non conformant file extension ERROR 1: unable to open
> database file: this file is a WAL-enabled database. It cannot be opened
> because it is presumably read-only or in a read-only directory. ERROR 1:
> Only read-only mode is supported for /vsicurl


What is strange here is that the same file works correctly if we execute
the command using the downloaded file.

Is there any way we can execute the command using a remote file?

Thanks so much in advance,
Kind regards!
-- 
*Moises Calzado*

Support Engineer

(US) +1 917 463 3232 | (ES) +34 911 165 823 | mcalz...@carto.com

___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Ahmet Temiz]

Hi,
How can I can create FeatureCollection  from :
```OGRMultiPolygon * poOGRMultiPolygon= (OGRMultiPolygon * )
unionCascaded;```


kind regards


-- 
Ahmet Temiz
Jeoloji Müh.
Afet ve Acil Durum Yönetimi Başkanlığı
Deprem  Dairesi Başkanlığı




Ahmet Temiz
Geological Eng.

Disaster and Emergency Management
of Presidency
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Ahmet Temiz]

Hi,

I need example c++ code for raster masking and clipping.
Do you have any?

Regards
-- 
Ahmet Temiz
Jeoloji Müh.
Afet ve Acil Durum Yönetimi Başkanlığı
Bilgi İşlem  Dairesi Başkanlığı-CBS Grubu




Ahmet Temiz
Geological Eng.
Information Systems - GIS Group
Disaster and Emergency Management
of Presidency
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Mateusz Loskot]

On 18 April 2018 at 13:57, Gane R  wrote:
> Hi all,
>
> Is there any way to disable one or all OGR formats when building on Linux.
> i.e.,
> when using ./configure
>
> Is there any way to disable certain raster formats when building on Linux

https://trac.osgeo.org/gdal/wiki/BuildingOnUnixWithMinimizedDrivers

-- 
Mateusz Loskot, http://mateusz.loskot.net
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Gane R]

Hi all,

Is there any way to disable one or all OGR formats when building on Linux.
i.e.,
when using ./configure

Is there any way to disable certain raster formats when building on Linux

Thanks
Ganesan
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Victor Chernetsky]

Hello,

Is there future GDAL release schedule available anywhere? I am actually
wondering about when is GDAL 2.2 is planned to be released? And are there
any more of 2.1.x releases ahead?

Thank you,
*Victor Chernetsky*
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Even Rouault]

On Tuesday 26 July 2016 16:37:14 Hema Yeedunuri wrote:
> aded images which uses wms services. I tried opening with
> gdalopen it is not working . the output format supported in geoserver is
> png / geotiff . but for a single band data the no. of bands i am getting
> are 3 . and null data is coming.

Hema,

You can change the number of bands returned by the WMS driver by setting the 
 element. See http://gdal.org/frmt_wms.html

Even

-- 
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] (no subject)

[Hema Yeedunuri]

hi all ,

I have build gdal-1.11.2 with libcurl support to read wms layers.it is
build sucessfully how can i read wms layers .I have geoserver on my
loalhost.I have loaded images which uses wms services. I tried opening with
gdalopen it is not working . the output format supported in geoserver is
png / geotiff . but for a single band data the no. of bands i am getting
are 3 . and null data is coming.


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

[gdal-dev] (no subject)

[Neil Oster]

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

[gdal-dev] (no subject)

[user gdal]

Hello,

I have a shape file which consists of many polygons (say 56).
Programmatically I wish to select (say, one poly.) from them. I proceeded
something like this.

OGRLayer *clippingshp; unsigned int NumPolygons;

while((tempfeature=inlayer->GetNextFeature( ))!=NULL)
{
  clippingshp->CreateFeature(tempfeature);
clippingshp->SetFeature(tempfeature);


  NumPolygons=clippingshp->GetFeatureCount( ); /* After this NumPolygons
should be 1 */

  /* But NumPolygons = 56, that is the whole file is selected */
  .
  .
  .
}

Can anyone tell me please how to select a feature (polygon in this context)
and copy it into an OGRLayer.


With many thanks,

Yours sincerely,
  B. Ramesh
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Константин Лапин]

Hi,
CPLSetConfigOption( "GDAL_HTTP_PROXY",proxy.toLocal8Bit().data());
CPLSetConfigOption( "GDAL_PROXY_AUTH", "NTLM" );
CPLSetConfigOption( "GDAL_HTTP_PROXYUSERPWD", " : " );
how to make "bypass proxi for local adresses" ?
С уважением,
Константин Лапин
lk...@mail.ru___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Dmitry Baryshnikov]

Hi Konstantin,

Can you explain your case. For example for WMS driver (TMS) 
GDAL_HTTP_PROXY config option works for me.


Best regards,
Dmitry

01.10.2015 14:39, Константин Лапин пишет:


Hi,

CPLSetConfigOption( "GDAL_HTTP_PROXY",proxy.toLocal8Bit().data());
CPLSetConfigOption( "GDAL_PROXY_AUTH", "NTLM" );
CPLSetConfigOption( "GDAL_HTTP_PROXYUSERPWD", " : " );

how to make "bypass proxi for local adresses" ?


С уважением,
Константин Лапин
lk...@mail.ru


___
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] (no subject)

[user gdal]

Dear Dr. Brad, Thanks for your kind and prompt reply. Below are my
answers/clarifications.

On Wed, 19 Aug 2015 04:01:32 PM user gdal wrote:
 Dear friends,
I have to programmatically clip (using Visual C++) with
 GDALWarp. I proceeded like this:

 GDALDataset *outimg;

 outimg-(GDALDataset *) GDALOpen(filename, GA_Update);
 /* I also tried, instead of the above, outimg=driver-Create(filename,
 NewCol, NewRow, Band, DataType, NULL); */

 GDALWarpOptions *WO=GDALCreateWarpOptions( );

 if(WO)
 {
   WO-hSrcDS=inimg; WO-nBandCount=Band;

What is inimg?
Ans. The bigger input image from which I want to draw subset (i.e., clip)
with the help of a vector layer myshp. It' (inimg) is of type GDALDataset
*. In the declaration of variables, it would appear as GDALDataset *inimg
Also, in WO-hCutline (pl. see below), I am passing myshp (which has
strictly polygon topology), of the type OGRLayer * (i.e., OGRLayer *myshp).
Of course the shape file and inimg have the same projection.


What (exactly) are you passing as filename?
Ans. char *filename (such as e:\sub\subimage.img. There are no other
characters in the 'filename')

Can you try posting an actual, compilable example?
Ans. There is a GUI from which I extract parameters which can give us a lot
of inconvenience. The program compiles without 'abusing' me. What I am
getting is a run time (or debug time) error.

--
GDALDataset *inimg, *outimg;

outimg-(GDALDataset *) GDALOpen(filename, GA_Update);
/* I also tried, instead of the above, outimg=driver-Create(filename,
NewCol, NewRow, Band, DataType, NULL); NewCol and NewRow, have values of
the range 30 each */

GDALWarpOptions *WO=GDALCreateWarpOptions( );

if(WO)
{
  WO-hSrcDS=inimg; WO-nBandCount=Band;

  // In a sample image I wish to use, Band =3, checked

  WO-hDstDS=(GDALDatasetH)(outimg); WO-hCutline=myshp;

  WO-panSrcBands=(int *) CPLMalloc(sizeof(int)*Band);
WO-panSrcBands[0]=Band;
  WO-panDstBands=(int *) CPLMalloc(sizeof(int)*Band);
WO-panDstBands[0]=Band;

  WO-pfnProgress=GDALTermProgress;

  GDALWarpOperation *operation=new GDALWarpOperation;

  if(operation)
  {
operation-Initialize(WO);
operation-ChunkAndWarpImage(0, 0, NewCol, NewRow);

GDALDestroyWarpOptions(WO); WO=NULL; delete operation; operation=NULL;

   }

   GDALClose((GDALDatasetH) outimg); outimg=NULL;

  }

--
Does the tutorial (just by itself, exactly as written) work on your machine?
Ans. To be honest I haven't checked. I, however, have gone through the
tutorial and after convincing myself that I have understood what to do, I
wrote the program.

__
__

With many thanks,
Yours sincerely,
 T.S. Prasad
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[user gdal]

Friends I have a raster image (img) file with a coordinate system. How to
get the information through program, the X and Y spatial resolutions?

Thanks,

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

Re: [gdal-dev] (no subject)

[Ramiro Marco Figuera]

Dear Ramesh,

I guess that using gdalinfo nameoffile should work.

Cheers,

Ramiro

On 07/23/2015 12:30 PM, user gdal wrote:
Friends I have a raster image (img) file with a coordinate system. How 
to get the information through program, the X and Y spatial resolutions?


Thanks,

B. Ramesh



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


--
Ramiro Marco Figuera
Department of Physics and Earth Sciences
Jacobs University Bremen
Campus Ring 1, 28759 Bremen, Germany
Office: Research III, Room 99b
Tel: +49 (0)421 200 3226

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

[gdal-dev] (no subject)

[Even Rouault]

Hi PSC,

Mikhail Gusev, one of our 2 GSoC student this year is about to submit a pull
request for the initial part of his - good - work on Geography Network Model.
See https://github.com/MikhanGusev/gdal/compare/OSGeo:trunk...trunk

This add a new set of API to the code base. Looking at our management rules
http://trac.osgeo.org/gdal/wiki/rfc1_pmc (When is Vote Required?: [...] Adding
substantial amounts of new code) and
http://trac.osgeo.org/gdal/wiki/rfc3_commiters (larger changes will require a
RFC approved by the PSC), it seems it would imply writing a RFC. We don't
generally do RFC for new drivers (that can add substantial amount of code), but
technically, this is not really a new driver either. Any opinion on how to
proceed on that matter ?

Best regards,

Even

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

[gdal-dev] (no subject)

[Jonathan Greenberg]

GDALers:

I've not had this error before when compiling GDAL for Linux boxes, but I'm
getting some sqlite error.  This is from:
./configure --prefix=/pathto/wherever --with=python
make

***
...
chmod a+x gdal-config
/bin/sh /projects/oa/shared/jgrn/code/src/gdal-1.10.1/libtool --mode=link
g++  gdalinfo.lo commonutils.lo
 /projects/oa/shared/jgrn/code/src/gdal-1.10.1/libgdal.la -o gdalinfo
libtool: link: g++ .libs/gdalinfo.o .libs/commonutils.o -o .libs/gdalinfo
 /projects/oa/shared/jgrn/code/src/gdal-1.10.1/.libs/libgdal.so
-L/usr/apps/oa/lib -L/usr/lib64/hdf -L/usr/apps/oa/lib/lib
-L/usr/local/openmpi-1.4-gcc/lib /usr/apps/oa/lib/libgeos_c.so
/usr/apps/oa/lib/libgeos.so /usr/apps/oa/lib/libsqlite3.so -lexpat -ljasper
-lnetcdf -lhdf5 /usr/apps/oa/lib/libmfhdf.so /usr/apps/oa/lib/libdf.so
/usr/apps/oa//lib/libjpeg.so /usr/apps/oa/lib/libpng15.so -lpthread -lrt
-lpcre -lcurl /usr/apps/oa/lib/libxml2.so -lm -ldl -lz  -Wl,-rpath
-Wl,/usr/apps/oa/lib -Wl,-rpath -Wl,/usr/apps/oa//lib
/projects/oa/shared/jgrn/code/src/gdal-1.10.1/.libs/libgdal.so: undefined
reference to `sqlite3_column_table_name'
collect2: ld returned 1 exit status
make[1]: *** [gdalinfo] Error 1
make[1]: Leaving directory
`/projects/oa/shared/jgrn/code/src/gdal-1.10.1/apps'
make: *** [apps-target] Error 2


***

I do have SQLite properly installed via the sqlite-autoconf-3080200
release, configured via:
./configure --prefix=/pathto/wherever
make
make install

Any ideas?

--j

-- 
Jonathan A. Greenberg, PhD
Assistant Professor
Global Environmental Analysis and Remote Sensing (GEARS) Laboratory
Department of Geography and Geographic Information Science
University of Illinois at Urbana-Champaign
259 Computing Applications Building, MC-150
605 East Springfield Avenue
Champaign, IL  61820-6371
Phone: 217-300-1924
http://www.geog.illinois.edu/~jgrn/
AIM: jgrn307, MSN: jgrn...@hotmail.com, Gchat: jgrn307, Skype: jgrn3007
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject) [sqlite configuration]

[Even Rouault]

Le jeudi 23 janvier 2014 17:09:41, Jonathan Greenberg a écrit :
 GDALers:
 
 I've not had this error before when compiling GDAL for Linux boxes, but I'm
 getting some sqlite error.  This is from:
 ./configure --prefix=/pathto/wherever --with=python
 make
 
 ***
 ...
 chmod a+x gdal-config
 /bin/sh /projects/oa/shared/jgrn/code/src/gdal-1.10.1/libtool --mode=link
 g++  gdalinfo.lo commonutils.lo
  /projects/oa/shared/jgrn/code/src/gdal-1.10.1/libgdal.la -o gdalinfo
 libtool: link: g++ .libs/gdalinfo.o .libs/commonutils.o -o .libs/gdalinfo
  /projects/oa/shared/jgrn/code/src/gdal-1.10.1/.libs/libgdal.so
 -L/usr/apps/oa/lib -L/usr/lib64/hdf -L/usr/apps/oa/lib/lib
 -L/usr/local/openmpi-1.4-gcc/lib /usr/apps/oa/lib/libgeos_c.so
 /usr/apps/oa/lib/libgeos.so /usr/apps/oa/lib/libsqlite3.so -lexpat -ljasper
 -lnetcdf -lhdf5 /usr/apps/oa/lib/libmfhdf.so /usr/apps/oa/lib/libdf.so
 /usr/apps/oa//lib/libjpeg.so /usr/apps/oa/lib/libpng15.so -lpthread -lrt
 -lpcre -lcurl /usr/apps/oa/lib/libxml2.so -lm -ldl -lz  -Wl,-rpath
 -Wl,/usr/apps/oa/lib -Wl,-rpath -Wl,/usr/apps/oa//lib
 /projects/oa/shared/jgrn/code/src/gdal-1.10.1/.libs/libgdal.so: undefined
 reference to `sqlite3_column_table_name'
 collect2: ld returned 1 exit status
 make[1]: *** [gdalinfo] Error 1
 make[1]: Leaving directory
 `/projects/oa/shared/jgrn/code/src/gdal-1.10.1/apps'
 make: *** [apps-target] Error 2
 
 
 ***
 
 I do have SQLite properly installed via the sqlite-autoconf-3080200
 release, configured via:
 ./configure --prefix=/pathto/wherever
 make
 make install
 
 Any ideas?

Jonathan,

The GDAL source code tries to use sqlite3_column_table_name only if it is 
available i the sqlite library. Indeed, not all builds of sqlite3 have 
sqlite3_column_table_name defined. So you should not get this error 
theroretically.

My hypothesis is that the test in ./configure that detects if 
sqlite3_column_table_name is available or not must link against another 
library than the one that GDAL is finally linking against. This is a bit weird 
however that this happens, but I'd encourage you checking that there are not 
several libsqlite3.so hanging around.

Even

-- 
Geospatial professional services
http://even.rouault.free.fr/services.html
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Phelippe Neveu]

To whom it may concern,

We have updates to the PCIDSK library that we would like to upstream to GDAL.

We have fixed numerous bugs, we have done many optimizations, and we also 
implemented a new binary tile directory.

How should we proceed with the upstream?

Thanks,

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

Re: [gdal-dev] (no subject)

[Frank Warmerdam]

Phelippe,

There are a few ways we could move forward on this.

I believe the apparent master version of the PCIDSK SDK from GDAL's point
of view is:

  http://svn.osgeo.org/gdal/sandbox/warmerdam/pcidsk

And I think the process I used before was to periodally migrate this into
the GDAL tree at:

  http://svn.osgeo.org/gdal/trunk/gdal/frmts/pcidsk/sdk/

Approaches:

1) You just drop a code dump on us, and existing developers integrate it.

2) We arrange commit access for you and you prepare an updated SDK in the
sandbox area in SVN.  We could then migrate that into the GDAL PCIDSK
driver itself when we are comfortable with it.

3) We arrange full commit access for you and you do your integration
directly in GDAL itself.

I prefer option (2) which avoids the likely breakage in GDAL that might
come from a wholesale replacement of the existing inside-GDAL code.  But
those who wish can still built against the standalone SDK for testing
purposes and so forth.  I'd imagine I'd be the GDAL developer to help move
things over eventually.

If you are ok with (2) then please create an OSGeo at:

  http://www.osgeo.org/osgeo_userid/

And I'll try to arrange sandbox-only commit access for you.

Best regards,
Frank



On Fri, Nov 8, 2013 at 12:19 PM, Phelippe Neveu ne...@pcigeomatics.comwrote:

  To whom it may concern,

  We have updates to the PCIDSK library that we would like to upstream to
 GDAL.

  We have fixed numerous bugs, we have done many optimizations, and we
 also implemented a new binary tile directory.

  How should we proceed with the upstream?

  Thanks,

  Phelippe

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




-- 
---+--
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 Software Developer
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Federico Jurio]

Hi all, I have a csv file named test.csv with the asociated test.vrt

*test.vrt*
OGRVRTDataSource
OGRVRTLayer name=test

SrcDataSourceC:\\Users\\FedericoJurio\\Desktop\\test.csv/SrcDataSource
GeometryTypewkbPoint/GeometryType
LayerSRSWGS84/LayerSRS
GeometryField encoding=PointFromColumns x=Longitude
y=Latitude/
/OGRVRTLayer
/OGRVRTDataSource

*test.csv*
Latitude,Longitude,Name
48.1,0.25,First point
49.2,1.1,Second point
47.5,0.75,Third point

In my code i can add the vrt file to the datasource successfully but when i
try to check the TestCapability(OLCRandomRead) from the correct layer
returns FALSE and i don't no why. Gdal can read the file from the console
with ogrinfo and everithing seems to be fine

*ogrinfo test.vrt*
INFO: Open of `test.vrt'
  using driver `VRT' successful.
1: test (Point)

Does anyone have any idea of why TestCapability returns false?

Thanks in advance,
Regards

-- 
Federico Jurio
SUR Emprendimientos Tecnológicos

Perú 345  Piso 5to Oficina B (C1067AAG)
Ciudad de Buenos Aires, Argentina
Tel. +54 (11) 4342-2976/84
federicoju...@suremptec.com.ar
www.suremptec.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Even Rouault]

Le vendredi 08 novembre 2013 22:41:07, Federico Jurio a écrit :
 Hi all, I have a csv file named test.csv with the asociated test.vrt
 
 *test.vrt*
 OGRVRTDataSource
 OGRVRTLayer name=test
 
 SrcDataSourceC:\\Users\\FedericoJurio\\Desktop\\test.csv/SrcDataSource
 GeometryTypewkbPoint/GeometryType
 LayerSRSWGS84/LayerSRS
 GeometryField encoding=PointFromColumns x=Longitude
 y=Latitude/
 /OGRVRTLayer
 /OGRVRTDataSource
 
 *test.csv*
 Latitude,Longitude,Name
 48.1,0.25,First point
 49.2,1.1,Second point
 47.5,0.75,Third point
 
 In my code i can add the vrt file to the datasource successfully but when i
 try to check the TestCapability(OLCRandomRead) from the correct layer
 returns FALSE and i don't no why. Gdal can read the file from the console
 with ogrinfo and everithing seems to be fine
 
 *ogrinfo test.vrt*
 INFO: Open of `test.vrt'
   using driver `VRT' successful.
 1: test (Point)
 
 Does anyone have any idea of why TestCapability returns false?

Quoting http://gdal.org/ogr/ogr__api_8h.html#a480adc8b839b04597f49583371d366fd 
:

OLCRandomRead / RandomRead: TRUE if the GetFeature() method is 
implemented in an optimized way for this layer, as opposed to the default 
implementation using ResetReading() and GetNextFeature() to find the requested 
feature id.

The organization of a CSV file doesn't allow *efficient* random access, hence 
TestCapability(OLCRandomRead) return FALSE. This doesn't prevent random 
reading from being possible. It will be just slow (on a file with several 
thousands of records... not a 3 line file ;-)). And ogrinfo just needs 
efficient 
sequential reading, which CSV allows.


 
 Thanks in advance,
 Regards

-- 
Geospatial professional services
http://even.rouault.free.fr/services.html
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Phelippe Neveu]

Hello Frank,

We would probably go with the approach #2. We will also have to merge any 
changes made in GDAL if any changes were done since our last update which was 2 
years ago I think.

I'm going on vacation for the next week, so I will get back to you on the 18th 
of November to arrange the access to the sandbox.

Regards,

Phelippe

From: fwarmer...@gmail.com fwarmer...@gmail.com on behalf of Frank Warmerdam 
warmer...@pobox.com
Sent: Friday, November 08, 2013 16:34
To: Phelippe Neveu
Cc: gdal-dev@lists.osgeo.org; David Stanley; Louis Burry
Subject: Re: [gdal-dev] (no subject)

Phelippe,

There are a few ways we could move forward on this.

I believe the apparent master version of the PCIDSK SDK from GDAL's point of 
view is:

  http://svn.osgeo.org/gdal/sandbox/warmerdam/pcidsk

And I think the process I used before was to periodally migrate this into the 
GDAL tree at:

  http://svn.osgeo.org/gdal/trunk/gdal/frmts/pcidsk/sdk/

Approaches:

1) You just drop a code dump on us, and existing developers integrate it.

2) We arrange commit access for you and you prepare an updated SDK in the 
sandbox area in SVN.  We could then migrate that into the GDAL PCIDSK driver 
itself when we are comfortable with it.

3) We arrange full commit access for you and you do your integration directly 
in GDAL itself.

I prefer option (2) which avoids the likely breakage in GDAL that might come 
from a wholesale replacement of the existing inside-GDAL code.  But those who 
wish can still built against the standalone SDK for testing purposes and so 
forth.  I'd imagine I'd be the GDAL developer to help move things over 
eventually.

If you are ok with (2) then please create an OSGeo at:

  http://www.osgeo.org/osgeo_userid/

And I'll try to arrange sandbox-only commit access for you.

Best regards,
Frank



On Fri, Nov 8, 2013 at 12:19 PM, Phelippe Neveu 
ne...@pcigeomatics.commailto:ne...@pcigeomatics.com wrote:
To whom it may concern,

We have updates to the PCIDSK library that we would like to upstream to GDAL.

We have fixed numerous bugs, we have done many optimizations, and we also 
implemented a new binary tile directory.

How should we proceed with the upstream?

Thanks,

Phelippe

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



--
---+--
I set the clouds in motion - turn up   | Frank Warmerdam, 
warmer...@pobox.commailto:warmer...@pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush| Geospatial Software Developer
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[user gdal]

Friends,
I sent this mail on Friday. Request for a remedy. If someone can send a
built-in project (needn't be with a high end code) using MFC in VC++ 2008
in a Win Zip or RAR, it will be even better.

The error  (pl. see No. 8) I got is
f:\sw\GDALexperiment\GDALexperiment\header\cpl_port.h: fatal error C1189:
#error :  Unexpected value for SIZEOF_UNSIGNED_LONG GDALexperiment.cpp

I wish to use MFC which is part of Visual C++ 2008 and GDAL.


1. I downloaded GDAL 'latest stable release'  as a zip file from
http://trac.osgeo.org/gdal/wiki/DownloadSource


2. Extracted the zip file. Compiled the project to get gdal110.dll + many
.lib files + many .h files


3. Created a folder in my project ('ref'). Copied gdal110.dll into the ref
folder


4. Likewise created two folders 'header' and 'lib'. Copied lib files from
the download to lib and h files
to header folders


5. In Tools|Options in Visual Studio in my project: accessed Projects and
Solutions: VC++  directories and mentioned header (Include Files), ref
(Reference Files) and lib (library files)


6. In CMyProjectdoc.h, included gdal_priv.h and cpl_conv.h just below #
pragma once


7. Declared a variable in the same file GDALDataset *indataset;


8. Compiled (without writing anything further) and got error(s): Unexpected
values...


Can anyone teach me where I went wrong and how to proceed correctly.



Thanking you,



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

Re: [gdal-dev] (no subject)

[Jukka Rahkonen]

Ahmet Temiz ahmettemiz88 at gmail.com writes:

 
 
 hello
 
 can you tell me what is wrong with following expression:
 (most likely regular expression)
 
 ogr2ogr  -where DETAY!='YRD_MUNHANI' or DETAY!='PAFTA_KENARI' or
DETAY!='KIYI' dene.shp h22a1top.shp


Hi,

Nothing wrong really, but because you use OR you will select everything
and that is probably not what you expect. I suppose your aim is to select
all features except these three: YRD_MUNHANI, PAFTA_KENARI, and KIYI.

DETAY=YRD_MUNHANI will be rejected by the first clause, but because it is
not PAFTA_KENARI it will be reselected by the next clause. Try with ANDs.
This is more understandable to read and seems to work also
ogrinfo -where DETAY not in ('YRD_MUNHANI','PAFTA_KENARI') wheretest.shp -al

If things go more complicated, use -sql instead of -where. If you select
also to use SQLite dialect by adding -dialect sqlite you can use standard
ANSI SQL.


-Jukka Rahkonen-


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

[gdal-dev] (no subject)

[Ahmet Temiz]

hello

can you tell me what is wrong with following expression:
(most likely regular expression)


ogr2ogr  -where DETAY!='YRD_MUNHANI' or DETAY!='PAFTA_KENARI' or
DETAY!='KIYI' dene.shp h22a1top.shp


regards

-- 
Ahmet Temiz
Jeoloji Müh.
Afet ve Acil Durum Yönetimi Başkanlığı
Planlama ve Zarar Azaltma Dairesi Başkanlığı




Ahmet Temiz
Geological Eng.
Information Systems - GIS Group
Disaster and Emergency Management
of Presidency
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Livneh Yehiyam]

Hello list
I'm using gdal on android.
I'm having problems reading data from tiff files larger than 2GB.
Smaller files works ok.
Is there a limit in android, or the android GDAL build? Is there some build 
flag to make it work?
Thanks
Yehiyam Livneh

Sent from my mobile



This message (including any attachments) issued by RAFAEL- ADVANCED DEFENSE 
SYSTEMS LTD. (hereinafter RAFAEL) contains confidential information intended 
for a specific individual and purpose, may constitute information that is 
privileged or confidential or otherwise protected from disclosure. If you are 
not the intended recipient, you should contact us immediately and thereafter 
delete this message from your system. You are hereby notified that any 
disclosure, copying, dissemination, distribution or forwarding of this message, 
or the taking of any action based on it, is strictly prohibited. If you have 
received this e-mail in error, please notify us immediately by e-mail 
mailto:law...@rafael.co.il and completely delete or destroy any and all 
electronic or other copies of the original message and any attachments thereof.
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[lzzgeo]

http://etrange-nuitenchantee.com/market/aGet%20Rich!%20Make%20Money%20Fast%20and%20Get%20Rich%20Quick

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

[gdal-dev] (no subject)

[Nikhil Sai Parupalli]

Hi All,

With the some useful tips and code from  gdal-dev@lists.osgeo.org  I could 
successfully complete a huge task to validate S58 checks .

Next up is developing a software where it creates a S57 data set .
Please give me any information pertaining to S57 data set creation

Thanks and Regards
Nikhil Sai Parupalli



Note: Do not print this email until and unless it is really required. Save 
paper , stay Green
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[SIVA RAMA KRISHNA]

s.r.kriis...@gmail.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Nikhil Sai Parupalli]

Hi All,

How do we get to know which S57 dataset has update sin it..


Thanks and Regards
Nikhil Sai Parupalli



Note: Do not print this email until and unless it is really required. Save 
paper , stay Green
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Dennis Burgess]

 

 



---
Dennis Burgess, Mikrotik Certified Trainer Author of Learn RouterOS-
Second Edition http://www.wlan1.com/product_p/mikrotik%20book-2.htm 
Link Technologies, Inc -- Mikrotik  WISP Support Services
Office: 314-735-0270 tel:314-735-0270  Website:
http://www.linktechs.net http://www.linktechs.net/  - Skype: linktechs
skype:linktechs?call 
-- Create Wireless Coverages with www.mywificoverage.com
http://www.mywificoverage.com/  -- 


--

 

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

[gdal-dev] (no subject)

[yehya imam]

.It seems to be the biggest on-line shop! It’s fantastic!  
http://phbalanse.no/friends.group.php?yfs=07r3
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Wright, Patrick]

smime.p7m
Description: S/MIME encrypted message
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Peisheng Zhao]

I hope you’ll enjoy after visiting this site... 
http://comitemarnebillard.com/friends_links.php?saolid=69te9
  ___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[kui yang]

gdal 1.8.0 has problem. when the filename is in Chinese
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Thomas Gratier]

Hello,

First time I come to this list.

I have a problem using FwTools 2.4.7 with VRT and CSV
I want to transform french geonames file to PostgreSQL with the geometry build

What I do to get this issue ?

Creation of a VRT called geonames.vrt containing


OGRVRTDataSource
OGRVRTLayer name=geonames
SrcDataSourcegeonames.csv/SrcDataSource
GeometryTypewkbPoint/GeometryType
LayerSRSEPSG:4326/LayerSRS
GeometryField encoding=PointFromColumns x=latitude y=longitude/
/OGRVRTLayer
/OGRVRTDataSource

or with attributes

OGRVRTDataSource
OGRVRTLayer name=geonames
SrcDataSourcegeonames.csv/SrcDataSource
Field name=latitude type=Real /
Field name=longitude type=Real /
Field name=geonameid type=Integer /
Field name=name type=String width=200 /
Field name=asciiname type=String width=200 /
Field name=alternatenames type=String width=5000 /
Field name=feature class type=String width=10 /
Field name=feature code type=String width=10 /
Field name=country code type=String width=2 /
Field name=cc2 type=String width=60 /
Field name=admin1 code type=String width=20 /
Field name=admin2 code type=String width=80 /
Field name=admin3 code type=String width=20 /
Field name=admin4 code type=String width=20 /
Field name=population type=Integer /
Field name=elevation type=Integer /
Field name=gtopo30 type=Integer /
Field name=timezone type=String width=100 /
Field name=modification date type=Date /
GeometryTypewkbPoint/GeometryType
LayerSRSEPSG:4326/LayerSRS
GeometryField encoding=PointFromColumns x=latitude y=longitude/
/OGRVRTLayer
/OGRVRTDataSource









Creation of a geonames.csv file with this structure ( file around 3 Mo
for all the database on http://raster.analysesig.net/geonames.zip)


latitude,longitude,geonameid,name,asciiname,alternatenames,feature
class,feature code,country code,cc2,admin1 code,admin2
code,admin3 code,admin4
code,population,elevation,gtopo30,timezone,modification
date
45.81667,1.76667,2967103,Peyrat-le-Château,Peyrat-le-Chateau,Peyrat,P,PPL,FR,,B1,87,872,87117,1140,,463,Europe/Paris,2006-04-29
43.16667,-0.91667,2967104,Soule,Soule,La
Soule,Soule,Xiberuea,Xiberüa,Zuberoa,L,RGN,FR,,97,64,,,0,,182,Europe/Paris,2006-04-07
45.16667,-0.68333,2967105,Île Boucheau,Ile Boucheau,Ile
Bouchaud,Ile du Grand Fagnard,Ile du Grand Fagnarg,Île Bouchaud,Île du
Grand Fagnard,Île du Grand
Fagnarg,T,ISLX,FR,FR,970,,-,Europe/Paris,1996-03-21
44.81667,6.11667,2967106,Val Godemard,Val Godemard,Val
Godemar,T,VAL,FR,FR,B8,05,,,0,,1054,Europe/Paris,1994-01-06
47.5,3.81667,2967107,Domecy-sur-le-Vault,Domecy-sur-le-Vault,Domecy,P,PPL,FR,,A1,89,892,89146,107,,303,Europe/Paris,2006-04-30
45.1,-0.7,2967108,Blaye,Blaye,Blaia,Blaye-et-Sainte-Luce,P,PPLA3,FR,,97,33,331,33058,5277,,-,Europe/Paris,2010-04-19
50.8,2.4,2967109,Zuytpeene,Zuytpeene,,P,PPL,FR,,B4,59,594,59669,483,,30,Europe/Paris,2006-04-30

Finally launching of

--
ogr2ogr -f PostgreSQL PG:host=localhost user=username
dbname=mydatabase password=my_password geonames.vrt

--
It stops in the middle of the import around line 51000 depending of
the use of the first or the second VRT example

I supposed I miss something but what ?

Thanks for any help

Regards

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

Re: [gdal-dev] (no subject)

[Even Rouault]

Thomas,

The CSV file is corrupted at line 51092.

43.21667,3.1,3018296,Fleury,Fleury,Fleury Bouches de 
l'Aude,Fleury-d'Aude,Perignan,Pérignan,P,PPL,FR,,A9,11,113,11145,2713,,60,Europe/Paris,2007-09-08

Notice the extra double-quote character between Fleury and Bouches de l'Aude. 
If you remove it, the process should go just fine

I also thing you should change the role of x and y in your VRT definition : x 
should be longitude and y latitude.

Best regards,

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

[gdal-dev] (no subject)

[Discourse Maps]

I am simply trying to open a raster image, manipulate in NumPy and then spit
it back out in GDAL.

I am having trouble setting null data...
I tried: ds.GetRasterBand(1).SetNoDataValue( - )
but nothing seems to work.

Any help is appreciated!


Below is my code:


#! /usr/bin/env python 

import os, sys
use_numeric = True

try:
from osgeo import ogr, gdal
from osgeo.gdalconst import *
import numpy
os.chdir(r'L:\users\gv\numPyGDAL\ospy_data4\ospy_data4')
use_numeric = False

except ImportError:
import ogr, gdal
from gdalconst import *
import Numeric
os.chdir(r'L:\users\gv\numPyGDAL\ospy_data4\ospy_data4')


# register all of the drivers
gdal.AllRegister()

# open the image
ds = gdal.Open('sequest', GA_ReadOnly)
if ds is None:
print 'Could not open image'
sys.exit(1)

band = ds.GetRasterBand(1))# 1-based index
# read data and add the value to the string
data = band.ReadAsArray()
print data
dims = data.shape 
nx=dims[1] 
ny=dims[0]
nxarray=numpy.zeros(nx)

dst_format = 'HFA' 
dst_datatype = gdal.GDT_Int16
dst_options = ['COMPRESS=LZW'] 
dst_file='sequest_final.img' 
dst_xsize = nx 
dst_ysize = ny 
dst_nbands = 1 

geoTransform = ds.GetGeoTransform()
proj = ds.GetProjection()

driver = gdal.GetDriverByName( dst_format ) 
dst_ds = driver.Create(dst_file, dst_xsize, dst_ysize, dst_nbands, 
dst_datatype, dst_options)

for j in range(ny):
nxarray = data[j,:]
nxarray.shape = (1,nx)
print nxarray
dst_ds.GetRasterBand(1).WriteArray(nxarray,0,j)
dst_ds.SetGeoTransform(geoTransform)
dst_ds.SetProjection(proj)
dst_ds = None 
_
Hotmail: Free, trusted and rich email service.
http://clk.atdmt.com/GBL/go/201469228/direct/01/___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Frank Warmerdam]

Discourse Maps wrote:

I am simply trying to open a raster image, manipulate in NumPy and then spit
it back out in GDAL.

I am having trouble setting null data...
I tried: ds.GetRasterBand(1).SetNoDataValue( - )
but nothing seems to work.

...

Dear Discourse Maps,

Presumably you should be setting the nodata value on the
first band of the dst_ds, not the source dataset.  But
I don't see any SetNoDataValue() call at all in your script.

I would suggest adding it right after the Create() call.

I would also note that the HFA driver does not have a COMPRESS=LZW
creation option, though you can use [ 'COMPRESSED=YES' ] to use
the grid compression scheme supported by the format.


dst_format = 'HFA'
dst_datatype = gdal.GDT_Int16
dst_options = ['COMPRESS=LZW']
dst_file='sequest_final.img'
dst_xsize = nx
dst_ysize = ny
dst_nbands = 1

geoTransform = ds.GetGeoTransform()
proj = ds.GetProjection()

driver = gdal.GetDriverByName( dst_format )
dst_ds = driver.Create(dst_file, dst_xsize, dst_ysize, dst_nbands,
dst_datatype, dst_options)

for j in range(ny):
nxarray = data[j,:]
nxarray.shape = (1,nx)
print nxarray
dst_ds.GetRasterBand(1).WriteArray(nxarray,0,j)
dst_ds.SetGeoTransform(geoTransform)
dst_ds.SetProjection(proj)
dst_ds = None


Your email client has reformatted the script so it is hard
to see the indentation.  But I would advise strongly against
doing SetGeoTransform() and SetProjection() many times on the
same dataset.

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] (no subject)

[Even Rouault]

Roger,

this is weird indeed as I use daily MSYS on trunk... I'd bet that you have  an
installed version of cpl_string.h in /usr/local/include that corresponds to the
1.6.3 version that didn't have the CPLIsUTF8 and CPLForceToASCII symbols. If
you've done a make install when building 1.6.3 before, it would have gone there.
And as -I/usr/local/include is before -I/home/s1155/gdal_svn/gdal/port on the
GCC command line, the old version would be picked up. I'd recommand you to use
--prefix when you configure to be able easily to switch between different
versions to avoid installing stuff in system directories.

Best regards,

Even

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

[gdal-dev] (no subject)

[Even Rouault]

Hi Tamas,

2 issues directly related to the big raster dimensions :

* Currently GDAL allocates an array of pointer to raster blocks of size
nBlocksPerRow x nBlocksPerColumn where nBlocksPerRow =
(nRasterXSize+nBlockXSize-1) / nBlockXSize and nBlocksPerColumn =
(nRasterYSize+nBlockYSize-1) / nBlockYSize.
When TileLevel = 20, nBlocksPerRow = nBlocksPerColumn = 1 048 576. This is too
big. Currently, we then try to use a second level of block array
(sub-blocking) as soon as nBlocksPerRow  SUBBLOCK_SIZE/2 (SUBBLOCK_SIZE=64).
This leads to allocating nSubBlocksPerRow x nSubBlocksPerColumn where
nSubBlocksPerRow = (nBlocksPerRow + SUBBLOCK_SIZE + 1)/SUBBLOCK_SIZE and 
nSubBlocksPerColumn = (nBlocksPerColumn + SUBBLOCK_SIZE + 1)/SUBBLOCK_SIZE.
In that case, nSubBlocksPerRow = nSubBlocksPerColumn = 16384, and
nSubBlocksPerRow x nSubBlocksPerColumn = 268 435 456 pointers, so 1 GB on a
32bit machine or 2 GB on a 64bit machine.

-- So we would need a three-level block cache. This would be feasable
(implementation detail of gdalrasterband.cpp).

Or I'm thinking that we could completely remove the array and use instead a hash
map that would map a block number to its raster block. Could cause a (small?)
slowdown of course but would be simpler to implement and maintain than a 3 level
array.

* At TileLevel = 23, nRasterXSize = 2 147 483 648, which is ahem just one byte
more than the largest signed int, so = -2147483648
-- The proper solution would be to promote nRasterXSize/nRasterYSize to be of
type GIntBig (64bit). We would also need to update the IReadBlock(),
IWriteBlock(), IRasterIO() interfaces (and probably many others) to use GIntBig.
So this is both a large API and ABI breakage, that IMHO looks more like a GDAL
2.0 project... An alternative would be to have both the existing 'int'
variables/interfaces and introduce new GIntBIG variables/interfaces, and migrate
GDAL internals (gcore/* alg/*) and the relevant drivers such as WMS to use the
GIntBig version. Could easily turn to be very messy in the end...

I cannot think of an easy workaround. An idea would be to split the WMS extent
into 4 pieces to go under the 2 147 486 647 pixel limit, but this is not really
doable without a small hack in the TMS minidriver as you would need to add an
offset to the x and y tile number requested to the server.

Best regards,

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

[gdal-dev] (no subject)

[Martin Geng]

I have a city map in Cassini_Soldner projection with a separate file containing 
the Soldner coordinates of the 4 corners.
I want to reproject it to WGS84 in order to render it e.g. with SharpMap and 
display some GPS positions.
I have tried it with some different parameter settings of gdalwarp but in every 
case the resulting image is distorted (the size of the generated tif is 
decreased to 15141, 6943).
Even if with the parameter -ts 15118 11339 the original size is adjusted 
SharpMap or QGis will display it distorted (IrfanView without distortion).

Is there another way to get a non distorted image or what did I wrong?

Martin

..\FWTools2.4.2\bingdalinfo pic.tif
Driver: GTiff/GeoTIFF
Files: pic.tif
Size is 15118, 11339
Coordinate System is `'
Metadata:
  TIFFTAG_SOFTWARE=IrfanView
  TIFFTAG_XRESOLUTION=600
  TIFFTAG_YRESOLUTION=600
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (0.0,0.0)
Lower Left  (0.0,11339.0)
Upper Right (15118.0,0.0)
Lower Right (15118.0,11339.0)
Center  ( 7559.0, 5669.5)
Band 1 Block=15118x1 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)
0: 0,0,0,255
1: 1,1,1,255


gcp mapXmapYpixelX  pixelY

upper left  27200   17200   0   0
upper right 30400   17200   15118   0
lower right 30400   14800   15118   -11339
lower left  27200   14800   0   -11339


Gcp mapXmapYpixelX  pixelY
---
upper left 13.4370322252183 52.4818154665202 0 0
upper right 13.4841338099227 52.4818815506509 15118 0
lower right 13.4842034283538 52.4603136448411 15118 -11339
lower left 13.4371248633278 52.4602476117909 0 -11339

First I use gdal_translate to add the corner pixel coordinates in WGS84 to the 
file as ground control points:

..\FWTools2.4.2\bingdal_translate pic.tif pic_gcp.tif
-gcp 0 0 13.4370322252183 52.4818154665202
-gcp 15118 0 13.4841338099227 52.4818815506509
-gcp 15118 11339 13.4842034283538 52.4603136448411
-gcp 0 11339 13.4371248633278 52.4602476117909

..\FWTools2.4.2\bingdalinfo pic_gcp.tif
Driver: GTiff/GeoTIFF
Files: pic_gcp.tif
Size is 15118, 11339
Coordinate System is `'
GCP Projection =
GCP[  0]: Id=1, Info= (0,0) - (13.4370322252183,52.4818154665202,0)
GCP[  1]: Id=2, Info= (15118,0) - (13.4841338099227,52.4818815506509,0)
GCP[  2]: Id=3, Info= (15118,11339) - (13.4842034283538,52.4603136448411,0)
GCP[  3]: Id=4, Info= (0,11339) - (13.4371248633278,52.4602476117909,0)
Metadata:
  TIFFTAG_SOFTWARE=IrfanView
  TIFFTAG_XRESOLUTION=600
  TIFFTAG_YRESOLUTION=600
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (0.0,0.0)
Lower Left  (0.0,11339.0)
Upper Right (15118.0,0.0)
Lower Right (15118.0,11339.0)
Center  ( 7559.0, 5669.5)
Band 1 Block=15118x1 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)
0: 0,0,0,255
1: 1,1,1,255

To perform the transformation I call gdalwarp
note:
-s_srs: By giving the ground control points in WGS84 coordinates the original 
Soldner projection has been discarded.
-tr: sets target resolution equal in x- and y- directions to get a non 
distorted picture

gdalwarp -s_srs EPSG:4326 -t_srs EPSG:4326 -tr 0.03115596289 
0.03115596289 pic_gcp.tif result2.tif
Copying color table from pic_gcp.tif to new file.
Creating output file that is 15141P x 6943L.
Processing input file pic_gcp.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.

..\FWTools2.4.2\bingdalinfo result2.tif
Driver: GTiff/GeoTIFF
Files: result2.tif
Size is 15141, 6943
Coordinate System is:
GEOGCS[WGS 84,
DATUM[WGS_1984,
SPHEROID[WGS 84,6378137,298.257223563,
AUTHORITY[EPSG,7030]],
AUTHORITY[EPSG,6326]],
PRIMEM[Greenwich,0],
UNIT[degree,0.0174532925199433],
AUTHORITY[EPSG,4326]]
Origin = (13.437037980137905,52.481881537880760)
Pixel Size = (0.03115596289,-0.03115596289)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  13.4370380,  52.4818815) ( 13d26'13.34E, 52d28'54.77N)
Lower Left  (  13.4370380,  52.4602500) ( 13d26'13.34E, 52d27'36.90N)
Upper Right (  13.4842112,  52.4818815) ( 13d29'3.16E, 52d28'54.77N)
Lower Right (  13.4842112,  52.4602500) ( 13d29'3.16E, 52d27'36.90N)
Center  (  13.4606246,  52.4710657) ( 13d27'38.25E, 52d28'15.84N)
Band 1 Block=15141x1 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)
0: 0,0,0,255
1: 1,1,1,255


-- 
Jetzt kostenlos herunterladen: Internet Explorer 8 und Mozilla Firefox 3 -
sicherer, schneller und einfacher! http://portal.gmx.net/de/go/atbrowser
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Riki Tiki]

@Frank and Hermann: Thanks a lot!

My other question is as follows:

1) gdal_translate -of AAIGrid a.tif a.asc

it creates the following: a.asc (ASC File), a.asc.aux (XML Document), and a 
(IDL project file)...

now to do the reverse i do the following:

2) gdal_translate -of GTiff a.asc a.tif

it creates: a.tif (TIF Image)

but when i open a.tif with ArcGIS it shows that its statistics are not 
calculated...

Can you tell me how to create the 2) with statistics already computed? So that 
when i open it in ArcGis for example I do not have to ask ArcGIS to compute the 
statistics? Thanks a lot!


===

@Joaquim:

when the AAIGrid is created its format is of the following form:

ncols250
nrows398
xllcornerX.Y
yllcornerZ.W
cellsize 30.863820105342
NODATA_value255
 255 255 255 255 255 255 255 255.


Now, I will be changing the for example this to the following

ncols250
nrows398
xllcornerX.Y
yllcornerZ.W
cellsize 30.863820105342
NODATA_value255
 255 255 2 255 1 255 0 255.

(this is only to show my intention...)


=

Best wishes all of you,
Riki





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

Re: [gdal-dev] (no subject)

[Frank Warmerdam]

2009/8/10 Riki Tiki mongoose_r...@mail.ru:


 Hi All,

 I was trying to write a c# program to copy the existing raster file to a new 
 raster file so that at the end i both files while displayed on tools like 
 arcmap will be exactly the same... i started as follows


 string src_file = @C:\raster1;
 Gdal.AllRegister();
 Dataset src_ds  = Gdal.Open(src_file, Access.GA_ReadOnly);

 Driver dr = Gdal.GetDriverByName(GTiff);
 string dest_file = @C:\raster_copy;
 Dataset dest_ds = dr.Create(dest_file, src_ds.RasterXSize, src_ds.RasterYSize,
                            1, DataType.GDT_Float32, null);



 now i do not know how to continue...? plz, can you help me?

 http://trac.osgeo.org/gdal/browser/trunk/gdal/swig/csharp/apps/GDALCreateCopy.cs

 i have already seen above one but that creates two different images when i 
 open in arcmap.

Riki,

In order to do better than GDALCreateCopy.cs you will need to
analyse why the source and destination images are different
in arcmap.  Is it a georeferencing issues?  Pixel types?

What is the source image format?  Perhaps it is an
arc binary grid?  These can have auxilary information
that does not appear in formats like GeoTIFF.   It might
be helpful to use gdalinfo on the source and destination
and see if there are any apparent differences.

If you want to see how to do something similar to
CreateCopy() using Create and direct methods you
can look through the default implementation of CreateCopy.

  http://trac.osgeo.org/gdal/browser/trunk/gdal/gcore/gdaldriver.cpp#L334

Good luck,
-- 
---+--
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] (no subject)

[Riki Tiki]

Hi All,

I was trying to write a c# program to copy the existing raster file to a new 
raster file so that at the end i both files while displayed on tools like 
arcmap will be exactly the same... i started as follows


string src_file = @C:\raster1;
Gdal.AllRegister();
Dataset src_ds  = Gdal.Open(src_file, Access.GA_ReadOnly);

Driver dr = Gdal.GetDriverByName(GTiff);
string dest_file = @C:\raster_copy;
Dataset dest_ds = dr.Create(dest_file, src_ds.RasterXSize, src_ds.RasterYSize,
1, DataType.GDT_Float32, null);



now i do not know how to continue...? plz, can you help me?

http://trac.osgeo.org/gdal/browser/trunk/gdal/swig/csharp/apps/GDALCreateCopy.cs

i have already seen above one but that creates two different images when i open 
in arcmap.

thanks a lot!
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Le minh NHAT]

Thank you

 

Le MINH NHAT  Dr. Eng. 
Disaster Management Center (DMC) 
Central Committe for Flood and Storm Control (CCFSC) 
Ministry of Agriculture and Rural Development (MARD) 
A4 No. 2 Ngoc Ha Str., Hanoi, Vietnam
Tel.: (84-4) 3 733 5686 / 3733 5805
Fax.: (84-4) 3733 6647 / 3733 5701  
Cel.: 0123 6073 898 
Email: minhn...@ccfsc.org.vn / nhatky...@gmail.com

 

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

[gdal-dev] (no subject)

[Aaron Brown]

Hi all,

I have to convert DTED Level 1 data into DTED Level 0 data, and was just 
wondering if anyone would be able to help?

Regards

_
View photos of singles in your area Click Here
http://dating.ninemsn.com.au/search/search.aspx?exec=gotp=qgc=2tr=1lage=18uage=55cl=14sl=0dist=50po=1do=2trackingid=1046138r2s=1_t=773166090_r=WLM_EndText___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Frank Warmerdam]

Sjur Kolberg wrote:
 
Hello, GDAL list;
 
Can anyone suggest what prevents the following GDALSuggestedWarpOutput 
from calculating an approximate geotransform? The resulting suggested 
GeoTransform is just a copy of the original GeoTransform.
 
Version is GDAL 1.6.0beta1 under Windows XP, I use Visual Studio 2005, 
this is a MFC application. poDataset is a 9900 band GRIB 1 data set 
from ECMWF. It appears to be read in OK, I can extract 
plausible metadata both from the dataset and from each band. 
 
   SourceWKT is (CString)

GEOGCS[Coordinate System imported from GRIB file,DATUM[unknown,
SPHEROID[Sphere,6367470,0]],PRIMEM[Greenwich,0],UNIT[degree,0.0174532925199433]]
   SourceWKT is from poDataset-GetProjectionRef()
 
   TargetWKT is (CString)

PROJCS[UTM 32 North (WGS84),GEOGCS[WGS 84,DATUM[WGS_1984,
SPHEROID[WGS 84,6378137,298.257223563,AUTHORITY[EPSG,7030]],
TOWGS84[0,0,0,0,0,0,0],AUTHORITY[EPSG,6326]],PRIMEM[Greenwich,
0,AUTHORITY[EPSG,8901]],UNIT[degree,0.0174532925199433,
AUTHORITY[EPSG,9108]],AXIS[Lat,NORTH],AXIS[Long,EAST],
AUTHORITY[EPSG,4326]],PROJECTION[Transverse_Mercator],
PARAMETER[latitude_of_origin,0],PARAMETER[central_meridian,9],
PARAMETER[scale_factor,0.9996],PARAMETER[false_easting,50],
PARAMETER[false_northing,0],UNIT[Meter,1]]
   TargetWKT is from 
ogrspatref.SetProjCS(msg);// CString msg is just a name

ogrspatref.SetWellKnownGeogCS( WGS84 );
if (ogrspatref.SetUTM( zone, i ) != OGRERR_NONE)  // zone==32, i==1 for 
northern hem.
 
 
And the code goes:



void *hTransformArg = GDALCreateGenImgProjTransformer( 
poDataset, sourceWKT, NULL, targetWKT, FALSE, 0, 1 );


CPLAssert( hTransformArg != NULL );

if (GDALSuggestedWarpOutput(pSrcDS, GDALGenImgProjTransform, 
hTransformArg, SugGeoTransform, nPixels, nLines )  != CE_None )


Sjur,

I tried this (test.wkt is the UTM WKT) and it seems to work plausibly:

testepsg -t 'GEOGCS[Coordinate System imported from GRIB 
file,DATUM[unknown,SPHEROID[Sphere,6367470,0]],PRIMEM[Greenwich,0],UNIT[degree,0.0174532925199433]]' 
test.wkt 4 65

OGRCT: Source: +proj=longlat +a=6367470 +b=6367470 +no_defs
OGRCT: Target: +proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
(4.00,65.00,0.00) - (264409.202311,7217779.091236,0.00)

So, the reprojection operation can apparently work.  The behavior you
are seeing appears to be congruent with reprojection doing nothing.
Is there any possibility that the PROJ library is not built into your
GDAL?

I see no other reason things should not work as you have constructed
them.  All I can see is to either debug through and see what is going
wrong, or to submit a ticket with info on how to reproduce the problem.

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] (no subject) GDALSuggestedWarpOutput() only returns the unaltered source geotransform

[Sjur Kolberg]

Yes, thanks, Frank. That was it.

Or even easier, I simply copied a proj.dll from another location into the 
directory where my program (and gdal.dll) resides.

Is it so that gdal looks for proj.dll without being linked with proj.lib, and 
without complaining if it's not found?

I'll build a gdal with debug info, allowing me to track.

Best regards and thanks,

Sjur :-)


 

-Original Message-
From: Frank Warmerdam [mailto:warmer...@pobox.com] 
Sent: Friday, May 15, 2009 4:34 PM
To: Sjur Kolberg
Cc: 'gdal-dev@lists.osgeo.org'
Subject: Re: [gdal-dev] (no subject)


Sjur,

I tried this (test.wkt is the UTM WKT) and it seems to work plausibly:

testepsg -t 'GEOGCS[Coordinate System imported from GRIB 
file,DATUM[unknown,SPHEROID[Sphere,6367470,0]],PRIMEM[Greenwich,0],UNIT[degree,0.0174532925199433]]'
 
test.wkt 4 65
OGRCT: Source: +proj=longlat +a=6367470 +b=6367470 +no_defs
OGRCT: Target: +proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
(4.00,65.00,0.00) - (264409.202311,7217779.091236,0.00)

So, the reprojection operation can apparently work.  The behavior you are 
seeing appears to be congruent with reprojection doing nothing.
Is there any possibility that the PROJ library is not built into your GDAL?

I see no other reason things should not work as you have constructed them.  All 
I can see is to either debug through and see what is going wrong, or to submit 
a ticket with info on how to reproduce the problem.

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



Sjur Kolberg wrote:
  
 Hello, GDAL list;
  
 Can anyone suggest what prevents the following GDALSuggestedWarpOutput 
 from calculating an approximate geotransform? The resulting suggested 
 GeoTransform is just a copy of the original GeoTransform.
  
 Version is GDAL 1.6.0beta1 under Windows XP, I use Visual Studio 2005, 
 this is a MFC application. poDataset is a 9900 band GRIB 1 data set 
 from ECMWF. It appears to be read in OK, I can extract plausible 
 metadata both from the dataset and from each band.
  
SourceWKT is (CString)
 GEOGCS[Coordinate System imported from GRIB file,DATUM[unknown, 
 SPHEROID[Sphere,6367470,0]],PRIMEM[Greenwich,0],UNIT[degree,0.0174532925199433]]
SourceWKT is from poDataset-GetProjectionRef()
  
TargetWKT is (CString)
 PROJCS[UTM 32 North (WGS84),GEOGCS[WGS 84,DATUM[WGS_1984, 
 SPHEROID[WGS 84,6378137,298.257223563,AUTHORITY[EPSG,7030]],
 TOWGS84[0,0,0,0,0,0,0],AUTHORITY[EPSG,6326]],PRIMEM[Greenwich,
 0,AUTHORITY[EPSG,8901]],UNIT[degree,0.0174532925199433,
 AUTHORITY[EPSG,9108]],AXIS[Lat,NORTH],AXIS[Long,EAST],
 AUTHORITY[EPSG,4326]],PROJECTION[Transverse_Mercator],
 PARAMETER[latitude_of_origin,0],PARAMETER[central_meridian,9],
 PARAMETER[scale_factor,0.9996],PARAMETER[false_easting,50],
 PARAMETER[false_northing,0],UNIT[Meter,1]]
TargetWKT is from 
 ogrspatref.SetProjCS(msg);// CString msg is just a 
 name
 ogrspatref.SetWellKnownGeogCS( WGS84 ); if (ogrspatref.SetUTM( zone, 
 i ) != OGRERR_NONE)  // zone==32, i==1 for northern hem.
  
  
 And the code goes:
 
 
 void *hTransformArg = GDALCreateGenImgProjTransformer( poDataset, 
 sourceWKT, NULL, targetWKT, FALSE, 0, 1 );
 
 CPLAssert( hTransformArg != NULL );
 
 if (GDALSuggestedWarpOutput(pSrcDS, GDALGenImgProjTransform, 
 hTransformArg, SugGeoTransform, nPixels, nLines )  != CE_None )
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Sjur Kolberg]

 
Hello, GDAL list;
 
Can anyone suggest what prevents the following GDALSuggestedWarpOutput 
from calculating an approximate geotransform? The resulting suggested 
GeoTransform is just a copy of the original GeoTransform.
 
Version is GDAL 1.6.0beta1 under Windows XP, I use Visual Studio 2005, 
this is a MFC application. poDataset is a 9900 band GRIB 1 data set 
from ECMWF. It appears to be read in OK, I can extract 
plausible metadata both from the dataset and from each band. 
 
   SourceWKT is (CString)
GEOGCS[Coordinate System imported from GRIB file,DATUM[unknown,
SPHEROID[Sphere,6367470,0]],PRIMEM[Greenwich,0],UNIT[degree,0.0174532925199433]]
   SourceWKT is from poDataset-GetProjectionRef()
 
   TargetWKT is (CString)
PROJCS[UTM 32 North (WGS84),GEOGCS[WGS 84,DATUM[WGS_1984,
SPHEROID[WGS 84,6378137,298.257223563,AUTHORITY[EPSG,7030]],
TOWGS84[0,0,0,0,0,0,0],AUTHORITY[EPSG,6326]],PRIMEM[Greenwich,
0,AUTHORITY[EPSG,8901]],UNIT[degree,0.0174532925199433,
AUTHORITY[EPSG,9108]],AXIS[Lat,NORTH],AXIS[Long,EAST],
AUTHORITY[EPSG,4326]],PROJECTION[Transverse_Mercator],
PARAMETER[latitude_of_origin,0],PARAMETER[central_meridian,9],
PARAMETER[scale_factor,0.9996],PARAMETER[false_easting,50],
PARAMETER[false_northing,0],UNIT[Meter,1]]
   TargetWKT is from 
ogrspatref.SetProjCS(msg);// CString msg is just a name
ogrspatref.SetWellKnownGeogCS( WGS84 );
if (ogrspatref.SetUTM( zone, i ) != OGRERR_NONE)  // zone==32, i==1 for 
northern hem.
 
 
And the code goes:


void *hTransformArg = GDALCreateGenImgProjTransformer( 
poDataset, sourceWKT, NULL, targetWKT, FALSE, 0, 1 );

CPLAssert( hTransformArg != NULL );

if (GDALSuggestedWarpOutput(pSrcDS, GDALGenImgProjTransform, 
hTransformArg, SugGeoTransform, nPixels, nLines )  != CE_None )

throw new GeneralExcpt(GDALSuggestedWarpOutput failed);




No exception is thown or other error, so everything 
seems OK, but the resulting SugGeoTransform is:

  [0] 4. double
  [1] 0.25000 double
  [2] 0.0 double
  [3] 65.000 double
  [4] 0.0 double
  [5] -0.25000 double
and nPixels, nLines are 33, 45 
This happens to be exactly the GeoTransform of poDataSet, in degrees.
 
 
I guess I'm missing something obvious, but I've spent some evenings now.
 
Thanks for any help.
 
Sjur K :-)
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] (no subject)

[Even Rouault]

I've not reviewed closely your code, but I'd suggest you to study the source 
code for gdalwarp, that probably contains the missing bit.

I think you should probably try using GDALCreateGenImgProjTransformer2() 
instead of GDALCreateGenImgProjTransformer(). That's what gdalwarp.cpp 
does...

http://trac.osgeo.org/gdal/browser/trunk/gdal/apps/gdalwarp.cpp


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

[gdal-dev] (no subject)

[Kolberg Sjur A]

Hello, list,
 
I have a GDAL warp operation creating only missing-values in the target
map, 
suggesting that the area I want to warp to is not covered by the
original image.
 
Trying to to trace the error, I find a strange behaviour of the 
GDALSuggestedWarpOutput() call: The resulting target GeoTransform array
has 
(almost) exactly the same values as the source GeoTransform, even though
the 
projections are different.
 
Using GDAL 1.6.0b1 under Windows, I have the following code
 
void *hTransformArg = GDALCreateGenImgProjTransformer( pSrc, 
 
(LPCTSTR) inWKT, 
 
NULL,
 
(LPCTSTR) outWKT, 
 
FALSE, 0, 1 );
CPLAssert( hTransformArg != NULL );
if (GDALSuggestedWarpOutput( pSrc, // GDALDataset
GDALGenImgProjTransform,
// function
hTransformArg, // void
pointer
adfDstGeoTransform, //
double *
nPixels, // int *
nLines ) // int *
!= CE_None ) 
{
AfxMessageBox(GDALSuggestedWarpOutput failed);
return NULL;
}
 
Prior to the GDALCreateGenImgProjTransformer() call, inWKT is
 
PROJCS[unnamed,GEOGCS[Normal Sphere (r=6370997),
DATUM[unknown,SPHEROID[sphere,6370997,0]],
PRIMEM[Greenwich,0],UNIT[degree,0.0174532925199433]],
PROJECTION[Polar_Stereographic],
PARAMETER[latitude_of_origin,60],
PARAMETER[central_meridian,0],
PARAMETER[scale_factor,1],
PARAMETER[false_easting,0],
PARAMETER[false_northing,0],
UNIT[Meter,1]]
 
whereas outWKT is
 
PROJCS[UTM 32 North (WGS84),PROJECTION[Transverse_Mercator],
PARAMETER[latitude_of_origin,0],
PARAMETER[central_meridian,9],
PARAMETER[scale_factor,0.9996],
PARAMETER[false_easting,50],
PARAMETER[false_northing,0],UNIT[Meter,1]]
 
The source GeoTransform matrix is:
- adfGeoTransform 0x0e3eaf58 double [6]
[0] 255106.000 double
[1] 1000.0 double
[2] 0.0 double
[3] -2511133.00 double
[4] 0.0 double
[5] -1000.0 double
 
while after the GDALSuggestedWarpOutput call, the destination
GeoTransform is:
- adfDstGeoTransform,6 0x042bdd18 double [6]
[0] 255106.000 double
[1] 999.89 double
[2] 0.0 double
[3] -2511133.00 double
[4] 0.0 double
[5] -999.89 double
which is a mere copy of the source numbers.
 
As I afterwards specify the transform I want for the warped image:
- adfDstGeoTransform,6 0x042bdd18 double [6]
[0] 33.000 double
[1] 1000.0 double
[2] 0.0 double
[3] 728.00 double
[4] 0.0 double
[5] -1000.0 double
 
...my specified UL-Y [3] (728) is pretty far from the suggested
-2511133.
Can this be a cause to my resulting image being empty?
This has been working up to recently, but I have been working on this
code as well as changing the GDAL version.
 
Any ideas?
Best regards,
 
Sjur K :-)
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] (no subject)

[Sarah Mulholland]

I need to write multiple images to a NITF file.  I notice on the gdal NITF page 
that gdal only writes one image or raster of data to a file.  On another page, 
I see that multiple images in a NITF can be read as subdatasets.  Is the 
documentation up to date that subdatasets in a NITF file are a read-only 
capability? In there any plan to allow the writing of subdatasets?

What's the usual approach for writing multiple images to a NITF from C and 
python?  If not gdal, is there some other package that people use?  I looked on 
the JITC web site, and it seems like most certified packages are parts of large 
systems, not just NITF libraries.

Any suggestions?  Thanks in advance.

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

[gdal-dev] (no subject)

[sara lucca]

Hi,

I'm an Italian university student and I'm working on LiDAR scanlations with 
GRASS and the data I'm using were in LAS format. 
One of the problem I've encountered is the conversion of the LAS file in file 
that GRASS can manage (ASCII).


I
saw that this format is not present in your supported format and I
would like to know if there is any project aimed at adding this format
in your vector format list.


Sorry for my bad english and thanks for the attention
Sara

_
Screensaver e sfondi… Scaricali dallo Spaces di Doretta!
https://cid-97b0d995232b954a.skydrive.live.com/home.aspx___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev