Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  605 1282 (N27E20 L=336)   213     73        0 ( 580/1261- 640/1312) -51.3 E21-E17 N24-N29
      2 1805 1222 (S04W56 L= 52)   384     48        1 (1779/1197-1830/1244)   0.0 W53-W59 S05-S03
      3 1530 1106 (S04W31 L= 27)  1662    449        1 (1414/ 911-1621/1239) -76.6 W26-W34 S16-N05
      4 1722 1187 (S04W47 L= 43)   514     98        1 (1670/1164-1762/1216) -23.2 W43-W51 S06-S01
      5  578 1175 (N21E23 L=332)   523    177        1 ( 482/1153- 648/1200)  -9.5 E30-E18 N18-N22
      6 1824 1109 (S11W55 L= 50)  1160    198        3 (1794/1023-1861/1196)  56.3 W50-W59 S15-S06
      7  926 1131 (N11E03 L=352)   348    114        0 ( 868/1110- 989/1156) -38.7 E06-W00 N08-N13
      8 1718 1060 (S11W44 L= 40)   448     97        1 (1695/1027-1747/1115)  53.1 W41-W46 S12-S08
      9 1870  857 (S27W58 L= 53)   880    130        1 (1826/ 829-1923/ 904) -23.2 W52-W66 S29-S23
     10 1121  799 (S12W00 L=356)   656    163        5 (1053/ 760-1185/ 832)   0.0 E04-W04 S13-S09
     11 1165  654 (S21E00 L=355)   756    304        0 (1072/ 592-1220/ 797)  36.9 E07-W05 S23-S13
     12 1709  674 (S34W40 L= 35)   478     50        1 (1685/ 646-1737/ 708) -50.2 W37-W42 S37-S31
     13  681  656 (S10E29 L=327)   408    147        2 ( 608/ 622- 735/ 699)   0.0 E34-E24 S11-S08
     14 1752  576 (S43W47 L= 42)  3209    235        9 (1689/ 468-1814/ 650)  61.4 W38-W58 S50-S37
     15  781  450 (S25E30 L=326)  2049    550        6 ( 626/ 255- 961/ 604)  22.3 E52-E13 S35-S19
     16  426  427 (S19E58 L=298)   393     57        0 ( 403/ 396- 439/ 452)  45.0 E60-E55 S20-S17
  

    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