Re: [gdal-dev] Linux installation

2011-11-10 Thread christian . mueller
Hi Steve, one of my customers is working with SLES11 on IBM z/Series and p/Series (NO Intel architecture). I have to administer some of these logical partitions. First, what hardware do you plan to use. Intel/AMD makes life easier. Second, SLES 11 is an enterprise edition with support you ha

Re: [gdal-dev] How would I build a custom pyramid ?

2011-08-22 Thread christian . mueller
nal data (0.17 and 0.33 are smaller than the original 0.5). So I would need to increase (instead of decrease) the number of pixels to create the lowest pyramid level. Or is it useless considering the way geoserver handles pyramid data ? Alain - Mail original ----- De: "chri

Re: [gdal-dev] How would I build a custom pyramid ?

2011-08-22 Thread christian . mueller
Zitat von Alain Camus : Hello list, I have a set of tiles at resolution 0.5 m/pxl. Based on this data, I'd like to create a pyramid that has zoom levels 0.17, 0.33, 0.67 and 1.7 (the data will only be viewed in geoserver at these levels, so I guess it will be faster if

Re: [gdal-dev] image quality with pyramids (with gdal_retile)

2011-03-15 Thread christian . mueller
Quoting Devrim Baris Acar : Hi, I am building a pyramid from a 650 MB geotiff file (rgb) , with the following command from fwtools 2.4.6 on a win32 box gdal_retile -levels 5 -pyramidOnly -targetDir pyramid bathymetry. I loaded these image files to geoserver and see that when I zoomed to the hi

Re: [gdal-dev] gdal_retile error

2011-03-08 Thread christian . mueller
Perhaps there is a resource limit for the number of the files in a directory. Try using -useDirForEachRow This builds another directory hierarchy with less files per directory. Quoting Vasile Craciunescu : Christian, The permissions are OK. 1-8 subdirs are created as soon as the command i

Re: [gdal-dev] gdal_retile error

2011-03-08 Thread christian . mueller
Perhaps a problem with the file permissions ? After creating the tiles for level 0, the utility wants to create a subdir called "1" and store the tiles for level 1 in this subdir. Check this with your user. Make mkdir 1 cd 1 touch test.tif rm test.tif Quoting Vasile Craciunescu : Dear al

[gdal-dev] Re: Performance improvment for gdal_retile

