[Thread Prev][Thread Next][Index]
[ferret_users] @FNR in 2D fills up x dimension first
- To: ferret users <ferret_users@xxxxxxxx>
- Subject: [ferret_users] @FNR in 2D fills up x dimension first
- From: Hella Riede <Hella.Riede@xxxxxx>
- Date: Thu, 28 Jun 2018 16:37:55 +0200
- Arc-authentication-results: i=3; mx.google.com; arc=pass (i=2); spf=softfail (google.com: domain of transitioning hella.riede@xxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=hella.riede@xxxxxx
- Arc-authentication-results: i=2; mx.google.com; arc=pass (i=1 spf=pass spfdomain=dwd.de); spf=softfail (google.com: domain of transitioning hella.riede@xxxxxx does not designate 65.127.216.228 as permitted sender) smtp.mailfrom=Hella.Riede@xxxxxx
- Arc-authentication-results: i=1; mx.google.com; spf=pass (google.com: domain of hella.riede@xxxxxx designates 141.38.3.245 as permitted sender) smtp.mailfrom=Hella.Riede@xxxxxx
- Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=message-id:content-language:mime-version:user-agent:date:subject :organization:from:to:arc-authentication-results :arc-message-signature:arc-authentication-results :arc-message-signature:arc-authentication-results; bh=QGMEJE2E570uVXIm+6CCwai4RuzfvY6xMTwrr9vxXdE=; b=VrGsGK5cEUFfsdQK/CkKEcRbO4GokMDKwYK7uE5+0i/Hd52BVGFx0PCu655HpBTN5F ahgMqDE5DTOBfR/+J34R9kEv+GHrI8k8iTHpHUfo+jGY6qMLUJ4d44Po3DGW10yMnewt d1lVpIoDu7/qhjHV6vrBcCpRTZZI/imZ7+iO+c6qG7UbMnqea2zWgFhU/9at9U8Y4Uvr zj22F+t68Rf2fM/GVXwBGDy0HZ3Q9qywLlKdQdFmjytL0YjBLsUvnsnSDUyOv2KwzRCf ROalDSGItAWSg8nLMt51+9+PUYvzqz+pOuWwbinv3NcC+iZqAm6gptPA5f6wJLhGTAXG E96A==
- Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=content-language:mime-version:user-agent:date:message-id:subject :organization:from:to:arc-authentication-results :arc-message-signature:arc-authentication-results; bh=QGMEJE2E570uVXIm+6CCwai4RuzfvY6xMTwrr9vxXdE=; b=SKa4nwEJUP4AXDT49EWtpIG7JufzNNzd8hHIuXTLai/MJi3LhLz22zlbRB66gGPfI0 M7APu/RwuDUuJ9jDeSRffv42DAM5OXzvAjsCc0CkFDsnwAyTPjJKa110x0lEiQGjp9VS XDk/Em4VJAZGPTHbbP0d9SpRPC67AthJe6TDyIlgk9RHVQb2mHUC2p6z322/RyE3jske Su+OO84cmzC3T0A+2dbH3Yia5Lh8UiMU4oHlzOgkaEvINOGMeLa66MNW9qhuxdi/77Pk N2DnXh7xMxOTd3TbO8hS7IQ5D+/zNBvz8sEOFvhFcjKrrNGAijk7qQorJcMIPWWGsgCM hudw==
- Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=content-language:mime-version:user-agent:date:message-id:subject :organization:from:to:arc-authentication-results; bh=QGMEJE2E570uVXIm+6CCwai4RuzfvY6xMTwrr9vxXdE=; b=o9ZBTVV1yxRdvCSHLCnsA9ygQ655RByz4R0bWxpOgzGWng69GfpIlFPGIH4Z1VJFuF BgVAYDpRguuYBBbvSdCPIwnnk6UoASqUcZj0tcDxZINfhhAcrLVucUCeM/21Vyi0NdOu mmtbx49oChXkflZOzC+uVDqWa6ahNINMdByohsmZu2UJqi1GDPUzWOTe95Zhe23Vr1Iu hDYrPg8ewvcwqEZ4JgTunuTCKUfWQJW0WKXXiMDaa5aCVG0flpbPAaZQ1ncBiDe2Jewo vIkHg3PAX52LiSa+2hjSJTelzNOm9dP6jR06HRZDiUqZIdb1uSn/2kKkMvKwzEdtec6+ bEzg==
- Arc-seal: i=3; a=rsa-sha256; t=1530196688; cv=pass; d=google.com; s=arc-20160816; b=n6l0Ie4M9SocdcJHBdzLa6AtqffSOvHdpkQMGfaepbBLWV0kK3fNFGKCI37kvusgM/ r4Z6yfoGAioj3W1QFhMrbMVTjT0UnFS9kM6X5Js1Cvm2rP0W+mAfHjm4LHkMp437os4s KTs8ybsmSe3DWPku+g33qyN3ROcb2V4Dlvn0+iQabdZ57rVyV71ZdKluYOQt6IN3GiUx LtbKVUC2NXIc74qWCsu4Xmz11HnyiJUxUX3C+VrlLjzYNV/t+Q3b/mD7MxtmZGnPNKVr Rz3XbnDStyR7rVf6eGEx2+JeSO0N7Q3h8Rh0fxcuJegtHesuJlvRV8qVy1OMwZnojkC9 rD/g==
- Arc-seal: i=2; a=rsa-sha256; t=1530196682; cv=pass; d=google.com; s=arc-20160816; b=Lm81U26sJCYqf0IC6a1QphSnkIDp4taIVcnw5JhIInUpqgaKMj8kY5KuuxCfekGLkp NZ6e5jfMJkr72Xb8EzzvwERpGz5aeQeWhMNmlQPQQo0N4ixASHVwBgCH0COTeuRNlPT0 LWyRHWPDBDZjMICvNz1e5bnv8WKZacV3pJ1CdlM4teK3/VDeYXtSz6KSVV4dVWeYNkHu ovfWreQpdFqfSwC5APkVWU+kGGhhYLTku7CjcRNZ3G7G23JqBqHPUTiNyvqCeI31/z+m RLunZWUco8lgUdzd0cBGE0qpW/8GCDDMiEFp04iMoAwtE3wpJb/segXSZmyN5mDdIE/l qwhA==
- Arc-seal: i=1; a=rsa-sha256; t=1530196676; cv=none; d=google.com; s=arc-20160816; b=MWAZoEzam9WAPVgRAPy6HCnGqI4xRcMwaYz5ecyYzQtgI6anpCDeIkhPlUkPzItYq/ VkWN8/iN2/POj2r/MT6cVQ0/tuDISntVBG0vBDDoBM42siIJ3nW09x2c++LjSHcKJFNg 4mQeNaBlvXx/SUUgP8KBJA0TYFxWl/Wc3Ho/3P0mgnqWp+k7vG1QvcUEbMYE8hd2xRaT SKZOZI2fjuiVq04Xg61Kfc28MNc2f7wc37H7mMYdOdGFg7jKBxyh9IXfDRgXyvjO9j1t j/9rCIZN1pwoke4odnQKfecJwT9Jj1cfOrI+uJabPCZxieHBFNmXcx3Jc3DlDnHefEIq DlHw==
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=sender:to:from:organization:subject:date:user-agent:mime-version :content-language:message-id:x-original-sender :x-original-authentication-results:precedence:mailing-list:list-id :list-post:list-help:list-archive; bh=QGMEJE2E570uVXIm+6CCwai4RuzfvY6xMTwrr9vxXdE=; b=fATgR5QLXhTzpYPKuMOpNYl0V/ZiCacuUsHTjjL8zsYmgRKo60KonDZpat52Zirg9o c2dY9hAqzPlvwd0UbR6hlhIpMqgcCKY+IsTJjZ+IxV1xeQo37Y8rnf/9Jf5Ndz+7NYKZ DURYYRJso9RxQtHYcXYRtxGrxFGTvXk+VWTzyYhwG1dTNhEa7v9Ba0daecKkYtgW4Z8+ QT3hs74GsK15aJFiqTmgXzaMMKVqI8+FJTOgTa/O5ft8W5KikV7tCsvcUWWbhfIhUqTp AZ6C6xNxmEr16CR8FCfqJ02AO1ZkV0TLBRXMMLbkY7pmsROvWkoSpp11zgYNts05vrAN KbgQ==
- List-archive: <https://groups.google.com/a/noaa.gov/group/ferret_users/>
- List-help: <https://support.google.com/a/noaa.gov/bin/topic.py?topic=25838>, <mailto:ferret_users+help@noaa.gov>
- List-id: <ferret_users.noaa.gov>
- List-post: <https://groups.google.com/a/noaa.gov/group/ferret_users/post>, <mailto:ferret_users@noaa.gov>
- Mailing-list: list ferret_users@xxxxxxxx; contact ferret_users+owners@xxxxxxxx
- Organization: Deutscher Wetterdienst
- Sender: owner-ferret_users@xxxxxxxx
- User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:52.0) Gecko/20100101 Thunderbird/52.8.0
Hello ferret community,
I am using @FNR in two dimensions and got an unexpected result (see
ferret script at end of this mail):
In the attached ferret_before_FNR.png, the original 2D grid cells are
shown (values inherited from dots by scat2grid, see below).
ferret_after_FNR.png shows the result after applying @FNR (fill with
nearest neighbor) in 2 dimensions.
As an example for the unexpected result, lets take the grid cell at
(x=8,y=5). Its value before @FNR was undefined. The nearest neighbor is
(x=8,y=8). Yet, the value is inherited from (x=2,y=4).
Is @FNR not supposed to be applied in 2D or is the logic behind
different from what I expected? It seems the missing points are always
filled up in x dimension first, and only remaining missing values are
then filled up in y direction ...
I attached a manual version of what I would have expected as
ferret_FNR_expected.png.
Thank you very much in advance for helping to clarify this.
Best regards,
Hella Riede
can var/all
can sym/all
let a = {1,1,8,8,2}
let b = {1,8,8,1,4}
let val = {2,4,6,8,10}
def ax/x=1:8:1 xax
def ax/y=1:8:1 yax
def sym lev = (1.5,10.5,1)
!! assign values from dots to grid cells
set win 1
let gridded = SCAT2GRIDGAUSS_XY(a,b,val,x[gx=xax],y[gy=yax],0.8,0.8,0.8,0)
shade/lev=($lev) gridded
!! plot original scattered values for comparison
go polymark poly/lev=($lev)/key/ov/line/fill a b val circle 2
!! fill up undefined values with nearest neighbor
set win 2
let nn = gridded[x=@FNR,y=@FNR]
shade/lev=($lev) nn
go polymark poly/lev=($lev)/key/ov/line/fill a b val circle 2
Attachment:
ferret_before_FNR.png
Description: PNG image
Attachment:
ferret_after_FNR.png
Description: PNG image
Attachment:
ferret_FNR_expected.png
Description: PNG image
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
Ferret
Privacy Policy | Disclaimer | Accessibility Statement