Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1072 1750 (N40E23 L= 61)   322     40        0 (1050/1733-1098/1771) -11.3 E25-E20 N38-N41
      2  505 1744 (N23E65 L= 19)   649     57        1 ( 475/1720- 530/1764)  -5.2 E72-E60 N21-N24
      3 1296 1498 (N29W02 L= 86)  1135    262        1 (1205/1394-1346/1581) -33.7 E05-W07 N23-N33
      4  963 1424 (N15E15 L= 69)   278    114        0 ( 942/1371- 988/1474) -74.1 E15-E13 N11-N18
      5 1079 1328 (N12E05 L= 79)   384    126        1 (1037/1274-1107/1366) -45.0 E08-E02 N09-N14
      6 1629 1232 (N24W31 L=115)   442     65        0 (1592/1201-1651/1258)   0.0 W28-W33 N21-N25
      7 1315 1183 (N11W12 L= 96)   810    311        1 (1197/1067-1368/1250) -15.9 W03-W17 N06-N13
      8 1604 1081 (N14W33 L=117)   294     59        0 (1588/1047-1620/1116)  90.0 W31-W33 N12-N16
      9 1702  980 (N12W43 L=127)   622      0       -1 (1682/ 895-1747/1054) -60.9 W40-W46 N07-N16
     10 1802  990 (N16W52 L=136)  1036    122        0 (1748/ 927-1855/1045)  -8.1 W45-W59 N14-N17
     11 1687  469 (S15W67 L=151)  1520    252        5 (1636/ 386-1728/ 558) -57.5 W60-W75 S20-S09
     12 1219  512 (S28W29 L=113)   433    103        0 (1162/ 497-1265/ 528)  33.7 W25-W31 S30-S25
     13  634  478 (S50E10 L= 74)   499     76        0 ( 607/ 454- 667/ 516) -18.4 E12-E06 S52-S47
     14  765  489 (S45W01 L= 85)   264     40        1 ( 746/ 467- 795/ 507)  63.4 W00-W02 S46-S42
     15 1202  362 (S38W37 L=122)   514    132        0 (1150/ 312-1256/ 403)  83.7 W36-W38 S42-S33
     16 1379  361 (S31W49 L=134)   585    130        0 (1322/ 338-1449/ 385)  49.4 W46-W53 S34-S27
     17 1447  309 (S32W60 L=144)   763    129        2 (1419/ 259-1475/ 361) -32.5 W55-W67 S34-S27
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y