Discussion :

[zenico64]

Bonjour à tous,

Je développe en WPF avec VisualStudio 2010 sur Windows 7 32bits.

Je réalise un programme dans lequel je dois faire du géo-processing.
Pour ce faire, je me suis tournée vers la solution GDAL.

J'ai donc télécharger la SDK à cet endroit :
http://www.gisinternals.com/sdk/





J'ai choisi les solutions surlignées


Je me retrouve donc avec les fichiers
gdal-111-1400-core.msi
release-1400-gdal-1-11-mapserver-6-4.zip

J'ai exécuté gdal-111-1400-core.msi et choisis une installation typique.
ce qui m'a créé le répertoire GDAL dans programfiles


J'extrait ensuite le fichier compressé release-1400-gdal-1-11-mapserver-6-4.zip dans le répertoire de l'instalation


une fois cela fait, j'ajoute à ma variable d'environnement path les valeurs :
;C:\Program Files\GDAL\csharp;C:\Program Files\GDAL

(Je ne me souviens plus d'où j'ai tiré cette procédure d'installation)

J'ajoute ensuite la bibliothèque ogr_csharp dans les références de mon projet pour pouvoir y faire appel.

Code C# :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
using OSGeo.OGR;
//...
var wkt1 = getWktFromArcGISObject(polygon1);
var geom1 = OSGeo.OGR.Geometry.CreateFromWkt(wkt1);
 
var wkt2 = getWktFromArcGISObject(polygon2);
var geom2 = OSGeo.OGR.Geometry.CreateFromWkt(wkt2);
 
return geom2.Intersect(geom1);
//...

Là mon application fonctionne à merveille et tout se passe bien.

Quand je la déploie sur un autre ordinateur avec un système 32bits, en répétant la procédure d’installation de GDAL tout se passe également à merveille.

Le problème survient quand j'essaie de déployer l'application sur un système 64bits. J'ai essayé plusieurs manipulations :

• Installer GDAL sur le système ciblé en 64 bits
• compiler l'application avec un .dll tiré d'un installation GDAL en 64 bits
• installer GDAL sur le système ciblé en 32 bits
• Autres essaies dont je ne me souviens pas la nature...

Rien ne marche et j'ai obtenus ces différents messages d'erreurs (pas dans l'ordre)
(dans le prochains message puisqu'il est interdit de joindre plus de 5 fichiers dans un message)

J'ai également vu quelque part ce message :

It depends on the platform target settings and the your current framework,
which kind of binaries should be used, see the table below:

Platform target Executing environment (Framework version) binary package should be used
/x86 x86 x86
/x86 x64 x86
/x64 x64 x64
/AnyCPU x86 x86
/AnyCPU x64 x64

All other combinations are not supported, and will probably result in
BadImageFormatException.

Peut être que l'installation que je fais de GDAL est une mauvaise solution. C'est cependant la seule solution que j'ai réussi à mettre en place.
Il est possible également que la solution réside dans la compilation de l'application sur un système 64 bits (Malheureusement je ne dispose pas d'un tel environnement de travail)

Merci à vous !

[zenico64]

[zenico64]

Afin de contourner mon problème, j'ai utilisé un autre moyen pour réaliser l'intersection entre deux formes géométriques :

avec System.Windows.Media

Code C# :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
 
        /// <summary>
        /// Vérifie avec exactitude mais de façon assez lente si deux polygones se chevauchent
        /// </summary>
        /// <param name="polygon1"></param>
        /// <param name="line"></param>
        /// <returns></returns>
        private static bool checkCompleteIntersectBetweenArcGISPolygonAndArcGISPolyline(ESRI.ArcGIS.Client.Geometry.Polygon polygon1, ESRI.ArcGIS.Client.Geometry.Polyline line)
        {
            var g1 = StreamGeometry.Parse(getPathFromArcGISPolygon(polygon1));
            var g2 = LineGeometry.Parse(getPathFromArcGISPolyLine(line));
            Geometry og2 = g2.GetWidenedPathGeometry(new Pen(Brushes.Black, 1.0));
 
            CombinedGeometry cg = new CombinedGeometry(GeometryCombineMode.Intersect, g1, og2);
            PathGeometry pg = cg.GetFlattenedPathGeometry();
            return pg.Figures.Count > 0;
        }
 
        /// <summary>
        /// Vérifie avec exactitude mais de façon assez lente si deux polygones se chevauchent
        /// </summary>
        /// <param name="polygon1"></param>
        /// <param name="polygon2"></param>
        /// <returns></returns>
        private static bool checkCompleteIntersectWithinArcGISPolygon(ESRI.ArcGIS.Client.Geometry.Polygon polygon1, ESRI.ArcGIS.Client.Geometry.Polygon polygon2)
        {
            var g1 = StreamGeometry.Parse(getPathFromArcGISPolygon(polygon1));
            var g2 = StreamGeometry.Parse(getPathFromArcGISPolygon(polygon2));
 
            CombinedGeometry cg = new CombinedGeometry(GeometryCombineMode.Intersect, g1, g2);
            PathGeometry pg = cg.GetFlattenedPathGeometry();
            return pg.Figures.Count > 0;
        }
 
        /// <summary>
        /// 
        /// </summary>
        /// <param name="polyLine"></param>
        /// <returns></returns>
        private static string getPathFromArcGISPolyLine(ESRI.ArcGIS.Client.Geometry.Polyline polyLine)
        {
            var provider = new CultureInfo("en");
            var path = "";
            foreach (var pathArcGis in polyLine.Paths)
            {
                var first = true;
                foreach (var pointArcGis in pathArcGis)
                {
                    if (first)
                    {
                        path += string.Format("M {0},{1}", pointArcGis.X.ToString(provider), pointArcGis.Y.ToString(provider));
                        first = false;
                        continue;
                    }
                    path += string.Format(" L {0},{1}", pointArcGis.X.ToString(provider), pointArcGis.Y.ToString(provider));
                }
            }
            return path;
        }
 
        /// <summary>
        /// 
        /// </summary>
        /// <param name="polygon1"></param>
        /// <returns></returns>
        private static string getPathFromArcGISPolygon(ESRI.ArcGIS.Client.Geometry.Polygon polygon1)
        {
            var provider = new CultureInfo("en");
            var path = "";
            foreach (var ring in polygon1.Rings)
            {
                var first = true;
                foreach (var pointArcGis in ring)
                {
                    if (first)
                    {
                        path += string.Format("M {0},{1}", pointArcGis.X.ToString(provider), pointArcGis.Y.ToString(provider));
                        first = false;
                        continue;
                    }
                    path += string.Format(" L {0},{1}", pointArcGis.X.ToString(provider), pointArcGis.Y.ToString(provider));
                }
                path += "z";
            }
            return path;
        }

Je n'ai donc pas vraiment résolu le problème mais je l'ai contourné...