Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  829 1707 (N35E25 L=196)   644    106        3 ( 787/1667- 874/1739)  51.3 E28-E22 N31-N37
      2  920 1585 (N27E15 L=206)   378    107        0 ( 889/1532- 960/1621)  74.1 E16-E13 N23-N30
      3  735 1207 (N01E20 L=202)   709    283        0 ( 692/1068- 776/1331) -60.3 E24-E15 S06-N08
      4 1588 1220 (N14W33 L=254)   278     66        2 (1563/1198-1617/1247) -33.7 W30-W34 N12-N15
      5 1887 1089 (N12W61 L=282)  2515    300        2 (1839/ 947-1924/1209) -72.6 W54-W66 N03-N19
      6  543  977 (S14E29 L=192)   915    318        1 ( 433/ 868- 663/1088)  47.0 E36-E22 S21-S06
      7  356  992 (S15E44 L=178)   742    284        0 ( 257/ 935- 455/1042)  20.6 E52-E36 S18-S11
      8 1194  822 (S16W13 L=234)   337    118        2 (1141/ 767-1219/ 864) -39.8 W09-W15 S18-S13
      9  481  657 (S35E34 L=187)   401     60        1 ( 451/ 634- 511/ 692) -31.0 E36-E31 S35-S32
     10  862  215 (S65W02 L=224)   541     67        1 ( 834/ 194- 894/ 243) -29.7 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