Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1658 1556 (N16W61 L=140)   417     65        1 (1626/1537-1683/1574)  -8.1 W57-W64 N14-N16
      2 1594 1455 (N13W49 L=128)  1291    190        2 (1531/1374-1659/1517) -79.7 W45-W51 N06-N17
      3  241 1460 (N52E66 L= 13)  1104     96        2 ( 192/1440- 288/1476)   0.0 E85-E55 N50-N53
      4  877 1377 (N29W01 L= 80)  1163    289        2 ( 753/1332-1014/1423) -35.5 E06-W09 N24-N34
      5 1244 1330 (N16W21 L=100)   679    363        0 (1092/1311-1410/1360) -20.2 W12-W31 N12-N20
      6  253  936 (N18E55 L= 24)  1869    387        3 ( 225/ 792- 301/1121)  75.3 E58-E49 N09-N28
      7 1880  998 (S21W60 L=139)   871     93        2 (1853/ 961-1906/1049) -90.0 W56-W64 S24-S18
      8 1374 1007 (S05W19 L= 98)   213     59        0 (1359/ 977-1388/1041)  56.3 W18-W20 S06-S03
      9 1681  855 (S23W37 L=116)   538    115        0 (1648/ 806-1714/ 904)  33.7 W33-W40 S24-S20
     10  527  321 (S24E60 L= 19)  1205    140        2 ( 496/ 268- 566/ 397)  26.6 E71-E53 S27-S18
     11  777  282 (S33E43 L= 37)   512     63        1 ( 747/ 255- 802/ 315)  23.2 E45-E38 S35-S31
  

    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