home / github / issue_comments

Menu
  • Search all tables
  • GraphQL API

issue_comments: 863117814

This data as json

html_url issue_url id node_id user created_at updated_at author_association body reactions performed_via_github_app issue
https://github.com/pydata/xarray/issues/2171#issuecomment-863117814 https://api.github.com/repos/pydata/xarray/issues/2171 863117814 MDEyOklzc3VlQ29tbWVudDg2MzExNzgxNA== 80449808 2021-06-17T10:17:53Z 2022-04-19T03:15:18Z NONE

i am not yet experimenting enough on python. My problem is that I wrote some code to draw a line using a colon at the beginning and at the end of my line I used numpy exclusively and since it is an xrarray array it takes me a long time. My data is all of dimensions lat, lon, SST, (1501,2001), and here is my code:

import xarray as xr
import glob
from numpy import *
import numpy as np
import matplotlib.pyplot as plt
from scipy.signal import get_window
import math
import pandas as pd
import numpy as np

filepaths_TEMP = "./ 2014 / MUR-SST_2014__Day_027_27-Jan-2014.nc"

data_sst = xr.open_dataset (filepaths_TEMP)
lon_sst = data_sst.variables ['lon'] [:]
lat_sst = data_sst.variables ['lat'] [:] #
sst = data_sst.variables ['SST'] [:]
#% to draw the line
latlim = [50, 50.5];
lonlim = [-10.2];
n1 = 1200; #% number of point along profile
delta_lon = -10 - (2);
delta_lat = 50 - 50.5;

#cd (input_path1)
#matfiles = dir (fullfile ('*. nc'));
# %%
Flat = latlim [1] + np.arange (0, n1) * delta_lat / n1
Plon = lonlim [1] + np.arange (0, n1) * delta_lon / n1
transect = np.zeros ((1200,1))
for ii in range (0, len (Flat)):

    lat_S1 = Flat [ii];
    lon_S1 = Plon [ii];
    a_1 = abs (lat_S1-lat_sst);
    a_2 = abs (lon_S1-lon_sst);
    a_3 = (a_1 + a_2);
    m_3 = a_3 [:]. min ()
    [px, py] = np.where (m_3 == a_3);
    px = px [0];
    py = py [0]; # clear a_1 a_2 a_3 m_3
    transect [ii] = sst [px, py];

My question is, isn't there a way to make it simpler and faster, especially thank you very much in advance.

{
    "total_count": 0,
    "+1": 0,
    "-1": 0,
    "laugh": 0,
    "hooray": 0,
    "confused": 0,
    "heart": 0,
    "rocket": 0,
    "eyes": 0
}
  325439138
Powered by Datasette · Queries took 0.668ms · About: xarray-datasette