在Mollweide上将lat，long转换为x，y
在Mollweide投影中，将纬度（lat）和经度（long）转换为x和y坐标的过程如下：
- Mollweide投影概念：Mollweide投影是一种等面积的伪圆柱投影，用于将地球表面的经纬度坐标转换为平面坐标。
- 转换公式：Mollweide投影的转换公式如下： x = (2 * sqrt(2) * long * cos(lat)) / pi y = sqrt(2) * sin(lat)
- 其中，x和y分别表示Mollweide投影平面上的坐标，lat表示纬度，long表示经度，pi表示圆周率。
- 分类和优势：Mollweide投影属于伪圆柱投影的一种，其优势在于保持了地球表面的面积比例，适用于需要保持地区面积比例的地图制作和分析。
- 应用场景：Mollweide投影常用于制作世界地图、气候分布图、人口密度图等需要保持面积比例的地理信息可视化应用。
- 腾讯云相关产品和产品介绍链接地址：腾讯云提供了丰富的云计算服务和解决方案，但在本回答中不涉及具体的产品推荐和链接地址。您可以访问腾讯云官方网站（https://cloud.tencent.com/）了解更多相关产品和服务信息。
请注意，以上回答仅涵盖了将lat和long转换为x和y坐标的基本概念和公式，如果需要更详细的解释或其他相关信息，请提供更具体的问题。
相关·内容
我正在编写一个程序，该程序将从用户获取一组经纬度坐标，将其转换为Mollweide投影地图的x&y坐标，然后报告这些坐标下像素的值(在本例中是噪声温度)。
我使用的地图/数据是，它是作为Mollweide投影图提供的。该数据采用.fits格式，是对408 MHz波段噪声的一次全天大测量。
根据，可以使用Newton迭代将经度/纬度转换为x/y地图坐标。我的程序中的迭代方案主要是基于维基百科页面和中的方法。
但是，我的程序似乎没有报告我输入的经度和纬度的正确值。我很大程度上怀疑这两种(或两种)因素中的一种是导致这个错误的原因：
我实现迭代方案的方式是不正确的，从而导致不正确的值被报告
浏览 2提问于2017-07-14得票数 1
我有一个3k +数据点分布在墨西哥湾北部的数据框架(这里我只提供6个)。我正在尝试为到海岸的距离(公里)创建一个新的变量。我有一个形状文件(gulf.shape)，我想使用，但我不清楚如何使用。
以下是一些数据
require(maptools)
require(sp)
library(rgdal)
library(lubridate)
df <- data.frame(Lat = c(26.84853, 28.38329, 28.00364,
29.53840, 29.32030, 26.81622, 25.28146),
浏览 1提问于2018-07-31得票数 1
1回答
我开始了一个“免费”开源项目，为地球海洋的pH创建一个新的数据集。
我从NOAA的开放数据集开始，用这些列创建了一个245万行数据集：
colnames(NOAA_NODC_OSD_SUR_pH_7to9)
[1] "Year" "Month" "Day" "Hour" "Lat" "Long" "Depth" "pH"
方法文档。
数据集.
我现在的目标是“合格”每一行(2.45米).为此，我需要计算Lat/Long的每个点到最近的海岸的距离。
浏览 4提问于2014-12-29得票数 6
我想根据一个名为'Dataset‘的变量为添加到地图的点设置自定义形状、大小和颜色。如果我将所有点的形状设置为相同的类型，我可以设置点的颜色，但我希望有一个包含更多信息的地图。当我运行这段代码时，所有的点都是黑色的圆圈。我遗漏了什么？
感谢大家的帮助和时间！！
下面是一个可重复使用的示例：
# Read in libraries
library(ggplot2)
library(maps)
library(maptools)
library(ggmap)
# Create mapping objects
world <-
浏览 18提问于2018-01-24得票数 1
回答已采纳
我正在尝试将下面代码的输出(这是一个的list列表)转换为一个数据格式，其中的列包含每个观察的简单功能集合 -one多边形。
以下是可重复例子的产出和预期产出：
# Reproducible example:
library(tidyverse)
library(sf)
library(magrittr)
# define radius for circle
radius <- 40
r <- units::set_units(radius, units::as_units("nmile"), mode = "standard"
浏览 1提问于2018-02-27得票数 7
回答已采纳
我有这个df：
d = {'name':'CompanyABCD' ,
'office_location':[{'office_x':'lat,long','office_y':'lat,long'}] ,
'total_employees':100}
test = pd.DataFrame(d)
测试结果如下：
name office_location tota
浏览 3提问于2020-04-17得票数 0
回答已采纳
我想制作一个带有标签的地图，但是我不想要多边形中的labes。有没有办法把标签放在多边形上？ 下面是一个可重现的例子： require(maps)
require(ggplot2)
require(ggrepel)
set.seed(1)
world <- map_data("world")
world <- world[world$region != "Antarctica",]
cities <- structure(list(city = c("Dallas", "Johannesburg", "
浏览 16提问于2020-06-24得票数 2
回答已采纳
2回答
我是地理空间python库和距离计算的新手，我想避免for循环，这可能是非常昂贵的。
我有纬度和经度的列表，我需要沿着轨迹顺序地创建另一个固定dist=5处的点的列表。
list_x = []
list_y = []
r_val = 5
long[0] = long_start
lat[0] = lat_start
for i,j in zip(range(len(long)), range(len(lat))):
for i,j in zip(range(len(long)), range(len(lat))):
long_cent = long[i]
浏览 0提问于2021-03-12得票数 0
1回答
我尝试计算一个节点在OpenStreetMap中的给定瓦片，这样我就可以手动将节点直接添加到current_nodes表中。如果我想在纬度X和经度Y处添加一个新节点，那么在插入到current_nodes表之前，我可以计算出所需的“平铺”值。 我正在尝试将当前的PostgreSQL函数转换为C#，但是它没有给出预期的输出。我不确定我是否正确地转换了它，或者我是否遗漏了什么。 //maptile_for_point(scaled_lat bigint, scaled_lon bigint, zoom integer)
DECLARE
lat CONSTANT DOUBLE PRECISIO
浏览 27提问于2020-10-21得票数 0
1回答
我有在R中创建球体的纬度和经度，我需要做的是根据相应的振幅给球体上色。在excel文件中有3列: lat、long和amplitude。所以每一行都是一个数据点。我得到的是一个统一的模式-它不能从振幅列中读取值。
test$x<-0.5*cos(test$lat)*cos(test$long)
test$y<-0.5*cos(test$lat)*sin(test$long)
test$z<-0.5*sin(test$lat)
test$color<-test$amplitude
ramp<- colorRamp(c("blue", "red
浏览 4提问于2019-08-23得票数 0
1回答
将csv ( 50k行)转换为json以供Django模板最终使用是相当慢的。我想知道我是否正确地转换了它，或者是否有更好的方法来做到这一点。
csv的前几行是：
tdate,lat,long,entity
3/6/2017,34.152568,-118.347831,x1
6/3/2015,34.069787,-118.384738,y1
1/1/2011,34.21377,-118.227904,x1
3/4/2013,33.81761,-118.070374,y1
我在视图中读取此csv，并以这种方式呈现请求：
def index(request):
df = pd.r
浏览 10提问于2018-01-05得票数 3
回答已采纳
我找到了一种使用ggplot2绘制Voronoi图解段的方法：
library(deldir)
library(ggplot2)
library(ggthemes)
set.seed(123)
df <- data.frame(lat = rnorm(20,39,10),long = rnorm(20,-98,15))
voronoi <- deldir(df$long, df$lat)
ggplot(data=df, aes(x=long,y=lat)) +
geom_segment(aes(x = x1, y = y1, xend = x2, yend = y2),siz
浏览 1提问于2014-06-19得票数 2
回答已采纳
1回答
我正在尝试使用vbscript获取从一个纬度/经度到另一个纬度/经度的rhumb线的标题。这类似于
我还试图将中的javascript代码转换为vbscript。
Javascript
var dPhi = Math.log(Math.tan(lat2/2+Math.PI/4)/Math.tan(lat1/2+Math.PI/4));
var q = dLat/dPhi;
if (!isFinite(q)) q = Math.cos(lat1.toRad());
// if dLon over 180° take shorter rhumb across 180° meridian:
if (
浏览 3提问于2014-01-30得票数 0
2回答
我有一张ggplot的地图(由geom_polygon、线、点和段组成)。原版为WGS，所有单位均为零下/伦度。所有额外的数据(对于线和点)也是以度为单位的。我想把整件事重新转换成等面积投影(Mollweide)。第一步(转换底层的rworldmap多边形)是相当无痛的，尽管有一些奇怪的水平线贯穿其中。但是，我不知道下一步该怎么做。正如我所说的，所有数据都是以度数表示的，我希望以度数来指定xlim/ylim、轴标签和(曲线)网格线。我是否需要把所有的东西都转换成Mollweide表，或者是否有可能在一个无痛的步骤中重新塑造我的最终地图对象？
到目前为止，这就是我要展示的奇怪的水平线。
libr
浏览 3提问于2015-01-12得票数 1
回答已采纳
我花了很多时间上网想办法来做这件事。通过exist文档，我编写了下面的代码，但是它没有显示正确的角度。有什么问题吗？
private double angleFromCoordinate(double lat1, double long1, double lat2,
double long2) {
double dLon = (long2 - long1);
double y = Math.sin(dLon) * Math.cos(lat2);
double x = Math.cos(lat1) * Math.s
浏览 2提问于2015-03-24得票数 2
回答已采纳
1回答
目前，我正在将我的熊猫数据框架中的纬度、经度坐标投影到笛卡尔平面上。因此，我有一种投影方法：
def convert_lat_long_xy(lat, lo):
return x, y
因此，它返回一个元组，我可以在我的dataframe上使用这个方法，如下所示：
df.apply(lambda x: convert_lat_long_xy(x.latitude, x.longitude), axis=1))
现在，我想要做的是在我的数据框架中创建两个额外的列，名为'x‘和'y’，以保存这些值。我知道我可以做这样的事：
df['proj'] = df.a
浏览 0提问于2019-05-10得票数 2
回答已采纳
1回答
我有两张桌子fdone, incident_one.
fdone包含四种颜色。
INC_KEY,LAT,LONG,FDID.
incident_one包含five颜色。
INC_KEY,LAT,LONG,STATE,FDID.
通过这个查询，我使fdone表LAT长2为4。
SELECT S.`LAT `,S.`LONG `,(S.`LAT `+5)as x,(S.`LONG `-5)as Y FROM `fdheader` S
然后我得到了
LAT，LONG，X，Y
然后使用联接查询。
我得检查一下，
incident_one LAT和LONG在LAT，LONG，X，Y.
我得查一下.，
浏览 2提问于2013-03-25得票数 0
我编写了一个函数来从一个lat, long服务器返回地址的ArcGIS。职能如下：
def get_lat_long(address):
coordinates = geocode_address(address)
lat = "{:.2f}".format(coordinates['y'])
long = "{:.2f}".format(coordinates['x'])
print(lat, long) # prints: 42.34 -71.15
return lat, long
正
浏览 2提问于2019-03-17得票数 1
回答已采纳
我正在学习matplotlib和映射函数Basemap，并试图在地图上绘制一个简单的lat/long坐标列表。然而，Basemap坐标转换为简单的x/y坐标是在地图上绘制的点，我不知道为什么。我的代码如下：
locationData = parseFile(locationFile)
fig = plt.figure(figsize=(8,8))
m = Basemap(projection='merc', resolution='h',
llcrnrlat=49, llcrnrlon=-13.5,
urcrnr
浏览 0提问于2018-03-07得票数 4
回答已采纳
1回答
我正在使用matplotlib basemap从熊猫数据中绘制一个轨迹。
简单地绘制lat、lon数据，如下所示：
map = Basemap(projection='stere',lon_0=lon,lat_0=lat, llcrnrlon=-65, llcrnrlat=55, urcrnrlon=50 , urcrnrlat=75)
map.drawcountries()
map.bluemarble()
#Place marker at Kangerlussuaq
x,y = map(lon, lat)
map.plot(x, y, 'bo', mar
浏览 0提问于2018-05-09得票数 2
回答已采纳
1回答
Python请求API
、、、、
我试图从NOAA API中检索一些数据，但是有一个错误我无法解决
location=[]
def find_xy(Name, lat, long):
api = url+str(lat)+','+str(long)
r = requests.get(api).json()
x = r['properties']['gridX']
y = r['properties']['gridY']
xy=(Name, str(lat), str(long), x, y)
浏览 1提问于2020-11-06得票数 0
1回答
我有一些数据如下：
>show(recruitment_info)
Centre Lat Long GroupA GroupB
1 CentreA 51.51770 -0.100400 907 47
2 CentreB 52.48947 -1.898575 1910 116
3 CentreC 51.45451 -2.587910 4419 277
我要绘制一张英国地图，然后为每个中心加分(在第1栏中用中心标记)。我还想要点大小来表示GroupA
浏览 2提问于2019-09-10得票数 2
回答已采纳
希望这是我天真地忽略的东西，但我不明白为什么我的pcolor图没有显示我的任何数据。我有2D lat/lon，一个时间步长为2D海冰浓度。任何提示都将不胜感激，这就是我在下面所尝试的：
from netCDF4 import Dataset
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
# ifile='upper_box.nc' # doesn't produce any plot
ifile='1979_sfc_o
浏览 0提问于2017-04-25得票数 2
回答已采纳
我正在尝试删除数据文件中的重复元素。这个DataFrame来自于计算给定的地理坐标列表之间的距离。正如您在下面的DataFrame中所看到的，数据是重复的，但我不能将索引设置为“dist”，因为在其他情况下，距离可能是0或1(重复)，然后重要数据将被丢弃。
import pandas as pd
df = pd.DataFrame({'Name_x':['a','b','c','d'],
'Name_y':['b','a','d'
浏览 0提问于2017-10-16得票数 0
回答已采纳
我有来自的地名的数据转储(csv)
我需要将它导入mysql，但是几何协调使用BNG (OSGB36)。Mysql是否有将这些协调器转换为wgs84 lat/long的功能，或者是否有其他的sql方法来实现这一点？
另一种选择可能是将其加载到postgis - postgis是否具有将BNG转换为lat/long的功能？也许我可以这样做，然后导出数据，并将其加载到mysql中？
在后地理信息系统中，我可以这样做(没有测试)
select AddGeometryColumn('locations', 'the_geom', 27700, 'POINT
浏览 1提问于2019-01-04得票数 0
我试图在一个开放的街道地图瓷砖上重叠点(来自dataframe df)。带点的图层虽然掩盖了OpenStreetMap地图的平铺。
library(OpenStreetMap)
library(tidyverse)
map <- openmap(c(lat = 46.2037+ 0.01, lon = 6.1439 - 0.01), c(lat = 46.2037 -
0.005, lon = 6.1439 + 0.01), type = "osm")
mapLatLon <- openproj(map)
df <- data.frame(&#
浏览 3提问于2019-12-02得票数 0
回答已采纳
我在将纬度和经度坐标转换为适合绘制到JPanel的X/Y坐标时遇到了困难。我可以将它们转换为适当的x/y值，但是这些值都非常相似，而且不足以真正绘制我想要绘制的地图。我使用的是Line2D的JPanel，但是我看到了使用多边形的例子，但是当我这样做的时候，地图看起来一点也不正确。
我转换的lat/long点都非常接近(都位于大学校园内)，所以我的x/y值也非常接近。我需要一种方法来解析它们，这样它们就可以用600x600像素的JPanel来表示。下面是我转换为XY及其值的一些要点的示例：
GL 43.130480 -77.631559
i2 43.130425 -77.631532
浏览 0提问于2016-04-20得票数 2
回答已采纳
<?php
header("Content-type:image/png");
$image1 = imagecreatefrompng('cylindrical_map.png');
$image2 = imagecreatefrompng('pin.png');
$long = 18.96;
$lat = 72.82;
$scale_x = imagesx($image1);
$scale_y = imagesy($image1);
$pt = getlocationcoords($lat, $long, $scale_x, $s
浏览 1提问于2013-12-31得票数 0
回答已采纳
2回答
我正在尝试将一个城市的经度/纬度位置转换为笛卡尔坐标，这样我就可以将其转换为像素位置，以便在pygame屏幕上显示，但我在弄清楚它时遇到了麻烦。
def ConvertPoint(self, Long, Lat, height, width):
self.x = Long
self.y = Lat
self.x = angleToPointX(self.x, self.y, 3959)
self.y = angleToPointY(self.x, self.y, 3959)
pygame.draw.circle(self.window, (255,0,0
浏览 0提问于2012-04-10得票数 2
回答已采纳
1回答
好的。所以我想应该是在谷歌地球上高亮显示纽约和北京…
我正在尝试弄清楚如何将点映射到3d基本体对象(球体)上，通过圆周通过任意方向获得它们的距离，通过直径获得它们的距离。这些点将是纬度和经度坐标。
现在，这就是我试图用来映射坐标(一个代码不可知的版本)：
x1 = radius * cos(long1) * cos(lat1);
y1 = radius * sin(long1) * cos(lat1);
z1 = radius * sin(lat1);
但我几乎可以肯定这是错的。我如何才能得到每个点的位置，并计算出它们在球体直径上的距离，以及它们围绕球体圆周的距离？
谢谢。
浏览 0提问于2009-12-05得票数 1
回答已采纳
1回答
与此非常相关的是，我正在尝试绘制一些世界区域，现在使用墨卡托投影，但在添加x和y限制时总是遇到麻烦：
ggplot(world, mapping = aes(x = long, y = lat, group = group)) +
geom_polygon(fill = "black", colour = "black") +
coord_map(projection = "mercator", xlim = c(-125, -30), ylim = c(-60, 35))
显然不是很好。当我使用coord_cartesian (按
浏览 0提问于2016-10-17得票数 4
1回答
我使用map包将一组经度/经度转换为地图投影(在本例中为lagrange)。接下来，我使用它来定义空间上的其他几个点。我现在想将这些新的点转换回lat/long，但mapproj没有提供这样的功能，我在其他地方也找不到它。
cord<-data.frame(x, y)
lagrange<-mapproject(cord$x, cord$y, projection="lagrange")
#Some
#clustering
#here
#Here are cluster centers, but in the form of lagrange projecti
浏览 1提问于2014-10-08得票数 1
1回答
嗨，我正在尝试绘制华盛顿金县的地图，因为我几乎没有需要放在地图上的数据点。下面是我使用的代码。 long <- c(47.5112,47.7210 ,47.3684,44)
lat <- c(-122.257, -122.319, -122.031,-120)
price <- c(287655,456355,662500,234563)
House <- data.frame(long, lat, price)
states <- map_data("state")
# west_coast <- states %>%
#
浏览 18提问于2019-09-28得票数 3
1回答
因此，我试图将值从具有纬度和经度的数据列中分离出来。
df <- data.frame(localidade = c(drugs_location))
df %>% separate(localidade, c("Lat", "Long"), sep = ",")
我使用上面的代码来分隔Lat和Long中的值，但是现在我需要使用这些值来绘制一个图，使用Lat作为x轴，Long作为y轴。
这是请求的dput(head(df))输出：
structure(list(localidade = c("(42.30676881,
浏览 2提问于2020-09-28得票数 0
回答已采纳
该脚本的目的是使用pylab对Mollweide投影中的一些数据进行等高线绘制。contourf有一个奇怪的行为。以下是代码
import numpy as np
import pylab as plt
ra = np.linspace(-np.pi, np.pi, 40)
dec= np.linspace(-np.pi/2, np.pi/2, 20)
X,Y = np.meshgrid(ra,dec)
Z = np.sin(X) * np.cos(X) * np.sin(Y) * np.cos(Y)
plt.figure()
ax = plt.subplot(111, projectio
浏览 0提问于2013-04-15得票数 2
回答已采纳
我已经实现了一些代码来创建一些代码来处理一个相对较小的位置的图像，比如平面，用于在存储的图像上的位置和传入的Lat/Long信息之间进行转换。
使用提供的公式，我编写了以下代码行
var vector = math.matrix(
[[x1],
[y1],
[x2],
[y2]]);
var matrix = math.matrix(
[[lat1,long1,1,0]
,[-long1,lat1,0,1]
,[lat2,long2,1,0]
,[-long2,lat2,0,1]]);
var solution = math.multiply(math.inv(matrix),vecto
浏览 0提问于2015-07-12得票数 2
回答已采纳
1回答
在自定义楼层图像上映射GPS坐标
、、、、
嗯，我有一张地板的图像。我需要做的是获取用户的当前位置并将其显示在楼层图像上。我知道我需要将经度和经度转换为像素，并将其指向自定义地图。
我关注了这个。但是我不能正确地理解它。一点也不接近。可能是因为它的建筑面积很小(不知道是否可以为小面积开发)。然后我看了看其他的问题。这仅仅是硬编码。我不想这么做。
因此，如果有人给我提供任何逻辑想法或代码样本或教程，这将是有帮助的。什么都行。提前感谢！！:-)
lat = 12.875538;
longitudeVal = 80.227112;
delta_lat = 0.0001; // Adjust Y
delta_long = 0.000
浏览 0提问于2012-08-16得票数 1
回答已采纳
从网上：
由于atan2返回范围-π. +π中的值(即-180°. + 180°)，将结果标准化到指南针方位(在0°. 360°范围内，将-ve值转换为180°. 360°)，转换为度，然后使用(θ+360) % 360，其中%是模块化的。
我正在使用VB6，我的代码如下所示
brng = (θ + 360) Mod 360
其中θ= -68
我得到的答案是292°，但预期答案是248°。
我是不是犯了个错误..？还是我错过了什么？请帮帮我。
更新：
我会进一步解释我的问题，
Dim b As Double
Dim x As Double
Dim y As Double
Dim D
浏览 4提问于2011-11-16得票数 0
1回答
我有一个数据帧，如下所示，我需要在循环中选择一行中的LONG_Y，LAT_Y和下一行中的LONG_X，LAT_X，我得到了这个for循环，但它给了我同一行的Y和X坐标。 df = pd.DataFrame({'ID':[11, 12, 13], 'LAT_X':[-47.401666667,-44.198333333,-44.580753], 'LONG_X':[-22.564722221999997,-19.967777778,-22.500296], 'LAT_Y':[-47.401666667, -45.461666666
浏览 22提问于2019-06-11得票数 0
假设我有一个.shp格式的地图，我用ggplot2绘制它，获得以下结果
然后我在上面画了几个点
ggplot() +
geom_polygon(aes(x=long, y=lat, group=group), data=map,
colour = "grey40", fill = "seagreen") +
coord_map("ortho", orientation=c(55, 0, 0)) +
stat_summary2d(mapping=aes(x=Longitude, y=Latitud
浏览 2提问于2013-08-19得票数 1
回答已采纳
基本上，我有一个POI的列表(名称、lat、long)，我想在UIView上绘制它们，相对于我当前的lat/long。我正在寻找一些将这些POI (lat/long)映射到UIView的最佳实践。
I不想使用MKMapView (不需要显示地图数据)。
我读到：
但我不知道如何从CLLocation到UIView上的(x，y)。我只想画出我现在的位置。因此，例如，如果我的屏幕表示一个20×30公里的区域，我如何将我的POI映射到它们相应的(x，y)坐标？
谢谢。
浏览 2提问于2012-12-15得票数 3
我是R的新手，正在尝试创建一张我的国家的地图，上面覆盖着点(代表城市，点按人口大小)。我已经用geom_sf()创建了地图，并将其存储为地图。
这是我试图复制的代码的一个示例，其中cities是一个包含列state、pop、lat和long (没有缺失值)的df：
map <- map +
geom_point(aes(x = cities$long[1], y = cities$lat[1], size = cities$pop[1])) +
geom_point(aes(x = cities$long[2], y = cities$lat[2], size = cities$
浏览 2提问于2020-05-30得票数 0
1回答
我在R中绘制了一张合唱图，它根据处于燃料贫困家庭的比例，使地方当局黯然失色。
#converts file into a readable format
options(scipen = 999)
#create a data frame from imported CSV file
mydata <- readr::read_csv("/Users/Desktop/R/Map/Fuel_Poverty.csv")
#import shape file
mymap <-st_read("Local_Authority_Districts__Decembe
浏览 11提问于2021-11-12得票数 0
回答已采纳
1回答
我正试图创建一个图例，说明情节是如何编码的。然而，这个传说并没有显现出来。
以下是当前的输出
这与目前的地块已经占据了整个空间有关吗？
下面是我创建情节的代码(并尝试添加一个图例)：
cols <- c("50 miles"="#B10026","100 miles"="#FC4E2A","150 miles"="#FEB24C", "200 miles"="#FFEDA0")
hurricane_plot <- ggplot(data = US_su
浏览 1提问于2019-09-15得票数 2
回答已采纳
2回答
这是我第一次使用javvscript，我正在开发一个基于jason数据的计算器，当我输入经度和纬度时，我得到了与经度和纬度相关联的网格代码--这是json数据的样子
myData =
{
"datagrid": [{"GRIDCODE":1735,"lat":35.91511132,"long":-5.401484264},
{"GRIDCODE":1805,"lat":35.90677815,"long":-5.468149593}
浏览 0提问于2020-04-28得票数 0
回答已采纳
我正在使用matlab中的m_map包来生成一些值的全局映射。
下面是一个示例：
lat2 = [-54.54,-33.15,-8.64,-1.3,5.8,15.35,30.7,32.86,37.54,40.25,...
43.84,45.95,47.65,49.38,51.04,53.33,54.70,56.33,58.59,60.02,60.84,...
62.09,63.15,64.15,66.06];
lon2 = [-68.02,-52.84,26.39,33.22,37.54,-83.84,90.66,119.30,89.34,-121.19,...
-77
浏览 5提问于2015-07-03得票数 2
回答已采纳
我在做一个星系团的散点图。
用sklearn.cluster.KMeans，我得到了和二维矩阵X一样长的列表.
正在运行
k = 5
df = pd.read_csv('data_latlong.csv')
lat = df['Lat'].values
long = df['Long'].values
X = np.matrix(list(zip(lat, long)))
kmeans = KMeans(n_clusters=k).fit(X)
plt.figure(figsize=(10, 10))
plt.scatter(X[:, 0],
浏览 1提问于2019-12-02得票数 0
回答已采纳
1回答
绘制R中的全局投影光栅
、、、
我正试图在等地球投影中绘制一个全局光栅(lon:-180- 180；lat -90- 90) (不重要--可能是Winkel Tripel或Robinson)，但两边的边界重复(见图)。我怎么才能避免这种情况？
和给出的答案只适用于在Mollweide中绘图，而不适用于其他投影。
下面是一个可重复的例子。
library(maptools)
library(raster)
data("wrld_simpl")
r <- raster(ncol=180, nrow=90)
r <- rasterize(wrld_simpl, r, field="UN"
浏览 1提问于2020-09-11得票数 0
回答已采纳
2回答
我一直在广泛地寻找一种从军械测量坐标转换为有效esri坐标的方法。我发现有相当多的页面可以转换为lat long (如果有点关闭)，但没有转换为esri的页面(我相信是utm)。
这是在python或JavaScript / actionscript等中使用的-我不太担心语法，更多的是对所涉及的数学知识的理解。
谢谢
伊恩
浏览 0提问于2011-09-04得票数 0
回答已采纳
1回答
我正在用geom_path绘制海岸线地图上的数据，我无法删除连接第一个和最后一个数据点的线。这方面的数据集相当大，但可以找到。
这个问题被报告并在上解决，尽管它对我的情况没有帮助。
LHplot <- ggplot(data = LH, aes(x = long, y = lat, group=group)) +
geom_path(aes(group=group), size = 1, color = "darkgrey") +
theme_bw() +
theme(axis.line.y=element_blank(),
axis.line.x
浏览 0提问于2018-11-22得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
