Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  872 1701 (N35E22 L=196)   665    125        3 ( 834/1659- 920/1730)  51.3 E24-E18 N31-N37
      2 1180 1393 (N18W05 L=223)   236     39        1 (1154/1372-1201/1417)  71.6 W03-W06 N16-N19
      3  792 1144 (S02E15 L=203)   386    115        0 ( 774/1086- 816/1205) -59.0 E16-E13 S05-N01
      4 1918 1141 (N16W64 L=283)  1564    255        1 (1870/1096-1952/1186) -38.7 W58-W69 N14-N18
      5  594  974 (S14E26 L=192)  1141    329        0 ( 477/ 864- 726/1088)  45.0 E33-E18 S21-S05
      6  385  982 (S15E41 L=177)   780    240        3 ( 291/ 928- 477/1030)  26.6 E48-E34 S19-S11
      7 1228  808 (S16W15 L=234)   473    169        0 (1161/ 751-1269/ 855) -35.5 W10-W18 S18-S13
      8 1262  727 (S20W19 L=237)   249     90        1 (1208/ 711-1306/ 747)  18.4 W15-W21 S21-S18
      9  518  648 (S35E31 L=187)   396     57        0 ( 490/ 624- 545/ 677) -36.9 E33-E28 S36-S33
     10  862  214 (S65W03 L=221)   527     57        1 ( 834/ 195- 891/ 240)  -7.1 E01-W07 S67-S62
  

    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