这个问题有点陈旧,但我的经验也是如此:差不多 30 年前我做过类似的事情。
这会在位图上产生简单的轮廓:
- 计算网格中每个点的场强(我假设您正在尝试根据点的质量绘制诸如引力场等高线之类的东西)。
- 用两种替代颜色为等高线之间的替代空间(您还没有得到)着色。例如如果等高线应该相隔 100 个单位(场强),则根据
ToInt(pixel_field_strength / 100) % 2 选择每个像素的颜色。
- 跟踪颜色边界的边缘以生成轮廓。例如,如果您的两种颜色是白色和黑色,则只保留与黑色像素相邻的白色像素。
如果您只对结果感兴趣,请使用 cmets 中建议的库。
纯粹出于怀旧之情,我找到了我原来的 BBC BASIC 代码。它仍然在 modelb(一个 BBC Micro 模拟器)上运行。
10 REM THIS COMES WITH NO WARRANTY!
20
30 REM Gravity field
40
50 MODE 1
60 PROCsetup
70 FOR Y%=300 TO 900 STEP 4
80 FOR X%=200 TO 800 STEP 4
90 R=LOG(FNforce(X%,Y%))
100 GCOL0,((R*10) MOD 2)+1
110 PLOT69,X%,Y%
120 NEXT
130 NEXT
140 PROCcontour
150 VDU19,1,0,0,0,0
160 VDU19,2,0,0,0,0
170 END
180 DEFPROCsetup
190 N%=5
200 DIM X%(N%),Y%(N%),M%(N%)
210 FOR P%=1 TO N%
220 READ X%(P%),Y%(P%),M%(P%)
230 NEXT
240 ENDPROC
250 DATA 625,625,1000000
260 DATA 425,725,1000000
270 DATA 475,425,1000000
280 DATA 375,575,1000000
290 DATA 725,525,1000000
300 DEFFNforce(X,Y)
310 P=0
320 FOR P%=1 TO N%
330 DX=X%(P%)-X:DY=Y%(P%)-Y
340 R=SQR(DX*DX+DY*DY)
350 P=P+M%(P%)/R
360 NEXT
370 =P
380 DEFPROCcontour
390 GCOL0,3
400 FOR Y%=300 TO 900 STEP 4
410 FOR X%=200 TO 800 STEP 4
420 IF POINT(X%,Y%)=1 AND (POINT(X%+4,Y%)=2 OR POINT(X%-4,Y%)=2 OR POINT(X%,Y%+4)=2 OR POINT(X%,Y%-4)=2) THEN PLOT69,X%,Y%
430 NEXT
440 NEXT
450 ENDPROC