Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  826 1801 (N40E36 L=328)   790     95        1 ( 781/1773- 881/1826)  -4.8 E41-E29 N37-N41
      2 1169 1706 (N40E05 L=358)  1121    123        3 (1116/1652-1209/1758) -41.2 E09-E01 N35-N43
      3 1066 1723 (N39E13 L=350)   617    118        1 (1006/1680-1111/1754) -12.5 E18-E09 N35-N40
      4  334 1598 (N18E67 L=297)   643     76        2 ( 311/1555- 361/1634) -63.4 E69-E64 N14-N20
      5  498 1422 (N08E41 L=322)   328     74        0 ( 467/1391- 529/1455)  53.1 E42-E39 N05-N10
      6 1115 1367 (N15E01 L=  2)   254     89        0 (1075/1333-1144/1384) -18.4 E03-W00 N13-N15
      7 1527 1344 (N23W26 L= 29)   250     60        0 (1501/1319-1554/1371) -26.6 W23-W27 N21-N23
      8  935 1345 (N10E11 L=352)   199     70        0 ( 904/1327- 968/1365)  14.0 E13-E09 N09-N11
      9  195 1149 (S11E62 L=301)  1516    195        6 ( 137/1073- 256/1220)  36.9 E69-E56 S15-S06
     10  303  990 (S20E49 L=315)   417    146        1 ( 289/ 928- 312/1076)  90.0 E50-E47 S23-S14
     11  265  778 (S33E54 L=309)   754    145        2 ( 220/ 733- 297/ 847)  38.7 E61-E49 S36-S28
     12 1607  735 (S11W43 L= 47)   618    200        0 (1587/ 640-1630/ 822) -65.6 W41-W47 S16-S05
     13 1293  650 (S23W25 L= 28)   958    268        1 (1263/ 524-1333/ 786) -61.9 W20-W28 S30-S15
     14 1682  690 (S11W51 L= 54)   458     96        2 (1657/ 648-1703/ 740) -38.7 W47-W54 S13-S08
     15 1542  493 (S26W50 L= 53)  1906    365        3 (1418/ 401-1671/ 589)  59.5 W43-W57 S34-S17
     16 1735  526 (S19W66 L= 70)   443     71        0 (1691/ 511-1768/ 540)  20.6 W62-W70 S20-S17
     17 1586  394 (S30W62 L= 65)  1101    157        2 (1517/ 355-1651/ 438)  56.3 W56-W65 S34-S25
  

    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