Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  410 1484 (N04E53 L= 37)   332     83        0 ( 399/1440- 422/1519) -59.0 E55-E50 N01-N06
      2 1796 1432 (N45W49 L=139)  1127    121        1 (1741/1397-1850/1456)  11.3 W40-W61 N41-N47
      3  579 1211 (S08E31 L= 59)   345     86        0 ( 552/1177- 610/1251)  78.7 E32-E29 S09-S04
      4 1752 1070 (N19W44 L=134)   593    159        1 (1712/1001-1804/1122) -15.3 W40-W51 N16-N21
      5  605  993 (S19E24 L= 66)   486    143        0 ( 545/ 945- 657/1030)  58.0 E26-E20 S23-S15
      6  429  986 (S24E37 L= 53)   573    202        1 ( 341/ 951- 509/1009)  12.1 E43-E29 S27-S21
      7  322  779 (S39E45 L= 45)   925    130        1 ( 277/ 731- 373/ 846)  39.3 E50-E39 S43-S33
      8 1314  518 (S25W35 L=125)  1367    234        5 (1240/ 451-1420/ 562)  63.4 W28-W40 S30-S20
      9 1175  406 (S36W33 L=123)   757    163        0 (1135/ 349-1218/ 490) -28.3 W26-W39 S38-S31
     10 1027  328 (S47W29 L=119)  1280    199        1 ( 973/ 277-1104/ 415) -16.4 W21-W38 S49-S41
  

    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