我试图对一组数据进行FFT处理(它们是10幅图片,中间画的是“T”),并计算给定区域的光谱幅度。
这是我的代码:
for i=1:10
filename = sprintf('T%d.GIF', i);
%fourier transform and showing
f = imread(filename); %read in image
z = fft2(double(f)); % do fourier transform
q = fftshift(z); % puts u=0,v=0 in
我有一幅有270个光谱波段的农业领域的高光谱光栅图像。我创建了一个多边形shapefile,它描述了哪些像素属于每个处理。有250个单独的多边形,每个多边形对应于每个处理的重复。我想找出每个多边形内所有像素的每个波段的平均像素值。 Image of raw hyperspectral data Image of polygons delineating treatments 我尝试在ArcGIS和QGIS中使用区域统计工具,但这两个工具一次只能在一个波段上运行统计。这样做270次似乎有点过分。 我还尝试使用ArcGIS中的分割栅格工具将栅格划分为对应于每个多边形的250个单独的栅格。分割栅格
我正在用R中的pavo包分析花色光谱,并且需要平均来自同一物种的所有光谱。光谱都在一个文件夹中,所以我需要分别阅读它们,并将它们平均起来,这样我就得到了每种植物的平均光谱。我对每个物种都有不同的光谱。
我如何指定要根据每个光谱的前四个字母来聚合光谱呢?在下面的例子中,我想把所有的多聚半乳糖光谱和所有的大分子光谱综合起来。在这个论坛上,我没有看到任何关于楼阁的帖子。提前感谢您的帮助。
我尝试使用函数aggspec,我没有设法根据每个频谱文件名的前几个字母进行聚合。
到目前为止,我只成功地读取了所有光谱,绘制了它们,并分别获得了每个文件所需的变量(色调、气密性、色度),但我需要按物种平均光谱。
l
我有一个RGB三胞胎的列表,我想用这样的方式来绘制它们,使它们形成一个类似光谱的东西。
我把它们换成了HSV,这似乎是人们推荐的。
from PIL import Image, ImageDraw
import colorsys
def make_rainbow_rgb(colors, width, height):
"""colors is an array of RGB tuples, with values between 0 and 255"""
img = Image.new("RGBA", (wid
我是R的新手,我开始在那里使用igraph进行图形可视化。下面的例子创建了一个由10个顶点组成的简单网络,并根据颜色值对它们进行着色(在本例中,为简单起见,我将其设置为与顶点的ids相同)。
library(igraph)
vertices <- 1:10
first <- 1:10
second <- c(2:10,1)
edges = cbind(first,second)
color = 1:10
net = graph_from_data_frame(edges,vertices=vertices ,directed=F )
V(net)$color = color