2010-12-18 Thread christian . mueller
Ah, things get complicated :-(. I am working on openSUSE 11.3, 32 Bit and Ubuntu 10.10 64 Bit. Additionally, I have SUSE Linux Enterprise on IBM pSeries and IBM zSeries (Mainframe). Some times I even use gdal on AIX. I compile against gdal 1.7, no debug mode. The important fact is that you

[gdal-dev] Performance improvment for gdal_retile

2010-12-18 Thread christian . mueller
Some last comment about the test case, the current source checks test_gdal_retile_2() ds.GetRasterBand(4).Checksum() != 35 ds.GetRasterBand(4).Checksum() != 35 If I run this on a Linux Intel 32 bit configuration, it has to be ds.GetRasterBand(4).Checksum() != 283 ds.GetRasterBand(4).Checksum()

[gdal-dev] Re: Performance improvement for gdal_retile

2010-12-16 Thread christian . mueller
Hi Franks, thanks for your quick answers. The patch is here http://trac.osgeo.org/gdal/ticket/3879 The test case test_gdal_retile.py is running and I checked the man page with doxygen. Please can you give me a ping after applying the patch to trunk. Since I am a geotools/geoserver develope

[gdal-dev] Re: Performance improvement for gdal_retile

2010-12-15 Thread christian . mueller
I have seen that Even has invested some time into gdal_retile, writing a test case on trunk and making some python specific changes to the utility itself. Since gdal_retile.py was my first python program I want to say thank you for that. Comparing 1.7 to trunk, I see that gdal_retile.py dif

[gdal-dev] Performance improvement for gdal_retile

2010-12-15 Thread christian . mueller
Hi, I am the maintainer of gdal_retile.py. A user and me developed an essential performance improvement for large images. What is the best way to get the python file and the dox file into trunk and into the 1.7 branch. I have an OSGeo Account, but no svn commit access, so I could open a

Re: [gdal-dev] towgs84 for epsg3035

2010-11-22 Thread christian . mueller
Why do you need toWGS84 params for converting wgs84 -->epsg:3035. I always thought toWGS84 params are needed for the other way around. But I am not an expert :-) ASFAIK the params are all 0. I use this definition PROJCS["ETRS89 / ETRS-LAEA", GEOGCS["ETRS89", DATUM["European_Terrestrial_Referen

Re: [gdal-dev] Re: [gdal_retile]reshape ecw

2010-10-22 Thread christian . mueller
In my experience, 128x128 is too small. 256x256 or 512x512 are good choices. I have no experience with GDALBUILDVRT, cannot assist here. Quoting Tacot : OK. So I try using this size for each tiles 128*128. The treatment is very long. So is it possible with GDALBUILDVRT, to translate virtua

Re: [gdal-dev] Re: [gdal_retile]reshape ecw

2010-10-20 Thread christian . mueller
The -ps parameter specifies the pixel width and pixel height of a tile. You must multiply these values. -ps 1000 1000 means 1000x1000 = 1 million pixels for each tile As a consequence, a tile needs 1 MB .. n MB , depending on your color model. The file size of the tiles my be smaller if co

Re: [gdal-dev] [gdal_retile]reshape ecw

2010-10-20 Thread christian . mueller
The error does not occur in the gdal_retile.py python code, it occurs in the gdal C code. Your tile size is really large producing images with 100 millions of pixel. Try a smaller pixel size, each pixel needs some bytes, dependent on your color model. Quoting Tacot : Hi, I try to cut

Re: [gdal-dev] tileindex performance

2010-10-14 Thread christian . mueller
Hi guys, I am the developer of the gdal_retile. Using gdal_retile on very large images only makes sense if you plan to import your tiles into a jdbc database. Serving images from so many files is not recommended. I tested on an AIX Box, tiling an 80 GB Erdas image. The performance decreas

Re: [gdal-dev] Re: Copy only georeference between two images

2010-09-13 Thread christian . mueller
I see your image is LZW compressed, AFAIK this is the same for png images. Try the following 1) use gdal_translate to create a png file and a corresponding world file from your shifted image. 2) Open the world file with an editor and adjust the values 3) use gdal_translate to create a geotiff

Re: [gdal-dev] Re: Copy only georeference between two images

2010-09-13 Thread christian . mueller
And what is the output of the second image ? Quoting canduc17 : My case is the first one: with gdalinfo on one of my images I obtain: Driver: GTiff/GeoTIFF Files: AVN_CLS57_ALAV2A127252720_20080614_20100906.tif Size is 8505, 8378 Coordinate System is: PROJCS["WGS 84 / UTM zone 31N", GEOGCS

Re: [gdal-dev] Copy only georeference between two images

2010-09-13 Thread christian . mueller
Do you have the georeference information included in your image file (geotiff for example) or do you have an extra world file containing 6 lines of numbers. In the second case, it is easy. Quoting canduc17 : One simple question: I have two identical images, exept for their their georefere

Re: [gdal-dev] Re: Creating new mosaic overview

2010-08-19 Thread christian . mueller
Look here, http://www.gdal.org/gdal_retile.html I think this should do your job. Cheers Quoting "Edi.Karadumi" : Yep but im creating an overview so i need a lowered resolution so i need the -ps parameter. I know that gdal_merge.py can take more than two files but the new mosaiced overview

Re: [gdal-dev] Raster/vector intersections in PostGIS WKT Raster

2010-06-30 Thread christian . mueller
No, I did not consider GDAL/Georaster. The reasons are 1) There is already a plugin (imagemosaic-jdbc) for fetching images from jdbc databases supporting Oracle Location Based Services, Postgis, DB2 spatial extender and jdbc databases without a spatial extension. I did not want to add a ne

Re: [gdal-dev] Raster/vector intersections in PostGIS WKT Raster

