svn-gvsig-desktop / tags / v1_1_Build_1001 / extensions / extGeoProcessing / src / com / iver / cit / gvsig / geoprocess / core / util / UnitUtils.java @ 11984
History | View | Annotate | Download (4.6 KB)
1 | 6108 | azabala | /*
|
---|---|---|---|
2 | * Created on 29-jun-2006
|
||
3 | *
|
||
4 | * gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana
|
||
5 | *
|
||
6 | * Copyright (C) 2004 IVER T.I. and Generalitat Valenciana.
|
||
7 | *
|
||
8 | * This program is free software; you can redistribute it and/or
|
||
9 | * modify it under the terms of the GNU General Public License
|
||
10 | * as published by the Free Software Foundation; either version 2
|
||
11 | * of the License, or (at your option) any later version.
|
||
12 | *
|
||
13 | * This program is distributed in the hope that it will be useful,
|
||
14 | * but WITHOUT ANY WARRANTY; without even the implied warranty of
|
||
15 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
||
16 | * GNU General Public License for more details.
|
||
17 | *
|
||
18 | * You should have received a copy of the GNU General Public License
|
||
19 | * along with this program; if not, write to the Free Software
|
||
20 | * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,USA.
|
||
21 | *
|
||
22 | * For more information, contact:
|
||
23 | *
|
||
24 | * Generalitat Valenciana
|
||
25 | * Conselleria d'Infraestructures i Transport
|
||
26 | * Av. Blasco Ib??ez, 50
|
||
27 | * 46010 VALENCIA
|
||
28 | * SPAIN
|
||
29 | *
|
||
30 | * +34 963862235
|
||
31 | * gvsig@gva.es
|
||
32 | * www.gvsig.gva.es
|
||
33 | *
|
||
34 | * or
|
||
35 | *
|
||
36 | * IVER T.I. S.A
|
||
37 | * Salamanca 50
|
||
38 | * 46005 Valencia
|
||
39 | * Spain
|
||
40 | *
|
||
41 | * +34 963163400
|
||
42 | * dac@iver.es
|
||
43 | */
|
||
44 | /* CVS MESSAGES:
|
||
45 | *
|
||
46 | * $Id$
|
||
47 | * $Log$
|
||
48 | 10552 | azabala | * Revision 1.3.2.5 2007-02-28 12:26:01 azabala
|
49 | * solved bug when view has a geodetic projection
|
||
50 | *
|
||
51 | * Revision 1.3.2.4 2007/01/12 13:25:53 caballero
|
||
52 | 9711 | caballero | * unidades de mapa y medida
|
53 | *
|
||
54 | * Revision 1.3.2.3 2006/12/22 10:15:07 caballero
|
||
55 | 9450 | caballero | * unidades de mapa y de medida
|
56 | *
|
||
57 | * Revision 1.3.2.2 2006/12/21 18:21:52 azabala
|
||
58 | 9425 | azabala | * Bug fixed: change getMapUnits for getDistanceUnits
|
59 | *
|
||
60 | * Revision 1.3.2.1 2006/11/15 00:08:17 jjdelcerro
|
||
61 | 8765 | jjdelcerro | * *** empty log message ***
|
62 | *
|
||
63 | * Revision 1.4 2006/09/15 10:42:54 caballero
|
||
64 | * extensibilidad de documentos
|
||
65 | *
|
||
66 | * Revision 1.3 2006/08/29 07:56:30 cesar
|
||
67 | 6880 | cesar | * Rename the *View* family of classes to *Window* (ie: SingletonView to SingletonWindow, ViewInfo to WindowInfo, etc)
|
68 | *
|
||
69 | * Revision 1.2 2006/08/29 07:21:09 cesar
|
||
70 | 6878 | cesar | * Rename com.iver.cit.gvsig.fmap.Fmap class to com.iver.cit.gvsig.fmap.MapContext
|
71 | *
|
||
72 | * Revision 1.1 2006/06/29 17:30:36 azabala
|
||
73 | 6108 | azabala | * first version in cvs
|
74 | *
|
||
75 | *
|
||
76 | */
|
||
77 | package com.iver.cit.gvsig.geoprocess.core.util; |
||
78 | |||
79 | import org.cresques.cts.IProjection; |
||
80 | import org.cresques.cts.gt2.CSUTM; |
||
81 | |||
82 | import com.iver.andami.PluginServices; |
||
83 | 6878 | cesar | import com.iver.cit.gvsig.fmap.MapContext; |
84 | 8765 | jjdelcerro | import com.iver.cit.gvsig.project.documents.view.gui.View; |
85 | 6108 | azabala | |
86 | /**
|
||
87 | * Utility methods to manage distances
|
||
88 | * (conversion from user entries to internal units,
|
||
89 | 9450 | caballero | * verifying if projection is planar or we are
|
90 | 6108 | azabala | * working with geographic coordinates, etc.)
|
91 | * @author azabala
|
||
92 | *
|
||
93 | */
|
||
94 | public class UnitUtils { |
||
95 | 9450 | caballero | |
96 | 6108 | azabala | public static final double EARTH_RADIUS = 6378137d; |
97 | |||
98 | /**
|
||
99 | * For computing with geodetic coordinates:
|
||
100 | * returns the angular measure (in radians)
|
||
101 | 9450 | caballero | * for a distance over the earth along a
|
102 | 6108 | azabala | * meridiam.
|
103 | * Because this consideration is an approximation,
|
||
104 | * we consideer the eart like an sphere (not an
|
||
105 | * ellipsoid)
|
||
106 | * @param dist distance in meters
|
||
107 | * @return arc of meridian whose length is the specified
|
||
108 | * distance
|
||
109 | */
|
||
110 | public static double toSexaAngularMeasure(double dist){ |
||
111 | /*
|
||
112 | * dist = 6378km(terrestrial radius) -> 2PI
|
||
113 | * passed distance -> incognite
|
||
114 | */
|
||
115 | return Math.toDegrees((2 * Math.PI * dist)/EARTH_RADIUS); |
||
116 | } |
||
117 | 9450 | caballero | |
118 | |||
119 | 6108 | azabala | |
120 | 9425 | azabala | /**
|
121 | * Converts a distance entered by user in the GUI in the same distance
|
||
122 | * in internal units (measure units)
|
||
123 | *
|
||
124 | * */
|
||
125 | 6108 | azabala | public static double getInInternalUnits(double userEntryDistance){ |
126 | 10552 | azabala | // VCN he modificado esto a como creo que deber?a de estar,
|
127 | //as? tiene en cuenta para calcular la distancia tanto las unidades
|
||
128 | //en las que se encuentra la cartograf?a como las unidades de medida
|
||
129 | //seleccionadas por el usuario.
|
||
130 | IProjection proj = ((View)PluginServices.
|
||
131 | getMDIManager(). |
||
132 | getActiveWindow()). |
||
133 | getMapControl(). |
||
134 | getViewPort(). |
||
135 | getProjection(); |
||
136 | |||
137 | 9450 | caballero | int distanceUnits = ((View)PluginServices. |
138 | getMDIManager(). |
||
139 | getActiveWindow()).getMapControl().getViewPort().getDistanceUnits(); |
||
140 | 10552 | azabala | |
141 | boolean isProjected = proj.isProjected();
|
||
142 | int mapUnits = -1; |
||
143 | if(isProjected){
|
||
144 | mapUnits = ((View)PluginServices.
|
||
145 | 6108 | azabala | getMDIManager(). |
146 | 9450 | caballero | getActiveWindow()).getMapControl().getViewPort().getMapUnits(); |
147 | 10552 | azabala | }else{
|
148 | mapUnits = 1;
|
||
149 | } |
||
150 | |||
151 | double distInInternalUnits = (userEntryDistance/MapContext.CHANGEM[mapUnits]) *
|
||
152 | MapContext.CHANGEM[distanceUnits]; |
||
153 | |||
154 | 6108 | azabala | /*
|
155 | * if layer's projections is in geographics
|
||
156 | coords, pass distance to a angular measure
|
||
157 | */
|
||
158 | 10552 | azabala | |
159 | if( (proj != null) && !(isProjected)){ |
||
160 | 6108 | azabala | distInInternalUnits = |
161 | 9450 | caballero | UnitUtils.toSexaAngularMeasure(distInInternalUnits); |
162 | 6108 | azabala | } |
163 | return distInInternalUnits;
|
||
164 | } |
||
165 | 9450 | caballero | |
166 | |||
167 | |||
168 | |||
169 | |||
170 | |||
171 | 6108 | azabala | } |