Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1408 1583 (N29W43 L=156)   313     35        2 (1388/1561-1430/1603)  21.8 W40-W45 N28-N30
      2  288 1089 (N25E51 L= 62)   797    124        0 ( 259/1030- 321/1138) -81.9 E54-E47 N21-N28
      3  357 1065 (N23E44 L= 69)   396     59        1 ( 337/1030- 379/1100)  53.1 E46-E42 N20-N24
      4  388  969 (N16E43 L= 70)   472    158        1 ( 356/ 902- 442/ 998)  33.7 E45-E38 N12-N18
      5 1595  927 (S13W33 L=146)   451    109        0 (1537/ 911-1645/ 942) -14.0 W28-W36 S14-S10
      6 1772  897 (S20W47 L=160)   362     74        0 (1756/ 857-1787/ 932)  36.9 W45-W49 S21-S17
      7 1061  562 (S21E09 L=104)  1937    471        0 ( 900/ 415-1182/ 725)  29.2 E23-W02 S27-S13
      8 1307  574 (S27W07 L=120)   510    193        0 (1259/ 507-1332/ 662)  45.0 W02-W10 S30-S21
      9  770  415 (S24E32 L= 80)  1801    377        1 ( 715/ 338- 878/ 556) -86.4 E37-E25 S31-S15
     10  690  280 (S33E48 L= 65)   638    107        0 ( 662/ 246- 715/ 322)  50.2 E53-E44 S35-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