2010-06-30 Thread christian . mueller
Hi, the documentation for the GeoRaster geotools plugin is here http://docs.codehaus.org/display/GEOTDOC/Oracle+GeoRaster+Plugin For geoserver look here http://docs.geoserver.org/stable/en/user/data/oraclegeoraster.html Technically , I investigated int the Oralce Java API. 2 problems occurre

Re: [gdal-dev] Raster/vector intersections in PostGIS WKT Raster

2010-06-29 Thread christian . mueller
Hi, superb job. Congratulations. I am a geotools/geoserver developer and one of my jobs is fetching image tiles or the corresponding pyramid tiles from jdbc data sources. Some weeks ago, I added support for Oracle Georaster and I want to support WKT Georaster too. wktraster makes it quit

Re: [gdal-dev] Reprojecting rasters

2010-06-08 Thread christian . mueller
I am not sure here, but this looks like each tile is projected for itself. My procedure always is 1) Build a large image without tiles 2) Reproject 3) Create Tiles/pyramids Try to avoid TILED=YES and reproject within the same gdalwarp call. Quoting Rob : Hi [Apols for the X-post] I'm tr

Re: [gdal-dev] Pretty Much Off Topic: Development IDEs

2010-06-01 Thread christian . mueller
Hi, I am also a fan of working in the command line. (vi is one of my friends) I use eclipse which is very powerful. Most of the time I am doing Java, but I also used eclipse for programming python and there is support for C++/C too. Java / python really makes your life easier. The only reas

Re: [gdal-dev] Re: gdal_retile and nodata

2010-05-25 Thread christian . mueller
Hmm, jpeg is not a lossless compression, LZW is lossless. It appears that we have a problem beyond our scope. I never use jpeg and prefer LZW or png. I would go with LZW, this is the safe way. 4 GB are not dramatically. I store my tiles/pyramids in a database prepared for geoserver, there

Re: [gdal-dev] Re: gdal_retile and nodata

2010-05-21 Thread christian . mueller
Any new tile is created according to the create options "-co Name=value" you can pass in as parameter. First an image is created, larege enough to hold the tile and its neighbors from level x-1 Second , the tiles from level x-1 are painted onto this image. Third, the scaling takes place.

Re: [gdal-dev] Re: gdal_retile and nodata

2010-05-21 Thread christian . mueller
AFAIK at the time I did the gdal_retile there was no concept of aux.xml files. Anyways, the utility does no image operation itself, we will see on Monday. Quoting acangi : After removing the GDAL_PAM_ENABLED=NO, I made a gdal_translate -a_nodata 0 on the files, this generated the aux.xml fi

Re: [gdal-dev] Re: gdal_retile and nodata

2010-05-21 Thread christian . mueller
You are right, give this a try. Quoting acangi : Maybe it's related to the fact I was preventing creation of *.aux.xml files with export GDAL_PAM_ENABLED=NO Currently investigating -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/gdal-retile-and-nodata-tp5079895p508

Re: [gdal-dev] gdal_retile and nodata

2010-05-20 Thread christian . mueller
gdal_retile use gdal core functions to rescale and I see no possibility to pass a nodata value to these functions. Have your tiff files a no data value in the meta info. Quoting acangi : When creating the level x, gdal_retile interpolates the tiles of level x-1. If not all data exist in le

Re: [gdal-dev] Re: GDAL_retile slowdown at next level

2010-05-11 Thread christian . mueller
Please keep me informed, perhaps we can do an improvement by adding a command line argument for the cache size. Quoting acangi : Thanks for the quick reply. I think it makes sense. First I'll try to move all data to a computing machine with a big hard disk, then if it's still slow, I'll try

Re: [gdal-dev] GDAL_retile slowdown at next level

2010-05-11 Thread christian . mueller
Hi, I am the developer of the gdal_retile.py I think the problem here is based on the following issue. To make a pyramid with tiles in the Dimension width * hight based and already tiled image with the same tile width and height, I have to load the base tile and all 8 neighbors touching the