RE: [gdal-dev] Re: ESRI file geodatabase support
Ragi, That is a good point, as always. Unfortunately I am not able to do that, but I can help in any other way... The zigGIS driver (I think ) is for reading/writing PostGIS using ArcMap, without having ArcSDE in the middle. As far as I know it does not read File geodatabases... I think for now you would have to reverse engineer the way fgdb's are written... given the effort, the discussion on whether to invest the time on fgbd or invest it on SL is at least interesting... and even fun for some of us. Duarte De: Ragi Burhum [mailto:r...@burhum.com] Enviada: quinta-feira, 17 de Junho de 2010 21:57 Para: gdal-dev@lists.osgeo.org Cc: Duarte Carreira; Eric Wolf; Peter J Halls; Matt Wilkie Assunto: RE: [gdal-dev] Re: ESRI file geodatabase support From: Duarte Carreira dcarre...@edia.ptmailto:dcarre...@edia.pt Subject: RE: [gdal-dev] Re: ESRI file geodatabase support To: Eric Wolf ebw...@gmail.commailto:ebw...@gmail.com, Peter J Halls p.ha...@york.ac.ukmailto:p.ha...@york.ac.uk Cc: gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org, Matt Wilkie map...@gmail.commailto:map...@gmail.com Well, if SpatiaLite offers some proper benefits and disseminates through all of the FOSS world, then it may get a strong enough push even for ESRI to pick it up. It happened before... (kml?) If SL would: 1) Be as fast as shapefile in production settings, desktop and webgis 2) Offer SQL support, spatial and otherwise, also through desktop tools like QGIS 3) Allow editing while serving (even if for 1 editor only) 4) Better support in QGIS than for shapefile (take advantage of Spatial SQL, all other functionality) 5) Same for MapServer, GeoServer, gvSIG, et al. 6) Allow easy managing of rasters inside the .db file, through QGIS 7) ??more ideas/requests?? Then it would be a very, very good contender... and the ball would be kicked to the other side. And it seems we're already there for some of the listed features. Duarte Or even better, instead of waiting and complaining, we could just write a GeoDatabase OGR (ESRI) Workspace ourselves. There are already examples of working ones out there http://svn.obtusesoft.com/core/trunk/ So instead of hoping and pleading for support, someone should just sit down and write it. I started one at one point (C++), got side tracked with other things. If anyone is interested in that source code, I would be happy to share that too. - Ragi ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Re: ESRI file geodatabase support
I happened to be looking at the ESRI ArcGIS 10 online documentation and found a page of acknowledgements, which includes an acknowledgement reference to SQLite. I've no idea what to make of this - maybe someone else has? Best wishes, Peter Duarte Carreira wrote: Ragi, That is a good point, as always. Unfortunately I am not able to do that, but I can help in any other way... The zigGIS driver (I think ) is for reading/writing PostGIS using ArcMap, without having ArcSDE in the middle. As far as I know it does not read File geodatabases... I think for now you would have to reverse engineer the way fgdb's are written... given the effort, the discussion on whether to invest the time on fgbd or invest it on SL is at least interesting... and even fun for some of us. Duarte De: Ragi Burhum [mailto:r...@burhum.com] Enviada: quinta-feira, 17 de Junho de 2010 21:57 Para: gdal-dev@lists.osgeo.org Cc: Duarte Carreira; Eric Wolf; Peter J Halls; Matt Wilkie Assunto: RE: [gdal-dev] Re: ESRI file geodatabase support From: Duarte Carreira dcarre...@edia.ptmailto:dcarre...@edia.pt Subject: RE: [gdal-dev] Re: ESRI file geodatabase support To: Eric Wolf ebw...@gmail.commailto:ebw...@gmail.com, Peter J Halls p.ha...@york.ac.ukmailto:p.ha...@york.ac.uk Cc: gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org, Matt Wilkie map...@gmail.commailto:map...@gmail.com Well, if SpatiaLite offers some proper benefits and disseminates through all of the FOSS world, then it may get a strong enough push even for ESRI to pick it up. It happened before... (kml?) If SL would: 1) Be as fast as shapefile in production settings, desktop and webgis 2) Offer SQL support, spatial and otherwise, also through desktop tools like QGIS 3) Allow editing while serving (even if for 1 editor only) 4) Better support in QGIS than for shapefile (take advantage of Spatial SQL, all other functionality) 5) Same for MapServer, GeoServer, gvSIG, et al. 6) Allow easy managing of rasters inside the .db file, through QGIS 7) ??more ideas/requests?? Then it would be a very, very good contender... and the ball would be kicked to the other side. And it seems we're already there for some of the listed features. Duarte Or even better, instead of waiting and complaining, we could just write a GeoDatabase OGR (ESRI) Workspace ourselves. There are already examples of working ones out there http://svn.obtusesoft.com/core/trunk/ So instead of hoping and pleading for support, someone should just sit down and write it. I started one at one point (C++), got side tracked with other things. If anyone is interested in that source code, I would be happy to share that too. - Ragi -- Peter J Halls, GIS Advisor, University of York Telephone: 01904 433806 Fax: 01904 433740 Snail mail: Computing Service, University of York, Heslington, York YO10 5DD This message has the status of a private and personal communication ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
RE: [gdal-dev] Re: ESRI file geodatabase support
Re-reading your email it seems I misunderstood what you meant... So using the zigGIS ArcMap/PostGIS provider one could adapt it to read OGR datasources, like QGIS does. You would then access SpatiaLite and potentially all other OGR formats from ArcMap. Good idea. Duarte De: Ragi Burhum [mailto:r...@burhum.com] Enviada: quinta-feira, 17 de Junho de 2010 21:57 Para: gdal-dev@lists.osgeo.org Cc: Duarte Carreira; Eric Wolf; Peter J Halls; Matt Wilkie Assunto: RE: [gdal-dev] Re: ESRI file geodatabase support From: Duarte Carreira dcarre...@edia.ptmailto:dcarre...@edia.pt Subject: RE: [gdal-dev] Re: ESRI file geodatabase support To: Eric Wolf ebw...@gmail.commailto:ebw...@gmail.com, Peter J Halls p.ha...@york.ac.ukmailto:p.ha...@york.ac.uk Cc: gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org, Matt Wilkie map...@gmail.commailto:map...@gmail.com Well, if SpatiaLite offers some proper benefits and disseminates through all of the FOSS world, then it may get a strong enough push even for ESRI to pick it up. It happened before... (kml?) If SL would: 1) Be as fast as shapefile in production settings, desktop and webgis 2) Offer SQL support, spatial and otherwise, also through desktop tools like QGIS 3) Allow editing while serving (even if for 1 editor only) 4) Better support in QGIS than for shapefile (take advantage of Spatial SQL, all other functionality) 5) Same for MapServer, GeoServer, gvSIG, et al. 6) Allow easy managing of rasters inside the .db file, through QGIS 7) ??more ideas/requests?? Then it would be a very, very good contender... and the ball would be kicked to the other side. And it seems we're already there for some of the listed features. Duarte Or even better, instead of waiting and complaining, we could just write a GeoDatabase OGR (ESRI) Workspace ourselves. There are already examples of working ones out there http://svn.obtusesoft.com/core/trunk/ So instead of hoping and pleading for support, someone should just sit down and write it. I started one at one point (C++), got side tracked with other things. If anyone is interested in that source code, I would be happy to share that too. - Ragi ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
RE: [gdal-dev] How to represent multi-dimensional array
Chris, I have some experience trying to get ArcGIS to work well with time series satellite imagery and 4D ocean models (e.g. HYCOM, ROMS). Similar to GDAL, it does not appear that multidimensionality was designed into the current versions of ArcGIS. Like GDAL, ArcGIS appears mainly designed to work with 2D images. Both support multi-band 2D images very well, and this can be used to at least access 3D data although there is not good support for understanding the value of the z coordinate. As far as I am aware, ESRI is not currently making significant investments in the area of multidimensional data. I am part of a working group initiated by ESRI and led by an ESRI program manager (Nawajish Noman) that is trying, in essence, to get the community of users who use both ArcGIS and netCDF to develop some Python geoprocessing tools for ArcGIS that provide more functionality than out-of-box tools already in ArcGIS. ESRI did not provide any of their own developers. I donated some dev time and prototyped a tool for them using the netcdf4-python package and others are carrying on from there. I do not believe this project will result in anything dramatically new but it might yield a few additional useful geoprocessing tools. I have reduced my participation after concluding that it makes more sense to focus on the development of my own project (see below). Here are some options I'm aware of with ArcGIS: 1. The Make NetCDF Raster Layer tool can represent 3D netCDF variables as multiband raster layers. Using the Dimension Values parameter, you can slice a 3D slab out of a 4D or higher-dimensional variable. This tool might be sufficient for your purposes. The trick is your netCDF has to meet a bunch of constraints for ArcGIS to recognize it. It has to have square cells. It has adhere to the CF or COARDS conventions (I forget which versions) and have all the right attributes. You can learn more about this in the ArcGIS Desktop Help under netCDF data, FAQs about. Also see netCDF data, animating and the other topics under netCDF data. 2. Under contract to NOAA, Applied Science Associates built a couple of tools that might be useful: the Environmental Data Connector (EDC) and TimeSlider Extension. Download from http://www.asascience.com/software/downloads/index.shtml, see other parts of the website for more info. EDC was built to download multidimensional OPeNDAP datasets into multiband rasters. TimeSlider is a UI extension to help with playback of time-series data. I think they can both work with netCDFs directly, not just OPeNDAP. Also look at ASA ViewNcDap, which can display multidimensional netCDF variables. (Although this last program is not an ArcGIS extension, and there are many other programs in the world that are also not ArcGIS extensions and that probably do as good or better a job.) 3. If you don't want to use netCDFs, you can fake multidimensionality for some scenarios by building a raster catalog with columns for the time and depth. For example, you can build an animation that selects slices by date as the animation plays through a range of dates. 4. My group, the Duke University Marine Geospatial Ecology Lab, is currently building 3D and 4D awareness into a collection of tools we publish, Marine Geospatial Ecology Tools (MGET, see http://code.nicholas.duke.edu/projects/mget), built in Python on GDAL and other FOSS packages. You can download MGET today, but the 3D and 4D awareness will not be there for several months so it probably won't be useful to you in the present form. When complete, these tools will allow you to run some common marine ecology scenarios directly from ArcGIS, such as building climatologies, sampling 3D and 4D datasets at points in time and space, identifying fronts, etc. Good luck, Jason -Original Message- From: gdal-dev-boun...@lists.osgeo.org [mailto:gdal-dev-boun...@lists.osgeo.org] On Behalf Of Christopher Barker Sent: Thursday, June 17, 2010 3:16 PM To: gdal-dev@lists.osgeo.org Subject: [gdal-dev] How to represent multi-dimensional array Hi folks, I have a dataset that is a 4-dimensional array of values: time,x,y,z We're currently using netcdf to store it, which is well suited to this kind of data. However, we also need to get it into a GIS (Arc in this case), and I'm trying to find a good way to do that. Both Arc and gdalinfo do strange things if I simply point them at the netcdf file. GDAL seems to (arbitrarily?) see it as a 9064 band data set, so it is taking particular slice (I think I'm getting a bunch of (time, z) 2-d bands. Anyway, I suspect that if I re-arrange the axis in the netcdf file, I might get something more reasonable, but my question is: What is a good format to ex[press this to a GIS system? I'm imaging multiple files, maybe geo-tiff, but how to I express time and elevation in a way that is natural to GIS? thanks for your thoughts, -Chris -- Christopher Barker, Ph.D. Oceanographer Emergency Response Division NOAA/NOS/ORR
Re: [gdal-dev] Re: ESRI file geodatabase support
And ... I think I may have found the file geodatabase documentation promised ... in the online ArecGIS 10 documentation, Administrator Library, Architecture of a geodatabase, Geodatabase XML. It would appear that, in addition to the XML Schema they published several years ago, the XML documentation has been augmented with a view to interoperability being achieved at the XML level. So, maybe no, they have not published the binary structure but they may have published enough to enable XML exchange. I've not read all the material - and I would need to think about much of it before I understood it. Whether there is enough there I cannot say ... I was looking at http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html Best wishes, Peter Duarte Carreira wrote: Ragi, That is a good point, as always. Unfortunately I am not able to do that, but I can help in any other way... The zigGIS driver (I think ) is for reading/writing PostGIS using ArcMap, without having ArcSDE in the middle. As far as I know it does not read File geodatabases... I think for now you would have to reverse engineer the way fgdb's are written... given the effort, the discussion on whether to invest the time on fgbd or invest it on SL is at least interesting... and even fun for some of us. Duarte De: Ragi Burhum [mailto:r...@burhum.com] Enviada: quinta-feira, 17 de Junho de 2010 21:57 Para: gdal-dev@lists.osgeo.org Cc: Duarte Carreira; Eric Wolf; Peter J Halls; Matt Wilkie Assunto: RE: [gdal-dev] Re: ESRI file geodatabase support From: Duarte Carreira dcarre...@edia.ptmailto:dcarre...@edia.pt Subject: RE: [gdal-dev] Re: ESRI file geodatabase support To: Eric Wolf ebw...@gmail.commailto:ebw...@gmail.com, Peter J Halls p.ha...@york.ac.ukmailto:p.ha...@york.ac.uk Cc: gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org, Matt Wilkie map...@gmail.commailto:map...@gmail.com Well, if SpatiaLite offers some proper benefits and disseminates through all of the FOSS world, then it may get a strong enough push even for ESRI to pick it up. It happened before... (kml?) If SL would: 1) Be as fast as shapefile in production settings, desktop and webgis 2) Offer SQL support, spatial and otherwise, also through desktop tools like QGIS 3) Allow editing while serving (even if for 1 editor only) 4) Better support in QGIS than for shapefile (take advantage of Spatial SQL, all other functionality) 5) Same for MapServer, GeoServer, gvSIG, et al. 6) Allow easy managing of rasters inside the .db file, through QGIS 7) ??more ideas/requests?? Then it would be a very, very good contender... and the ball would be kicked to the other side. And it seems we're already there for some of the listed features. Duarte Or even better, instead of waiting and complaining, we could just write a GeoDatabase OGR (ESRI) Workspace ourselves. There are already examples of working ones out there http://svn.obtusesoft.com/core/trunk/ So instead of hoping and pleading for support, someone should just sit down and write it. I started one at one point (C++), got side tracked with other things. If anyone is interested in that source code, I would be happy to share that too. - Ragi -- Peter J Halls, GIS Advisor, University of York Telephone: 01904 433806 Fax: 01904 433740 Snail mail: Computing Service, University of York, Heslington, York YO10 5DD This message has the status of a private and personal communication ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] gdal_warp
Hi Guys I'm reprojecting some rasters from OSGB (EPSG:27700) to WGS84 (EPSG:4326) and am running into a problem. The rasters are reprojecting nicely, but I am left with black pixels around the image where the reprojected image doesnt fit squarely within the rectangular dimensions of the file. I posted about this issue a few weeks ago, and then realised that I should be marking that colour as transparent when rendering using GeoServer, and ho! everything was OK. However, I am now reprojecting some greyscale rasters, and by making the black pixels on the outside transparent, I am losing the black details of the raster image that I want to display! Is there a way of selecting the colour that gdal uses to populate these dead pixels around the side? MTIA Rob ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Re: ESRI file geodatabase support
Could well be the case. P. Jason Roberts wrote: Peter, Although I have not seen the acknowledgement you mention, I have noticed that ArcGIS 9.x appears to create SQLite databases in the user's temp directory during certain operations. Perhaps they are just using it internally to implement various processing. Jason -Original Message- From: gdal-dev-boun...@lists.osgeo.org [mailto:gdal-dev-boun...@lists.osgeo.org] On Behalf Of Peter J Halls Sent: Friday, June 18, 2010 9:09 AM To: Duarte Carreira Cc: gdal-dev@lists.osgeo.org; Matt Wilkie; Ragi Burhum Subject: Re: [gdal-dev] Re: ESRI file geodatabase support I happened to be looking at the ESRI ArcGIS 10 online documentation and found a page of acknowledgements, which includes an acknowledgement reference to SQLite. I've no idea what to make of this - maybe someone else has? Best wishes, Peter Duarte Carreira wrote: Ragi, That is a good point, as always. Unfortunately I am not able to do that, but I can help in any other way... The zigGIS driver (I think ) is for reading/writing PostGIS using ArcMap, without having ArcSDE in the middle. As far as I know it does not read File geodatabases... I think for now you would have to reverse engineer the way fgdb's are written... given the effort, the discussion on whether to invest the time on fgbd or invest it on SL is at least interesting... and even fun for some of us. Duarte De: Ragi Burhum [mailto:r...@burhum.com] Enviada: quinta-feira, 17 de Junho de 2010 21:57 Para: gdal-dev@lists.osgeo.org Cc: Duarte Carreira; Eric Wolf; Peter J Halls; Matt Wilkie Assunto: RE: [gdal-dev] Re: ESRI file geodatabase support From: Duarte Carreira dcarre...@edia.ptmailto:dcarre...@edia.pt Subject: RE: [gdal-dev] Re: ESRI file geodatabase support To: Eric Wolf ebw...@gmail.commailto:ebw...@gmail.com, Peter J Halls p.ha...@york.ac.ukmailto:p.ha...@york.ac.uk Cc: gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org, Matt Wilkie map...@gmail.commailto:map...@gmail.com Well, if SpatiaLite offers some proper benefits and disseminates through all of the FOSS world, then it may get a strong enough push even for ESRI to pick it up. It happened before... (kml?) If SL would: 1) Be as fast as shapefile in production settings, desktop and webgis 2) Offer SQL support, spatial and otherwise, also through desktop tools like QGIS 3) Allow editing while serving (even if for 1 editor only) 4) Better support in QGIS than for shapefile (take advantage of Spatial SQL, all other functionality) 5) Same for MapServer, GeoServer, gvSIG, et al. 6) Allow easy managing of rasters inside the .db file, through QGIS 7) ??more ideas/requests?? Then it would be a very, very good contender... and the ball would be kicked to the other side. And it seems we're already there for some of the listed features. Duarte Or even better, instead of waiting and complaining, we could just write a GeoDatabase OGR (ESRI) Workspace ourselves. There are already examples of working ones out there http://svn.obtusesoft.com/core/trunk/ So instead of hoping and pleading for support, someone should just sit down and write it. I started one at one point (C++), got side tracked with other things. If anyone is interested in that source code, I would be happy to share that too. - Ragi -- Peter J Halls, GIS Advisor, University of York Telephone: 01904 433806 Fax: 01904 433740 Snail mail: Computing Service, University of York, Heslington, York YO10 5DD This message has the status of a private and personal communication ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] gdal_warp
Rob, You can try the -dstalpha option in gdalwarp to create an alpha band. ( http://www.gdal.org/gdalwarp.html) On Fri, Jun 18, 2010 at 7:22 PM, Rob g...@vanbooth.com wrote: Hi Guys I'm reprojecting some rasters from OSGB (EPSG:27700) to WGS84 (EPSG:4326) and am running into a problem. The rasters are reprojecting nicely, but I am left with black pixels around the image where the reprojected image doesnt fit squarely within the rectangular dimensions of the file. I posted about this issue a few weeks ago, and then realised that I should be marking that colour as transparent when rendering using GeoServer, and ho! everything was OK. However, I am now reprojecting some greyscale rasters, and by making the black pixels on the outside transparent, I am losing the black details of the raster image that I want to display! Is there a way of selecting the colour that gdal uses to populate these dead pixels around the side? MTIA Rob ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Best regards, Chaitanya kumar CH. /tʃaɪθənjə/ /kʊmɑr/ +91-9494447584 17.2416N 80.1426E ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
RE: [gdal-dev] Re: ESRI file geodatabase support
Just a quick note: SQLite is supported through the FME add-on aka Interoperability Extension to ArcMap. Duarte -Mensagem original- De: Peter J Halls [mailto:p.ha...@york.ac.uk] Enviada: sexta-feira, 18 de Junho de 2010 14:49 Para: Duarte Carreira Cc: Ragi Burhum; gdal-dev@lists.osgeo.org; Eric Wolf; Matt Wilkie Assunto: Re: [gdal-dev] Re: ESRI file geodatabase support And ... I think I may have found the file geodatabase documentation promised ... in the online ArecGIS 10 documentation, Administrator Library, Architecture of a geodatabase, Geodatabase XML. It would appear that, in addition to the XML Schema they published several years ago, the XML documentation has been augmented with a view to interoperability being achieved at the XML level. So, maybe no, they have not published the binary structure but they may have published enough to enable XML exchange. I've not read all the material - and I would need to think about much of it before I understood it. Whether there is enough there I cannot say ... I was looking at http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html Best wishes, Peter Duarte Carreira wrote: Ragi, That is a good point, as always. Unfortunately I am not able to do that, but I can help in any other way... The zigGIS driver (I think ) is for reading/writing PostGIS using ArcMap, without having ArcSDE in the middle. As far as I know it does not read File geodatabases... I think for now you would have to reverse engineer the way fgdb's are written... given the effort, the discussion on whether to invest the time on fgbd or invest it on SL is at least interesting... and even fun for some of us. Duarte De: Ragi Burhum [mailto:r...@burhum.com] Enviada: quinta-feira, 17 de Junho de 2010 21:57 Para: gdal-dev@lists.osgeo.org Cc: Duarte Carreira; Eric Wolf; Peter J Halls; Matt Wilkie Assunto: RE: [gdal-dev] Re: ESRI file geodatabase support From: Duarte Carreira dcarre...@edia.ptmailto:dcarre...@edia.pt Subject: RE: [gdal-dev] Re: ESRI file geodatabase support To: Eric Wolf ebw...@gmail.commailto:ebw...@gmail.com, Peter J Halls p.ha...@york.ac.ukmailto:p.ha...@york.ac.uk Cc: gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org gdal-dev@lists.osgeo.orgmailto:gdal-dev@lists.osgeo.org, Matt Wilkie map...@gmail.commailto:map...@gmail.com Well, if SpatiaLite offers some proper benefits and disseminates through all of the FOSS world, then it may get a strong enough push even for ESRI to pick it up. It happened before... (kml?) If SL would: 1) Be as fast as shapefile in production settings, desktop and webgis 2) Offer SQL support, spatial and otherwise, also through desktop tools like QGIS 3) Allow editing while serving (even if for 1 editor only) 4) Better support in QGIS than for shapefile (take advantage of Spatial SQL, all other functionality) 5) Same for MapServer, GeoServer, gvSIG, et al. 6) Allow easy managing of rasters inside the .db file, through QGIS 7) ??more ideas/requests?? Then it would be a very, very good contender... and the ball would be kicked to the other side. And it seems we're already there for some of the listed features. Duarte Or even better, instead of waiting and complaining, we could just write a GeoDatabase OGR (ESRI) Workspace ourselves. There are already examples of working ones out there http://svn.obtusesoft.com/core/trunk/ So instead of hoping and pleading for support, someone should just sit down and write it. I started one at one point (C++), got side tracked with other things. If anyone is interested in that source code, I would be happy to share that too. - Ragi -- Peter J Halls, GIS Advisor, University of York Telephone: 01904 433806 Fax: 01904 433740 Snail mail: Computing Service, University of York, Heslington, York YO10 5DD This message has the status of a private and personal communication ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] gdal_warp
Thanks for the replies so far, but I'm obviously missing a trick as it doesnt seem to be working. The gdalinfo on my TIF looks like this. E:\rastersgdalinfo tq00ne.tif Driver: GTiff/GeoTIFF Files: tq00ne.tif tq00ne.tfw Size is 7874, 7874 Coordinate System is `' Origin = (505000.000,11.000) Pixel Size = (0.635001270002540,-0.635001270002540) Metadata: TIFFTAG_IMAGEDESCRIPTION=1:1 TILE TQ00NE TIFFTAG_DATETIME=2008:07:14 23:22:37 TIFFTAG_COPYRIGHT=ORDNANCE SURVEY CROWN COPYRIGHT 2008 TIFFTAG_XRESOLUTION=400 TIFFTAG_YRESOLUTION=400 TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) Image Structure Metadata: COMPRESSION=LZW INTERLEAVE=BAND Corner Coordinates: Upper Left ( 505000.000, 11.000) Lower Left ( 505000.000, 105000.000) Upper Right ( 51.000, 11.000) Lower Right ( 51.000, 105000.000) Center ( 507500.000, 107500.000) Band 1 Block=7874x8 Type=Byte, ColorInterp=Palette Image Structure Metadata: NBITS=4 Color Table (RGB with 16 entries) 0: 0,0,0,255 1: 107,107,107,255 2: 90,90,90,255 3: 220,220,220,255 4: 209,209,209,255 5: 243,243,243,255 6: 190,190,190,255 7: 243,243,243,255 8: 190,190,190,255 9: 255,255,255,255 10: 0,0,0,255 11: 0,0,0,255 12: 0,0,0,255 13: 0,0,0,255 14: 0,0,0,255 15: 0,0,0,255 But whenever I try to do a gdalwarp with the using -dstalpha I get an error message and an unsuable TIF file as output E:\rastersgdalwarp -of GTiff -s_srs EPSG:27700 -t_srs EPSG:4326 -co TFW=YES -co TILED=YES -co COMPRESS=LZW -dstalpha TQ00NE.TIF TQwarpa.TIF Copying color table from TQ00NE.TIF to new file. Creating output file that is 9687P x 6132L. ERROR 6: SetColorTable() not supported for multi-sample TIFF files. Processing input file TQ00NE.TIF. 0...10...20...30...40...50...60...70...80...90...100 - done. E:\rastersgdalwarp -of GTiff -s_srs EPSG:27700 -t_srs EPSG:4326 -co TFW=YES -co TILED=YES -co COMPRESS=LZW -dstnodata 12 -dstalpha TQ00NE.TIF TQwarpa.TIF Copying color table from TQ00NE.TIF to new file. Creating output file that is 9687P x 6132L. ERROR 6: SetColorTable() not supported for multi-sample TIFF files. Processing input file TQ00NE.TIF. 0...10...20...30...40...50...60...70...80...90...100 - done. Any more suggestions? Thanks Rob On 18 June 2010 15:06, Chaitanya kumar CH chaitanya...@gmail.com wrote: Rob, You can try the -dstalpha option in gdalwarp to create an alpha band. ( http://www.gdal.org/gdalwarp.html) On Fri, Jun 18, 2010 at 7:22 PM, Rob g...@vanbooth.com wrote: Hi Guys I'm reprojecting some rasters from OSGB (EPSG:27700) to WGS84 (EPSG:4326) and am running into a problem. The rasters are reprojecting nicely, but I am left with black pixels around the image where the reprojected image doesnt fit squarely within the rectangular dimensions of the file. I posted about this issue a few weeks ago, and then realised that I should be marking that colour as transparent when rendering using GeoServer, and ho! everything was OK. However, I am now reprojecting some greyscale rasters, and by making the black pixels on the outside transparent, I am losing the black details of the raster image that I want to display! Is there a way of selecting the colour that gdal uses to populate these dead pixels around the side? MTIA Rob ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Best regards, Chaitanya kumar CH. /tʃaɪθənjə/ /kʊmɑr/ +91-9494447584 17.2416N 80.1426E ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Re: ESRI file geodatabase support
This is exactly what I have been talking about. On Fri, Jun 18, 2010 at 5:37 AM, Duarte Carreira dcarre...@edia.pt wrote: Re-reading your email it seems I misunderstood what you meant… So using the zigGIS ArcMap/PostGIS provider one could adapt it to read OGR datasources, like QGIS does. You would then access SpatiaLite and potentially all other OGR formats from ArcMap. Good idea. Duarte ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Re: ESRI file geodatabase support
Right. And you have to fork another $3000 (or $1500?) **per seat** for one of these licenses **on top** of what you already paid for ArcGIS desktop. On Fri, Jun 18, 2010 at 6:16 AM, Duarte Carreira dcarre...@edia.pt wrote: Just a quick note: SQLite is supported through the FME add-on aka Interoperability Extension to ArcMap. Duarte ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Re: ESRI file geodatabase support
This is meant to be an export/import type of thing. So you won't be able to open something in ArcMap and start manipulating it live using the ArcMap/ArcCatalog tools. The other approach I am proposing would be OGR as a seamless datasource that integrates with ArcGIS at the data access layer level. It would be transparent to the other tools. - Ragi On Fri, Jun 18, 2010 at 5:49 AM, Peter J Halls p.ha...@york.ac.uk wrote: And ... I think I may have found the file geodatabase documentation promised ... in the online ArecGIS 10 documentation, Administrator Library, Architecture of a geodatabase, Geodatabase XML. It would appear that, in addition to the XML Schema they published several years ago, the XML documentation has been augmented with a view to interoperability being achieved at the XML level. So, maybe no, they have not published the binary structure but they may have published enough to enable XML exchange. I've not read all the material - and I would need to think about much of it before I understood it. Whether there is enough there I cannot say ... I was looking at http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html Best wishes, Peter ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
RE: [gdal-dev] Re: ESRI file geodatabase support
Well, I'm sold. What you need to do it?? Duarte De: Ragi Burhum [mailto:r...@burhum.com] Enviada: sexta-feira, 18 de Junho de 2010 16:31 Para: Peter J Halls Cc: Duarte Carreira; gdal-dev@lists.osgeo.org; Eric Wolf; Matt Wilkie Assunto: Re: [gdal-dev] Re: ESRI file geodatabase support This is meant to be an export/import type of thing. So you won't be able to open something in ArcMap and start manipulating it live using the ArcMap/ArcCatalog tools. The other approach I am proposing would be OGR as a seamless datasource that integrates with ArcGIS at the data access layer level. It would be transparent to the other tools. - Ragi On Fri, Jun 18, 2010 at 5:49 AM, Peter J Halls p.ha...@york.ac.ukmailto:p.ha...@york.ac.uk wrote: And ... I think I may have found the file geodatabase documentation promised ... in the online ArecGIS 10 documentation, Administrator Library, Architecture of a geodatabase, Geodatabase XML. It would appear that, in addition to the XML Schema they published several years ago, the XML documentation has been augmented with a view to interoperability being achieved at the XML level. So, maybe no, they have not published the binary structure but they may have published enough to enable XML exchange. I've not read all the material - and I would need to think about much of it before I understood it. Whether there is enough there I cannot say ... I was looking at http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html Best wishes, Peter ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Problem with pixel/line coordinates
Thanks for your explanation, Frank! Now it works fine On Thu, 17 Jun 2010 13:59:46 -0400 Frank Warmerdam warmer...@pobox.com wrote: Alexander, I believe the addition of 0.5 in the above is incorrect. In the simple, non-rotated case, all values from geoTransform[0] to geoTransform[0] + geoTransform[1] would be on the 1st pixel (ie. pX = 0). As you have coded it, when you are half way across the pixel you are switching into the next one. Keep in mind that (geoTransform[0],geoTransform[3]) is the top left corner of the top left pixel - not the center. def pixelToMap( pX, pY, geoTransform ): mX, mY = applyGeoTransform( pX, pY, geoTransform ) return mX, mY Conversely, here if pX and pY are coming in as integer rather than floating point values, then you will likely want to add half a pixel before transforming so that you get the geoeferenced location of the center of the pixel rather than the upper left corner. -- Alexander Bruy ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] gdal on OS X
Hi List, I'm trying to get gdal working on OS X 10.6. I'm have tried a variety of ways with varying success. I installed GDAL 1.7 Complete from http://www.kyngchaos.com/software:frameworks If I go to /Library/Frameworks/GDAL.framework/unix/bin I can seem to run the unix binaries. For example, ./gdal_info provides usage instructions. I cannot run the Python programs. If I launch my Python and type from osgeo import gdal it doesn't work because I run the EPD Python not the OS X default. That command *does* work, and the Python binaries do run (or at least print usage info) if I comment out the EPD section in my .bash_profile and run using the system Python. Any ideas how to run GDAL with EPD? My second attempt was to build gdal from source. export MACOSX_DEPLOYMENT_TARGET=10.5 export ARCHFLAGS=-arch i386 ./configure --prefix=/Users/mankoff/local/gdal-1.7.2 --with-macosx-framework --with-python CFLAGS=-m32 CPPFLAGS=-m32 CXXFLAGS=-m32 make The make command crashes with /bin/sh /Users/mankoff/local/src/gdal-1.7.2/libtool --mode=link g++ gdalinfo.lo /Users/mankoff/local/src/gdal-1.7.2/libgdal.la -o gdalinfo libtool: link: g++ .libs/gdalinfo.o -o .libs/gdalinfo /Users/mankoff/local/src/gdal-1.7.2/.libs/libgdal.dylib -L/usr/lib -lsqlite3 -lexpat -lz -ldl -lcurl ld: warning: in .libs/gdalinfo.o, file was built for i386 which is not the architecture being linked (x86_64) ld: warning: in /Users/mankoff/local/src/gdal-1.7.2/.libs/libgdal.dylib, file was built for i386 which is not the architecture being linked (x86_64) Undefined symbols: _main, referenced from: start in crt1.10.5.o ld: symbol(s) not found collect2: ld returned 1 exit status make[1]: *** [gdalinfo] Error 1 make: *** [apps-target] Error 2 Clearly a 32/64 bit issue. I think I'm missing some simple link flag to tell it what architecture I'm aiming for. Does anyone know what it is? I'd love to type fink install gdal but for some reason it depends on and also installs qt, python, etc. There seem to be an awful lot of dependencies that the raw source compilation doesn't use. Thanks for any help you can provide, -k. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] gdal on OS X
On Jun 18, 2010, at 2:42 PM, Ken Mankoff wrote: Hi List, I'm trying to get gdal working on OS X 10.6. I'm have tried a variety of ways with varying success. I installed GDAL 1.7 Complete from http://www.kyngchaos.com/software:frameworks If I go to /Library/Frameworks/GDAL.framework/unix/bin I can seem to run the unix binaries. For example, ./gdal_info provides usage instructions. I cannot run the Python programs. If I launch my Python and type from osgeo import gdal it doesn't work because I run the EPD Python not the OS X default. That command *does* work, and the Python binaries do run (or at least print usage info) if I comment out the EPD section in my .bash_profile and run using the system Python. Any ideas how to run GDAL with EPD? If the EPD Python you have installed is the same major version as the system (2.6), this should work: try copying the gdal path file from /Library/Python/2.6/site-packages/gdal.pth to the user site-packages folder of the EPD python. My GDAL python binaries are built for Python 2.6, and may (probably?) not work on Python 2.5 or 2.7. Definitely (?) not Python 3. - William Kyngesburye kyngchaos*at*kyngchaos*dot*com http://www.kyngchaos.com/ [Trillian] What are you supposed to do WITH a maniacally depressed robot? [Marvin] You think you have problems? What are you supposed to do if you ARE a maniacally depressed robot? No, don't try and answer, I'm 50,000 times more intelligent than you and even I don't know the answer... - HitchHiker's Guide to the Galaxy ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] gdal on OS X
On Fri, 18 Jun 2010, William Kyngesburye wrote: On Jun 18, 2010, at 2:42 PM, Ken Mankoff wrote: Hi List, I'm trying to get gdal working on OS X 10.6. I'm have tried a variety of ways with varying success. I installed GDAL 1.7 Complete from http://www.kyngchaos.com/software:frameworks If I go to /Library/Frameworks/GDAL.framework/unix/bin I can seem to run the unix binaries. For example, ./gdal_info provides usage instructions. I cannot run the Python programs. If I launch my Python and type from osgeo import gdal it doesn't work because I run the EPD Python not the OS X default. That command *does* work, and the Python binaries do run (or at least print usage info) if I comment out the EPD section in my .bash_profile and run using the system Python. Any ideas how to run GDAL with EPD? If the EPD Python you have installed is the same major version as the system (2.6), this should work: try copying the gdal path file from /Library/Python/2.6/site-packages/gdal.pth to the user site-packages folder of the EPD python. My GDAL python binaries are built for Python 2.6, and may (probably?) not work on Python 2.5 or 2.7. Definitely (?) not Python 3. Yup, this works for me: cp /Library/Python/2.6/site-packages/gdal.pth /Library/Frameworks/Python.framework/Versions/Current/lib/python2.6/site-packages/ Thanks, -k. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev