source: trunk/zoo-project/zoo-kernel/service_internal_ms.c @ 557

Last change on this file since 557 was 550, checked in by djay, 10 years ago

Add otb2zcfg and OTB applications support without observer by now. Fix issue with maxOccurs and multiple downloaded value for the same input.
  • Property svn:eol-style set to native
  • Property svn:mime-type set to text/x-csrc
File size: 34.4 KB
Line 
1/**
2 * Author : Gérald FENOY
3 *
4 *  Copyright 2010-2011 Fondazione Edmund Mach. All rights reserved.
5 *
6 * Permission is hereby granted, free of charge, to any person obtaining a copy
7 * of this software and associated documentation files (the "Software"), to deal
8 * in the Software without restriction, including without limitation the rights
9 * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
10 * copies of the Software, and to permit persons to whom the Software is
11 * furnished to do so, subject to the following conditions:
12 *
13 * The above copyright notice and this permission notice shall be included in
14 * all copies or substantial portions of the Software.
15 *
16 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
17 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
18 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
19 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
20 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
21 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
22 * THE SOFTWARE.
23 */
24
25#ifdef USE_MS
26#ifndef WIN32
27#define CLASS class
28#else
29#define CLASS _class
30#endif
31#include "service_internal_ms.h"
32
33/**
34 * Map composed by a main.cfg maps name as key and the corresponding
35 * MapServer Mafile Metadata name to use
36 * see doc from here :
37 *  - http://mapserver.org/ogc/wms_server.html
38 *  - http://mapserver.org/ogc/wfs_server.html
39 *  - http://mapserver.org/ogc/wcs_server.html
40 */
41map* getCorrespondance(){
42  map* res=createMap("encoding","ows_encoding");
43  addToMap(res,"abstract","ows_abstract");
44  addToMap(res,"title","ows_title");
45  addToMap(res,"keywords","ows_keywordlist");
46  addToMap(res,"fees","ows_fees");
47  addToMap(res,"accessConstraints","ows_accessconstraints");
48  addToMap(res,"providerName","ows_attribution_title");
49  addToMap(res,"providerSite","ows_service_onlineresource");
50  addToMap(res,"individualName","ows_contactperson");
51  addToMap(res,"positionName","ows_contactposition");
52  addToMap(res,"providerName","ows_contactorganization");
53  addToMap(res,"role","ows_role");
54  addToMap(res,"addressType","ows_addresstype");
55  addToMap(res,"addressCity","ows_city");
56  addToMap(res,"addressDeliveryPoint","ows_address");
57  addToMap(res,"addressPostalCode","ows_postcode");
58  addToMap(res,"addressAdministrativeArea","ows_stateorprovince");
59  addToMap(res,"addressCountry","ows_country");
60  addToMap(res,"phoneVoice","ows_contactvoicetelephone");
61  addToMap(res,"phoneFacsimile","ows_contactfacsimiletelephone");
62  addToMap(res,"addressElectronicMailAddress","ows_contactelectronicmailaddress");
63  // Missing Madatory Informations
64  addToMap(res,"hoursOfService","ows_hoursofservice");
65  addToMap(res,"contactInstructions","ows_contactinstructions");
66  return res;
67}
68
69void setMapSize(maps* output,double minx,double miny,double maxx,double maxy){
70  double maxWidth=640;
71  double maxHeight=480;
72  double deltaX=maxx-minx;
73  double deltaY=maxy-miny;
74  double qWidth;
75  qWidth=maxWidth/deltaX;
76  double qHeight;
77  qHeight=maxHeight/deltaY;
78#ifdef DEBUGMS
79  fprintf(stderr,"deltaX : %.15f \ndeltaY : %.15f\n",deltaX,deltaY);
80  fprintf(stderr,"qWidth : %.15f \nqHeight : %.15f\n",qWidth,qHeight);
81#endif
82
83  double width=deltaX*qWidth;
84  double height=height=deltaY*qWidth;
85  if(deltaX<deltaY){
86    width=deltaX*qHeight;
87    height=deltaY*qHeight;
88  }
89  if(height<0)
90    height=-height;
91  if(width<0)
92    width=-width;
93  char sWidth[1024];
94  char sHeight[1024];
95  sprintf(sWidth,"%.3f",width);
96  sprintf(sHeight,"%.3f",height);
97#ifdef DEBUGMS
98  fprintf(stderr,"sWidth : %.15f \nsHeight : %.15f\n",sWidth,sHeight);
99#endif
100  if(output!=NULL){
101    addToMap(output->content,"width",sWidth);
102    addToMap(output->content,"height",sHeight);
103  }
104}
105
106void setReferenceUrl(maps* m,maps* tmpI){
107  outputMapfile(m,tmpI);
108  map *msUrl=getMapFromMaps(m,"main","mapserverAddress");
109  map *msOgcVersion=getMapFromMaps(m,"main","msOgcVersion");
110  map *dataPath=getMapFromMaps(m,"main","dataPath");
111  map *sid=getMapFromMaps(m,"lenv","usid");
112  map* format=getMap(tmpI->content,"mimeType");
113  map* rformat=getMap(tmpI->content,"requestedMimeType");
114  map* width=getMap(tmpI->content,"width");
115  map* height=getMap(tmpI->content,"height");
116  map* protoMap=getMap(tmpI->content,"msOgc");
117  map* versionMap=getMap(tmpI->content,"msOgcVersion");
118  char options[3][5][25]={
119    {"WMS","1.3.0","GetMap","layers=%s","wms_extent"},
120    {"WFS","1.1.0","GetFeature","typename=%s","wcs_extent"},
121    {"WCS","1.1.0","GetCoverage","coverage=%s","wcs_extent"}
122  };
123  int proto=0;
124  if(rformat==NULL){
125    rformat=getMap(tmpI->content,"mimeType");
126  }
127  if(strncasecmp(rformat->value,"text/xml",8)==0)
128    proto=1;
129  if(strncasecmp(rformat->value,"image/tiff",10)==0)
130    proto=2;
131  if(protoMap!=NULL){
132    if(strncasecmp(protoMap->value,"WMS",3)==0)
133      proto=0;
134    else{
135      if(strncasecmp(protoMap->value,"WFS",3)==0)
136        proto=1;
137      else 
138        proto=2;
139    }
140  }
141  char *protoVersion=options[proto][1];
142  if(proto==1){
143    if(msOgcVersion!=NULL)
144      protoVersion=msOgcVersion->value;
145    if(versionMap!=NULL)
146      protoVersion=versionMap->value;
147  }
148
149  map* extent=getMap(tmpI->content,options[proto][4]);
150  map* crs=getMap(tmpI->content,"crs");
151  int hasCRS=1;
152  if(crs==NULL){
153    crs=getMapFromMaps(m,"main","crs");
154    if(crs==NULL){
155      crs=createMap("crs","epsg:4326");
156      hasCRS=0;
157    }
158  }
159  char layers[128];
160  sprintf(layers,options[proto][3],tmpI->name);
161
162  char* webService_url=(char*)malloc((strlen(msUrl->value)+strlen(format->value)+strlen(tmpI->name)+strlen(width->value)+strlen(height->value)+strlen(extent->value)+256)*sizeof(char));
163
164  if(proto>0){
165    sprintf(webService_url,
166            "%s?map=%s/%s_%s.map&request=%s&service=%s&version=%s&%s&format=%s&bbox=%s&crs=%s",
167            msUrl->value,
168            dataPath->value,
169            tmpI->name,
170            sid->value,
171            options[proto][2],
172            options[proto][0],
173            protoVersion,
174            layers,
175            rformat->value,
176            extent->value,
177            crs->value
178            );
179  }
180  else{
181    sprintf(webService_url,
182            "%s?map=%s/%s_%s.map&request=%s&service=%s&version=%s&%s&width=%s&height=%s&format=%s&bbox=%s&crs=%s",
183            msUrl->value,
184            dataPath->value,
185            tmpI->name,
186            sid->value,
187            options[proto][2],
188            options[proto][0],
189            protoVersion,
190            layers,
191            width->value,
192            height->value,
193            rformat->value,
194            extent->value,
195            crs->value
196            );
197  }
198  if(hasCRS==0){
199    freeMap(&crs);
200    free(crs);
201  }
202  addToMap(tmpI->content,"Reference",webService_url);
203  free(webService_url);
204}
205
206/**
207 * Set projection using Authority Code and Name if available or fallback to
208 * proj4 definition if available or fallback to default EPSG:4326
209 */
210void setSrsInformations(maps* output,mapObj* m,layerObj* myLayer,
211                        char* pszProjection){
212  OGRSpatialReferenceH  hSRS;
213  map* msSrs=NULL;
214  hSRS = OSRNewSpatialReference(NULL);
215  if( pszProjection!=NULL && strlen(pszProjection)>1){
216    if(OSRImportFromWkt( hSRS, &pszProjection ) == CE_None ){
217      char *proj4Str=NULL;
218      if(OSRGetAuthorityName(hSRS,NULL)!=NULL && 
219         OSRGetAuthorityCode(hSRS,NULL)!=NULL){
220        char tmpSrs[20];
221        sprintf(tmpSrs,"%s:%s",
222                OSRGetAuthorityName(hSRS,NULL),OSRGetAuthorityCode(hSRS,NULL));
223        msLoadProjectionStringEPSG(&m->projection,tmpSrs);
224        msLoadProjectionStringEPSG(&myLayer->projection,tmpSrs);
225       
226        char tmpSrss[256];
227        sprintf(tmpSrss,"EPSG:4326 EPSG:900913 %s",tmpSrs);
228       
229        msInsertHashTable(&(m->web.metadata), "ows_srs", tmpSrss);
230        msInsertHashTable(&(myLayer->metadata), "ows_srs", tmpSrss);
231       
232#ifdef DEBUGMS
233        fprintf(stderr,"isGeo %b\n\n",OSRIsGeographic(hSRS)==TRUE);
234#endif
235        if(output!=NULL){
236          if(OSRIsGeographic(hSRS)==TRUE)
237            addToMap(output->content,"crs_isGeographic","true");
238          else
239            addToMap(output->content,"crs_isGeographic","false");
240          addToMap(output->content,"crs",tmpSrs);
241        }
242      }
243      else{
244        OSRExportToProj4(hSRS,&proj4Str);
245        if(proj4Str!=NULL && strlen(proj4Str)>0){
246#ifdef DEBUGMS
247          fprintf(stderr,"PROJ (%s)\n",proj4Str);
248#endif
249          msLoadProjectionString(&(m->projection),proj4Str);
250          msLoadProjectionString(&(myLayer->projection),proj4Str);
251          if(output!=NULL){ 
252            if(OSRIsGeographic(hSRS)==TRUE)
253              addToMap(output->content,"crs_isGeographic","true");
254            else
255              addToMap(output->content,"crs_isGeographic","false");
256          }
257          free(proj4Str);
258        }
259        else{
260          msLoadProjectionStringEPSG(&m->projection,"EPSG:4326");
261          msLoadProjectionStringEPSG(&myLayer->projection,"EPSG:4326");
262          if(output!=NULL){
263            addToMap(output->content,"crs_isGeographic","true");
264          }
265        }
266        if(output!=NULL){
267          addToMap(output->content,"crs","EPSG:4326");
268          addToMap(output->content,"real_extent","true");
269        }
270        msInsertHashTable(&(m->web.metadata),"ows_srs", "EPSG:4326 EPSG:900913");
271        msInsertHashTable(&(myLayer->metadata),"ows_srs","EPSG:4326 EPSG:900913");
272      }
273    }
274  }
275  else{
276    if(output!=NULL){
277      msSrs=getMap(output->content,"msSrs");
278    }
279    if(msSrs!=NULL){
280      msLoadProjectionStringEPSG(&m->projection,msSrs->value);
281      msLoadProjectionStringEPSG(&myLayer->projection,msSrs->value);
282      char tmpSrs[128];
283      sprintf(tmpSrs,"%s EPSG:4326 EPSG:900913",msSrs->value);
284      msInsertHashTable(&(m->web.metadata),"ows_srs",tmpSrs);
285      msInsertHashTable(&(myLayer->metadata),"ows_srs",tmpSrs);
286    }else{
287      msLoadProjectionStringEPSG(&m->projection,"EPSG:4326");
288      msLoadProjectionStringEPSG(&myLayer->projection,"EPSG:4326");
289      msInsertHashTable(&(m->web.metadata),"ows_srs","EPSG:4326 EPSG:900913");
290      msInsertHashTable(&(myLayer->metadata),"ows_srs","EPSG:4326 EPSG:900913");
291    }
292    if(output!=NULL){
293      addToMap(output->content,"crs",msSrs->value);
294      addToMap(output->content,"crs_isGeographic","true");
295    }
296  }
297
298  OSRDestroySpatialReference( hSRS );
299}
300
301void setMsExtent(maps* output,mapObj* m,layerObj* myLayer,
302                 double minX,double minY,double maxX,double maxY){
303  msMapSetExtent(m,minX,minY,maxX,maxY);
304#ifdef DEBUGMS
305  fprintf(stderr,"Extent %.15f %.15f %.15f %.15f\n",minX,minY,maxX,maxY);
306#endif
307  char tmpExtent[1024];
308  sprintf(tmpExtent,"%.15f %.15f %.15f %.15f",minX,minY,maxX,maxY);
309#ifdef DEBUGMS
310  fprintf(stderr,"Extent %s\n",tmpExtent);
311#endif
312  msInsertHashTable(&(myLayer->metadata), "ows_extent", tmpExtent);
313 
314  if(output!=NULL){
315    map* test=getMap(output->content,"real_extent");
316    if(test!=NULL){
317      pointObj min, max;
318      projectionObj tempSrs;
319      min.x = m->extent.minx;
320      min.y = m->extent.miny;
321      max.x = m->extent.maxx;
322      max.y = m->extent.maxy;
323      char tmpSrsStr[1024];
324      msInitProjection(&tempSrs);
325      msLoadProjectionStringEPSG(&tempSrs,"EPSG:4326");
326
327      msProjectPoint(&(m->projection),&tempSrs,&min);
328      msProjectPoint(&m->projection,&tempSrs,&max);
329     
330      sprintf(tmpExtent,"%.3f,%.3f,%.3f,%.3f",min.y,min.x,max.y,max.x);
331      map* isGeo=getMap(output->content,"crs_isGeographic");
332#ifdef DEBUGMS
333      fprintf(stderr,"isGeo = %s\n",isGeo->value);
334#endif
335      if(isGeo!=NULL && strcasecmp("true",isGeo->value)==0)
336        sprintf(tmpExtent,"%f,%f,%f,%f", minY,minX, maxY, maxX);
337      addToMap(output->content,"wms_extent",tmpExtent);
338      sprintf(tmpSrsStr,"%.3f,%.3f,%.3f,%.3f",min.x,min.y,max.x,max.y);
339      addToMap(output->content,"wcs_extent",tmpExtent);
340    }else{
341      sprintf(tmpExtent,"%f,%f,%f,%f",minX, minY, maxX, maxY);
342      map* isGeo=getMap(output->content,"crs_isGeographic");
343      if(isGeo!=NULL){
344#ifdef DEBUGMS
345        fprintf(stderr,"isGeo = %s\n",isGeo->value);
346#endif
347        if(isGeo!=NULL && strcasecmp("true",isGeo->value)==0)
348          sprintf(tmpExtent,"%f,%f,%f,%f", minY,minX, maxY, maxX);
349      }
350      addToMap(output->content,"wms_extent",tmpExtent); 
351      sprintf(tmpExtent,"%.3f,%.3f,%.3f,%.3f",minX,minY,maxX,maxY);
352      addToMap(output->content,"wcs_extent",tmpExtent);
353    }
354  }
355
356  setMapSize(output,minX,minY,maxX,maxY);
357}
358
359int tryOgr(maps* conf,maps* output,mapObj* m){
360
361  map* tmpMap=getMap(output->content,"storage");
362  char *pszDataSource=tmpMap->value;
363
364  /**
365   * Try to open the DataSource using OGR
366   */
367  OGRRegisterAll();
368  /**
369   * Try to load the file as ZIP
370   */
371
372  OGRDataSourceH poDS1 = NULL;
373  OGRSFDriverH *poDriver1 = NULL;
374  char *dsName=(char*)malloc((8+strlen(pszDataSource)+1)*sizeof(char));
375  char *odsName=zStrdup(pszDataSource);
376  char *sdsName=zStrdup(pszDataSource);
377  char *demo=strstr(odsName,".");
378  sdsName[strlen(sdsName)-(strlen(demo)-1)]='d';
379  sdsName[strlen(sdsName)-(strlen(demo)-2)]='i';
380  sdsName[strlen(sdsName)-(strlen(demo)-3)]='r';
381  sdsName[strlen(sdsName)-(strlen(demo)-4)]=0;
382
383  odsName[strlen(odsName)-(strlen(demo)-1)]='z';
384  odsName[strlen(odsName)-(strlen(demo)-2)]='i';
385  odsName[strlen(odsName)-(strlen(demo)-3)]='p';
386  odsName[strlen(odsName)-(strlen(demo)-4)]=0;
387  sprintf(dsName,"/vsizip/%s",odsName);
388
389#ifdef DEBUGMS
390  fprintf(stderr,"Try loading %s, %s, %s\n",dsName,odsName,dsName);
391#endif
392
393  FILE* file = fopen(pszDataSource, "rb");
394  FILE* fileZ = fopen(odsName, "wb");
395  free(odsName);
396  fseek(file, 0, SEEK_END);
397  unsigned long fileLen=ftell(file);
398  fseek(file, 0, SEEK_SET);
399  char *buffer=(char *)malloc(fileLen+1);
400  fread(buffer, fileLen, 1, file);
401  fwrite(buffer,fileLen, 1, fileZ);
402  fclose(file);
403  fclose(fileZ);
404  free(buffer);
405#ifdef DEBUGMS
406  fprintf(stderr,"Try loading %s",dsName);
407#endif
408  poDS1 = OGROpen( dsName, FALSE, poDriver1 );
409  if( poDS1 == NULL ){
410    fprintf(stderr,"Unable to access the DataSource as ZIP File\n");
411    setMapInMaps(conf,"lenv","message","Unable to open datasource in read only mode");
412    OGR_DS_Destroy(poDS1);
413  }else{
414#ifdef DEBUGMS
415    fprintf(stderr,"The DataSource is a  ZIP File\n");
416#endif
417    char** demo=VSIReadDir(dsName);
418    int i=0;
419    zMkdir(sdsName
420#ifndef WIN32
421                ,S_IRWXU | S_IRGRP | S_IXGRP | S_IROTH | S_IXOTH
422#endif
423                );
424    while(demo[i]!=NULL){
425#ifdef DEBUGMS
426      fprintf(stderr,"ZIP File content : %s\n",demo[i]);
427#endif
428      char *tmpDs=(char*)malloc((strlen(dsName)+strlen(demo[i])+2)*sizeof(char));
429      sprintf(tmpDs,"%s/%s",dsName,demo[i]);
430      fprintf(stderr,"read : %s\n",tmpDs);
431     
432      VSILFILE* vsif=VSIFOpenL(tmpDs,"rb");
433#ifdef DEBUGMS
434      fprintf(stderr,"open : %s\n",tmpDs);
435#endif
436      VSIFSeekL(vsif,0,SEEK_END);
437      vsi_l_offset size=VSIFTellL(vsif);
438#ifdef DEBUGMS
439      fprintf(stderr,"size : %d\n",size);
440#endif
441      VSIFSeekL(vsif,0,SEEK_SET);
442      char *vsifcontent=(char*) malloc(((int)size+1)*sizeof(char));
443      VSIFReadL(vsifcontent,1,(size_t)size,vsif);
444      char *fpath=(char*) malloc((strlen(sdsName)+strlen(demo[1])+2)*sizeof(char));
445      sprintf(fpath,"%s/%s",sdsName,demo[i]);
446      int f=zOpen(fpath,O_WRONLY|O_CREAT,S_IRUSR|S_IWUSR|S_IRGRP|S_IWGRP|S_IROTH|S_IWOTH);
447      zWrite(f,vsifcontent,(int)size);
448      close(f);
449      chmod(fpath,S_IRWXU | S_IRGRP | S_IXGRP | S_IROTH);
450      char* tmpP=strstr(fpath,".shp");
451      if(tmpP==NULL)
452        tmpP=strstr(fpath,".SHP");
453      if(tmpP!=NULL){
454#ifdef DEBUGMS
455        fprintf(stderr,"*** DEBUG %s\n",strstr(tmpP,"."));
456#endif
457        if( strcmp(tmpP,".shp")==0 || strcmp(tmpP,".SHP")==0 ){
458          tmpMap=getMap(output->content,"storage");
459          free(tmpMap->value);
460          tmpMap->value=(char*) malloc((strlen(fpath)+1)*sizeof(char));
461          sprintf(tmpMap->value,"%s",fpath);
462          pszDataSource=tmpMap->value;
463#ifdef DEBUGMS
464          fprintf(stderr,"*** DEBUG %s\n",pszDataSource);
465#endif
466        }
467      }
468      VSIFCloseL(vsif);
469      i++;
470    }
471  }
472  free(sdsName);
473  free(dsName);
474
475  OGRDataSourceH poDS = NULL;
476  OGRSFDriverH *poDriver = NULL;
477  poDS = OGROpen( pszDataSource, FALSE, poDriver );
478  if( poDS == NULL ){
479#ifdef DEBUGMS
480    fprintf(stderr,"Unable to access the DataSource %s\n",pszDataSource);
481#endif
482    setMapInMaps(conf,"lenv","message","Unable to open datasource in read only mode");
483    OGR_DS_Destroy(poDS);
484    OGRCleanupAll();
485#ifdef DEBUGMS
486    fprintf(stderr,"Unable to access the DataSource, exit! \n"); 
487#endif
488    return -1;
489  }
490
491  int iLayer = 0;
492  for( iLayer=0; iLayer < OGR_DS_GetLayerCount(poDS); iLayer++ ){
493    OGRLayerH poLayer = OGR_DS_GetLayer(poDS,iLayer);
494
495    if( poLayer == NULL ){
496#ifdef DEBUGMS
497      fprintf(stderr,"Unable to access the DataSource Layer \n");
498#endif
499      setMapInMaps(conf,"lenv","message","Unable to open datasource in read only mode");
500      return -1;
501    }
502
503    /**
504     * Add a new layer set name, data
505     */
506    if(msGrowMapLayers(m)==NULL){
507      return -1;
508    }
509    if(initLayer((m->layers[m->numlayers]), m) == -1){
510      return -1;
511    }
512
513    layerObj* myLayer=m->layers[m->numlayers];
514#ifdef DEBUGMS
515    dumpMaps(output);
516#endif
517    myLayer->name = zStrdup(output->name);
518    myLayer->tileitem=NULL;
519    myLayer->data = zStrdup(OGR_L_GetName(poLayer));
520    myLayer->connection = zStrdup(pszDataSource);
521    myLayer->index = m->numlayers;
522    myLayer->dump = MS_TRUE;
523    myLayer->status = MS_ON;
524    msConnectLayer(myLayer,MS_OGR,pszDataSource);
525
526    /**
527     * Detect the Geometry Type or use Polygon
528     */
529    if(OGR_L_GetGeomType(poLayer) != wkbUnknown){
530      switch(OGR_L_GetGeomType(poLayer)){
531      case wkbPoint:
532      case wkbMultiPoint:
533      case wkbPoint25D:
534      case wkbMultiPoint25D:
535#ifdef DEBUGMS
536        fprintf(stderr,"POINT DataSource Layer \n");
537#endif
538        myLayer->type = MS_LAYER_POINT;
539        break;
540      case wkbLineString :
541      case wkbMultiLineString :
542      case wkbLineString25D:
543      case wkbMultiLineString25D:
544#ifdef DEBUGMS
545        fprintf(stderr,"LINE DataSource Layer \n");
546#endif
547        myLayer->type = MS_LAYER_LINE;
548        break;
549      case wkbPolygon:
550      case wkbMultiPolygon:
551      case wkbPolygon25D:
552      case wkbMultiPolygon25D:
553#ifdef DEBUGMS
554        fprintf(stderr,"POLYGON DataSource Layer \n");
555#endif
556        myLayer->type = MS_LAYER_POLYGON;
557        break;
558      default:
559        myLayer->type = MS_LAYER_POLYGON;
560        break;
561      }
562    }else
563      myLayer->type = MS_LAYER_POLYGON;
564
565    /**
566     * Detect spatial reference or use WGS84
567     */
568    OGRSpatialReferenceH srs=OGR_L_GetSpatialRef(poLayer);
569    if(srs!=NULL){
570      char *wkt=NULL;
571      OSRExportToWkt(srs,&wkt);
572      setSrsInformations(output,m,myLayer,wkt);
573      free(wkt);
574    }
575    else{
576      addToMap(output->content,"crs","EPSG:4326");
577      addToMap(output->content,"crs_isGeographic","true");
578      msLoadProjectionStringEPSG(&m->projection,"EPSG:4326");
579      msInsertHashTable(&(m->web.metadata), "ows_srs", "EPSG:4326 EPSG:900913");
580      msInsertHashTable(&(myLayer->metadata), "ows_srs", "EPSG:4326 EPSG:900913");
581    }
582
583    OGREnvelope oExt;
584    if (OGR_L_GetExtent(poLayer,&oExt, TRUE) == OGRERR_NONE){
585      setMsExtent(output,m,myLayer,oExt.MinX, oExt.MinY, oExt.MaxX, oExt.MaxY);
586    }
587 
588    /**
589     * Detect the FID column or use the first attribute field as FID
590     */
591    char *fid=(char*)OGR_L_GetFIDColumn(poLayer);
592    if(strlen(fid)==0){
593      OGRFeatureDefnH def=OGR_L_GetLayerDefn(poLayer);
594      int fIndex=0;
595      for(fIndex=0;fIndex<OGR_FD_GetFieldCount(def);fIndex++){
596        OGRFieldDefnH fdef=OGR_FD_GetFieldDefn(def,fIndex);
597        fid=(char*)OGR_Fld_GetNameRef(fdef);
598        break;
599      }
600    }
601    msInsertHashTable(&(myLayer->metadata), "gml_featureid", fid);
602    msInsertHashTable(&(myLayer->metadata), "gml_include_items", "all");
603    msInsertHashTable(&(myLayer->metadata), "ows_name", output->name);
604    map* tmpMap=getMap(output->content,"title");
605    if(tmpMap!=NULL)
606      msInsertHashTable(&(myLayer->metadata), "ows_title", tmpMap->value);
607    else
608      msInsertHashTable(&(myLayer->metadata), "ows_title", "Default Title");
609   
610    if(msGrowLayerClasses(myLayer) == NULL)
611      return -1;
612    if(initClass((myLayer->CLASS[myLayer->numclasses])) == -1)
613      return -1;
614    myLayer->CLASS[myLayer->numclasses]->type = myLayer->type;
615    if(msGrowClassStyles(myLayer->CLASS[myLayer->numclasses]) == NULL)
616      return -1;
617    if(initStyle(myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]) == -1)
618      return -1;
619    /**
620     * Apply msStyle else fallback to the default style
621     */
622    tmpMap=getMap(output->content,"msStyle");
623    if(tmpMap!=NULL)
624      msUpdateStyleFromString(myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles],tmpMap->value,0);
625    else{
626      /**
627       * Set style
628       */
629      myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.red=125;
630      myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.green=125;
631      myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.blue=255;
632      myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->outlinecolor.red=80;
633      myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->outlinecolor.green=80;
634      myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->outlinecolor.blue=80;
635     
636      /**
637       * Set specific style depending on type
638       */
639      if(myLayer->type == MS_LAYER_POLYGON)
640        myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->width=3;
641      if(myLayer->type == MS_LAYER_LINE){
642        myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->width=3;
643        myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->outlinewidth=1.5;
644      }
645      if(myLayer->type == MS_LAYER_POINT){
646        myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->symbol=1;
647        myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->size=15;
648      }
649     
650    }
651    myLayer->CLASS[myLayer->numclasses]->numstyles++;
652    myLayer->numclasses++;
653   
654    m->layerorder[m->numlayers] = m->numlayers;
655    m->numlayers++;
656
657  }
658
659  OGR_DS_Destroy(poDS);
660  OGRCleanupAll();
661
662  return 1;
663}
664
665
666int tryGdal(maps* conf,maps* output,mapObj* m){
667  map* tmpMap=getMap(output->content,"storage");
668  char *pszFilename=tmpMap->value;
669  GDALDatasetH hDataset;
670  GDALRasterBandH hBand;
671  double adfGeoTransform[6];
672  int i, iBand;
673 
674  /**
675   * Try to open the DataSource using GDAL
676   */
677  GDALAllRegister();
678  hDataset = GDALOpen( pszFilename, GA_ReadOnly );
679  if( hDataset == NULL ){
680#ifdef DEBUGMS
681    fprintf(stderr,"Unable to access the DataSource %s \n",pszFilename);
682#endif
683    setMapInMaps(conf,"lenv","message","gdalinfo failed - unable to open");
684    GDALDestroyDriverManager();
685    return -1;
686  }
687#ifdef DEBUGMS
688  fprintf(stderr,"Accessing the DataSource %s %d\n",pszFilename,__LINE__);
689#endif
690
691  /**
692   * Add a new layer set name, data
693   */
694  if(msGrowMapLayers(m)==NULL){
695    return -1;
696  }
697  if(initLayer((m->layers[m->numlayers]), m) == -1){
698    return -1;
699  }
700  m->layers[m->numlayers]->index=m->numlayers;
701
702  layerObj* myLayer=m->layers[m->numlayers];
703  myLayer->name = zStrdup(output->name);
704  myLayer->tileitem=NULL;
705  myLayer->data = zStrdup(pszFilename);
706  myLayer->index = m->numlayers;
707  myLayer->dump = MS_TRUE;
708  myLayer->status = MS_ON;
709  myLayer->type = MS_LAYER_RASTER;
710
711  char *title=output->name;
712  tmpMap=getMap(output->content,"title");
713  if(tmpMap!=NULL)
714    title=tmpMap->value;
715  char *abstract=output->name;
716  tmpMap=getMap(output->content,"abstract");
717  if(tmpMap!=NULL)
718    abstract=tmpMap->value;
719
720  msInsertHashTable(&(myLayer->metadata), "ows_label", title);
721  msInsertHashTable(&(myLayer->metadata), "ows_title", title);
722  msInsertHashTable(&(myLayer->metadata), "ows_abstract", abstract);
723  msInsertHashTable(&(myLayer->metadata), "ows_rangeset_name", output->name);
724  msInsertHashTable(&(myLayer->metadata), "ows_rangeset_label", title);
725
726  /**
727   * Set Map Size to the raster size
728   */
729  m->width=GDALGetRasterXSize( hDataset );
730  m->height=GDALGetRasterYSize( hDataset );
731 
732  /**
733   * Set projection using Authority Code and Name if available or fallback to
734   * proj4 definition if available or fallback to default EPSG:4326
735   */
736  const char *tRef=GDALGetProjectionRef( hDataset );
737  if( tRef != NULL && strlen(tRef)>0 ){
738    char *pszProjection;
739    pszProjection = (char *) GDALGetProjectionRef( hDataset );
740    //#ifdef DEBUGMS
741    fprintf(stderr,"Accessing the DataSource %s\n",GDALGetProjectionRef( hDataset ));
742    //#endif
743    setSrsInformations(output,m,myLayer,pszProjection);
744  }else{
745    fprintf(stderr,"NO SRS FOUND ! %s\n",GDALGetProjectionRef( hDataset ));   
746  }
747
748
749  /**
750   * Set extent
751   */
752  if( GDALGetGeoTransform( hDataset, adfGeoTransform ) == CE_None ){
753    if( adfGeoTransform[2] == 0.0 && adfGeoTransform[4] == 0.0 ){
754
755      double minX = adfGeoTransform[0]
756        + adfGeoTransform[2] * GDALGetRasterYSize(hDataset);
757      double minY = adfGeoTransform[3]
758        + adfGeoTransform[5] * GDALGetRasterYSize(hDataset);
759
760      double maxX = adfGeoTransform[0]
761        + adfGeoTransform[1] * GDALGetRasterXSize(hDataset);
762      double maxY = adfGeoTransform[3]
763        + adfGeoTransform[4] * GDALGetRasterXSize(hDataset);
764
765       setMsExtent(output,m,myLayer,minX,minY,maxX,maxY);
766
767    }
768  }
769
770  /**
771   * Extract information about available bands to set the bandcount and the
772   * processing directive
773   */
774  char nBands[2];
775  int nBandsI=GDALGetRasterCount( hDataset );
776  sprintf(nBands,"%d",GDALGetRasterCount( hDataset ));
777  msInsertHashTable(&(myLayer->metadata), "ows_bandcount", nBands);
778  if(nBandsI>=3)
779    msLayerAddProcessing(myLayer,"BANDS=1,2,3");
780  else if(nBandsI>=2)
781    msLayerAddProcessing(myLayer,"BANDS=1,2");
782  else
783    msLayerAddProcessing(myLayer,"BANDS=1");
784
785  /**
786   * Name available Bands
787   */
788  char lBands[6];
789  char *nameBands=NULL;
790  for( iBand = 0; iBand < nBandsI; iBand++ ){
791    sprintf(lBands,"Band%d",iBand+1);
792    if(nameBands==NULL){
793      nameBands=(char*)malloc((strlen(lBands)+1)*sizeof(char));
794      sprintf(nameBands,"%s",lBands);
795    }else{
796      if(iBand<4){
797        char *tmpS=zStrdup(nameBands);
798        nameBands=(char*)realloc(nameBands,(strlen(nameBands)+strlen(lBands)+1)*sizeof(char));
799        sprintf(nameBands,"%s %s",tmpS,lBands);
800        free(tmpS);
801      }
802    }
803  }
804  msInsertHashTable(&(myLayer->metadata), "ows_bandnames", nameBands);
805 
806  /**
807   * Loops over metadata informations to setup specific informations
808   */
809  for( iBand = 0; iBand < nBandsI; iBand++ ){
810    //int         bGotNodata;//, bSuccess;
811    double      adfCMinMax[2]/*, dfNoData*/;
812    //int         nBlockXSize, nBlockYSize, nMaskFlags;
813    //double      /*dfMean, dfStdDev*/;
814    hBand = GDALGetRasterBand( hDataset, iBand+1 );
815
816    CPLErrorReset();
817    GDALComputeRasterMinMax( hBand, FALSE, adfCMinMax );
818    char tmpN[21];
819    sprintf(tmpN,"Band%d",iBand+1);
820    if (CPLGetLastErrorType() == CE_None){
821      char tmpMm[100];
822      sprintf(tmpMm,"%.3f %.3f",adfCMinMax[0],adfCMinMax[1]);
823      char tmpI[21];
824      sprintf(tmpI,"%s_interval",tmpN);
825      msInsertHashTable(&(myLayer->metadata), tmpI, tmpMm);
826
827      map* test=getMap(output->content,"msClassify");
828      if(test!=NULL && strncasecmp(test->value,"true",4)==0){
829        /**
830         * Classify one band raster pixel value using regular interval
831         */
832        int _tmpColors[10][3]={
833          {102,153,204},
834          {51,102,153},
835          {102,102,204},
836          {51,204,0},
837          {153,255,102},
838          {204,255,102},
839          {102,204,153},
840          {255,69,64},
841          {255,192,115},
842          {255,201,115}
843        };
844         
845        if(nBandsI==1){
846          double delta=adfCMinMax[1]-adfCMinMax[0];
847          double interval=delta/10;
848          double cstep=adfCMinMax[0];
849          for(i=0;i<10;i++){
850            /**
851             * Create a new class
852             */
853            if(msGrowLayerClasses(myLayer) == NULL)
854              return -1;
855            if(initClass((myLayer->CLASS[myLayer->numclasses])) == -1)
856              return -1;
857            myLayer->CLASS[myLayer->numclasses]->type = myLayer->type;
858            if(msGrowClassStyles(myLayer->CLASS[myLayer->numclasses]) == NULL)
859              return -1;
860            if(initStyle(myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]) == -1)
861              return -1;
862           
863            /**
864             * Set class name
865             */
866            char className[7];
867            sprintf(className,"class%d",i);
868            myLayer->CLASS[myLayer->numclasses]->name=zStrdup(className);
869           
870            /**
871             * Set expression
872             */
873            char expression[1024];
874            if(i+1<10)
875              sprintf(expression,"([pixel]>=%.3f AND [pixel]<%.3f)",cstep,cstep+interval);
876            else
877              sprintf(expression,"([pixel]>=%.3f AND [pixel]<=%.3f)",cstep,cstep+interval);
878            msLoadExpressionString(&myLayer->CLASS[myLayer->numclasses]->expression,expression);
879           
880            /**
881             * Set color
882             */
883            myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.red=_tmpColors[i][0];
884            myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.green=_tmpColors[i][1];
885            myLayer->CLASS[myLayer->numclasses]->styles[myLayer->CLASS[myLayer->numclasses]->numstyles]->color.blue=_tmpColors[i][2];
886            cstep+=interval;
887            myLayer->CLASS[myLayer->numclasses]->numstyles++;
888            myLayer->numclasses++;
889           
890          }
891         
892          char tmpMm[100];
893          sprintf(tmpMm,"%.3f %.3f",adfCMinMax[0],adfCMinMax[1]);
894         
895        }
896      }
897    }
898    if( strlen(GDALGetRasterUnitType(hBand)) > 0 ){
899      char tmpU[21];
900      sprintf(tmpU,"%s_band_uom",tmpN);
901      msInsertHashTable(&(myLayer->metadata), tmpU, GDALGetRasterUnitType(hBand));
902    }
903
904  }
905
906  m->layerorder[m->numlayers] = m->numlayers;
907  m->numlayers++;
908  GDALClose( hDataset );
909  GDALDestroyDriverManager();
910  CPLCleanupTLS();
911  return 1;
912}
913
914/**
915 * Create a MapFile for WMS, WFS or WCS Service output
916 */
917void outputMapfile(maps* conf,maps* outputs){
918
919  /**
920   * Firs store the value on disk
921   */
922  map* mime=getMap(outputs->content,"mimeType");
923  char *ext="data";
924  if(mime!=NULL)
925    if(strncasecmp(mime->value,"application/json",16)==0)
926      ext="json";
927
928  map* tmpMap=getMapFromMaps(conf,"main","dataPath");
929  map* sidMap=getMapFromMaps(conf,"lenv","usid");
930  char *pszDataSource=(char*)malloc((strlen(tmpMap->value)+strlen(sidMap->value)+strlen(outputs->name)+17)*sizeof(char));
931  sprintf(pszDataSource,"%s/ZOO_DATA_%s_%s.%s",tmpMap->value,outputs->name,sidMap->value,ext); 
932  int f=zOpen(pszDataSource,O_WRONLY|O_CREAT,S_IRUSR|S_IWUSR|S_IRGRP|S_IWGRP|S_IROTH|S_IWOTH);
933  map *gfile=getMap(outputs->content,"generated_file");
934  if(gfile!=NULL){
935    readGeneratedFile(conf,outputs->content,gfile->value);         
936  }
937  map* sizeMap=getMap(outputs->content,"size");
938  map* vData=getMap(outputs->content,"value");
939  if(sizeMap!=NULL){
940    zWrite(f,vData->value,atoi(sizeMap->value)*sizeof(char));
941  }
942  else{
943    zWrite(f,vData->value,(strlen(vData->value)+1)*sizeof(char));
944  }
945  close(f);
946  addToMap(outputs->content,"storage",pszDataSource);
947  free(pszDataSource);
948
949  /*
950   * Create an empty map, set name, default size and extent
951   */
952  mapObj *myMap=msNewMapObj();
953  free(myMap->name);
954  myMap->name=zStrdup("ZOO-Project_WXS_Server");
955  msMapSetSize(myMap,2048,2048);
956  msMapSetExtent(myMap,-1,-1,1,1);
957 
958  /*
959   * Set imagepath and imageurl using tmpPath and tmpUrl from main.cfg
960   */
961  map *tmp1=getMapFromMaps(conf,"main","tmpPath");
962  myMap->web.imagepath=zStrdup(tmp1->value);
963  tmp1=getMapFromMaps(conf,"main","tmpUrl");
964  myMap->web.imageurl=zStrdup(tmp1->value);
965 
966  /*
967   * Define supported output formats
968   */
969  outputFormatObj *o1=msCreateDefaultOutputFormat(NULL,"AGG/PNG","png");
970  o1->imagemode=MS_IMAGEMODE_RGBA;
971  o1->transparent=MS_TRUE;
972  o1->inmapfile=MS_TRUE;
973  msAppendOutputFormat(myMap,msCloneOutputFormat(o1));
974  msFreeOutputFormat(o1);
975
976#ifdef USE_KML
977  outputFormatObj *o2=msCreateDefaultOutputFormat(NULL,"KML","kml");
978  o2->inmapfile=MS_TRUE; 
979  msAppendOutputFormat(myMap,msCloneOutputFormat(o2));
980  msFreeOutputFormat(o2);
981#endif
982
983  outputFormatObj *o3=msCreateDefaultOutputFormat(NULL,"GDAL/GTiff","tiff");
984  if(!o3)
985    fprintf(stderr,"Unable to initialize GDAL driver !\n");
986  else{
987    o3->imagemode=MS_IMAGEMODE_BYTE;
988    o3->inmapfile=MS_TRUE; 
989    msAppendOutputFormat(myMap,msCloneOutputFormat(o3));
990    msFreeOutputFormat(o3);
991  }
992
993  outputFormatObj *o4=msCreateDefaultOutputFormat(NULL,"GDAL/AAIGRID","grd");
994  if(!o4)
995    fprintf(stderr,"Unable to initialize GDAL driver !\n");
996  else{
997    o4->imagemode=MS_IMAGEMODE_INT16;
998    o4->inmapfile=MS_TRUE; 
999    msAppendOutputFormat(myMap,msCloneOutputFormat(o4));
1000    msFreeOutputFormat(o4);
1001  }
1002
1003#ifdef USE_CAIRO
1004  outputFormatObj *o5=msCreateDefaultOutputFormat(NULL,"CAIRO/PNG","cairopng");
1005  if(!o5)
1006    fprintf(stderr,"Unable to initialize CAIRO driver !\n");
1007  else{
1008    o5->imagemode=MS_IMAGEMODE_RGBA;
1009    o5->transparent=MS_TRUE;
1010    o5->inmapfile=MS_TRUE;
1011    msAppendOutputFormat(myMap,msCloneOutputFormat(o5));
1012    msFreeOutputFormat(o5);
1013  }
1014#endif
1015
1016  /*
1017   * Set default projection to EPSG:4326
1018   */
1019  msLoadProjectionStringEPSG(&myMap->projection,"EPSG:4326");
1020  myMap->transparent=1;
1021
1022  /**
1023   * Set metadata extracted from main.cfg file maps
1024   */
1025  maps* cursor=conf;
1026  map* correspondance=getCorrespondance();
1027  while(cursor!=NULL){
1028    map* _cursor=cursor->content;
1029    map* vMap;
1030    while(_cursor!=NULL){
1031      if((vMap=getMap(correspondance,_cursor->name))!=NULL){
1032        if (msInsertHashTable(&(myMap->web.metadata), vMap->value, _cursor->value) == NULL){
1033#ifdef DEBUGMS
1034          fprintf(stderr,"Unable to add metadata");
1035#endif
1036          return;
1037        }
1038      }
1039      _cursor=_cursor->next;
1040    }
1041    cursor=cursor->next;
1042  }
1043  freeMap(&correspondance);
1044  free(correspondance);
1045
1046  /*
1047   * Set mapserver PROJ_LIB/GDAL_DATA or any other config parameter from
1048   * the main.cfg [mapserver] section if any
1049   */
1050  maps *tmp3=getMaps(conf,"mapserver");
1051  if(tmp3!=NULL){
1052    map* tmp4=tmp3->content;
1053    while(tmp4!=NULL){
1054      msSetConfigOption(myMap,tmp4->name,tmp4->value);
1055      tmp4=tmp4->next;
1056    }
1057  }
1058
1059  /**
1060   * Set a ows_rootlayer_title, 
1061   */
1062  if (msInsertHashTable(&(myMap->web.metadata), "ows_rootlayer_name", "ZOO_Project_Layer") == NULL){
1063#ifdef DEBUGMS
1064    fprintf(stderr,"Unable to add metadata");
1065#endif
1066    return;
1067  }
1068  if (msInsertHashTable(&(myMap->web.metadata), "ows_rootlayer_title", "ZOO_Project_Layer") == NULL){
1069#ifdef DEBUGMS
1070    fprintf(stderr,"Unable to add metadata");
1071#endif
1072    return;
1073  }
1074
1075  /**
1076   * Enable all the WXS requests using ows_enable_request
1077   * see http://mapserver.org/trunk/development/rfc/ms-rfc-67.html
1078   */
1079  if (msInsertHashTable(&(myMap->web.metadata), "ows_enable_request", "*") == NULL){
1080#ifdef DEBUGMS
1081    fprintf(stderr,"Unable to add metadata");
1082#endif
1083    return;
1084  }
1085  msInsertHashTable(&(myMap->web.metadata), "ows_srs", "EPSG:4326");
1086
1087  if(tryOgr(conf,outputs,myMap)<0)
1088    if(tryGdal(conf,outputs,myMap)<0)
1089      return ;
1090
1091  tmp1=getMapFromMaps(conf,"main","dataPath");
1092  char *tmpPath=(char*)malloc((13+strlen(tmp1->value))*sizeof(char));
1093  sprintf(tmpPath,"%s/symbols.sym",tmp1->value);
1094  msInitSymbolSet(&myMap->symbolset);
1095  myMap->symbolset.filename=zStrdup(tmpPath);
1096  free(tmpPath);
1097
1098  map* sid=getMapFromMaps(conf,"lenv","usid");
1099  char *mapPath=
1100    (char*)malloc((7+strlen(sid->value)+strlen(outputs->name)+strlen(tmp1->value))*sizeof(char));
1101  sprintf(mapPath,"%s/%s_%s.map",tmp1->value,outputs->name,sid->value);
1102  msSaveMap(myMap,mapPath);
1103  free(mapPath);
1104  msGDALCleanup();
1105  msFreeMap(myMap);
1106}
1107
1108#endif
Note: See TracBrowser for help on using the repository browser.

Download now !

Search

Context Navigation

Main Navigation

ZOO Sponsors

http://www.zoo-project.org/trac/chrome/site/img/geolabs-logo.pnghttp://www.zoo-project.org/trac/chrome/site/img/neogeo-logo.png http://www.zoo-project.org/trac/chrome/site/img/apptech-logo.png http://www.zoo-project.org/trac/chrome/site/img/3liz-logo.png http://www.zoo-project.org/trac/chrome/site/img/gateway-logo.png

Become a sponsor !

Knowledge partners

http://www.zoo-project.org/trac/chrome/site/img/ocu-logo.png http://www.zoo-project.org/trac/chrome/site/img/gucas-logo.png http://www.zoo-project.org/trac/chrome/site/img/polimi-logo.png http://www.zoo-project.org/trac/chrome/site/img/fem-logo.png http://www.zoo-project.org/trac/chrome/site/img/supsi-logo.png http://www.zoo-project.org/trac/chrome/site/img/cumtb-logo.png

Become a knowledge partner

Related links

http://zoo-project.org/img/ogclogo.png http://zoo-project.org/img/osgeologo.png