Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1408 1587 (N30W43 L=156)   328     38        1 (1382/1562-1430/1608)  11.3 W41-W46 N28-N30
      2  292 1089 (N25E51 L= 62)   797    115        3 ( 260/1029- 326/1138) -74.1 E54-E47 N21-N28
      3  363 1069 (N23E44 L= 69)   379     70        1 ( 343/1035- 385/1102)  45.0 E45-E41 N20-N24
      4  396  972 (N16E42 L= 71)   496    150        2 ( 357/ 903- 454/1000)  39.8 E45-E37 N12-N18
      5 1601  930 (S12W33 L=146)   386    104        1 (1545/ 913-1646/ 945) -15.9 W29-W36 S14-S10
      6 1775  895 (S20W48 L=160)   339     79        0 (1761/ 853-1790/ 931)  53.1 W45-W49 S22-S18
      7 1064  564 (S21E09 L=104)  1790    550        2 ( 900/ 419-1186/ 726)  27.5 E23-W02 S27-S13
      8 1327  623 (S24W09 L=122)   231     77        1 (1316/ 580-1339/ 663)  59.0 W07-W10 S26-S21
      9  774  414 (S24E32 L= 81)  1661    369        1 ( 722/ 338- 881/ 553) -82.4 E36-E25 S31-S15
     10  693  280 (S33E47 L= 65)   617    111        0 ( 665/ 246- 717/ 324)  45.0 E52-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