Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  782 1802 (N45E33 L=261)  1098    123        3 ( 728/1776- 852/1825)  -4.1 E38-E24 N42-N47
      2 1118 1703 (N40E01 L=292)   306     77        0 (1078/1677-1157/1721)  45.0 E03-W01 N37-N42
      3  441 1483 (N19E46 L=248)  1503    267        1 ( 377/1398- 534/1584) -29.7 E56-E35 N13-N25
      4 1695 1292 (N19W43 L=336)   788    174        1 (1629/1243-1778/1345) -13.0 W37-W51 N17-N21
      5 1390 1190 (N08W21 L=314)   230     38        0 (1361/1176-1420/1208)   0.0 W18-W22 N07-N09
      6 1884  886 (S02W66 L=359)  1145    132        2 (1834/ 843-1941/ 917)  18.4 W60-W75 S05-N00
      7 1119  837 (S16W08 L=301)   212     60        0 (1087/ 819-1153/ 853)   0.0 W05-W09 S16-S14
      8  634  701 (S28E23 L=270)  1185    368        0 ( 559/ 555- 740/ 839)  55.0 E30-E16 S38-S18
      9 1776  601 (S21W65 L=358)   899    119        2 (1725/ 565-1823/ 623)  15.3 W58-W73 S22-S18
  

    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