Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  380 1120 (S11E43 L=296)   545    112        1 ( 362/1059- 404/1178) -66.8 E44-E41 S14-S07
      2 1602 1086 (N08W34 L= 13)   369    109        0 (1585/1028-1613/1151) -82.9 W33-W35 N04-N12
      3  293 1066 (S15E50 L=289)  1208    343        4 ( 184/1025- 370/1118)   0.0 E60-E42 S17-S11
      4  120  920 (S25E67 L=272)  1634    241        0 (  94/ 837- 148/ 976)  69.4 E69-E63 S29-S21
      5 1771  933 (N04W51 L= 29)   641    101        3 (1748/ 883-1792/ 975) -51.3 W48-W53 N01-N06
      6 1302  728 (S18W23 L=  1)   420    146        0 (1220/ 707-1359/ 768)  32.0 W17-W25 S20-S14
      7 1616  719 (S12W44 L= 23)   538    152        0 (1559/ 688-1671/ 740)   0.0 W39-W49 S13-S10
      8  432  628 (S41E37 L=302)   796    322        1 ( 346/ 564- 516/ 721)  37.4 E46-E29 S46-S33
      9 1693  592 (S17W56 L= 34)   595    157        1 (1653/ 551-1734/ 635) -33.7 W50-W58 S19-S14
     10 1482  470 (S29W44 L= 23)  3869    747       11 (1227/ 435-1724/ 529)  16.7 W26-W66 S35-S20
  

    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