Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1580 1522 (N27W50 L= 84)   568     79        0 (1540/1478-1609/1563) -80.5 W48-W51 N23-N29
      2  966 1412 (N31W02 L= 37)   302    112        0 ( 911/1398-1023/1435)   8.1 E01-W06 N30-N31
      3 1802 1383 (N13W65 L=100)   590    159        0 (1769/1362-1835/1414)   0.0 W61-W71 N12-N14
      4  551 1313 (N30E29 L=  6)  1718    374        2 ( 398/1266- 710/1369)   2.3 E41-E16 N27-N32
      5 1128 1179 (N14W09 L= 43)   687    169        5 (1065/1145-1204/1220) -40.6 W05-W12 N10-N16
      6  521  711 (S07E37 L=358)  1906    350        5 ( 474/ 567- 564/ 869)  68.7 E44-E33 S15-N03
      7  772  758 (S07E19 L= 15)   622    148        1 ( 689/ 723- 843/ 784) -24.0 E24-E15 S09-S04
      8 1614  700 (S24W35 L= 70)   836    149        1 (1550/ 649-1670/ 738)  18.4 W30-W40 S26-S21
      9  702  461 (S25E31 L=  3)   489    171        1 ( 656/ 402- 744/ 510)  35.0 E36-E26 S29-S21
     10  555  305 (S36E51 L=343)  1038    207        1 ( 506/ 252- 594/ 375)  71.6 E55-E49 S41-S29
  

    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