Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  708 1658 (N25E40 L=223)  2052    402        5 ( 579/1540- 852/1739)  66.0 E43-E32 N15-N33
      2  893 1640 (N29E26 L=236)   860    160        3 ( 848/1559- 925/1702) -45.0 E31-E21 N24-N33
      3 1351 1393 (N26W10 L=272)  1890    480        3 (1242/1169-1405/1563) -51.3 W00-W17 N12-N35
      4 1253 1181 (N11W09 L=271)   203     79        0 (1222/1153-1279/1213) -26.6 W06-W10 N09-N11
      5 1655 1017 (N12W39 L=301)   332    111        0 (1642/ 960-1666/1073) -74.1 W37-W39 N08-N15
      6  309 1045 (S19E48 L=214)   495    112        0 ( 263/1009- 359/1072)  32.0 E52-E43 S22-S16
      7  686  908 (S19E17 L=245)   396    190        1 ( 649/ 850- 737/ 965) -81.9 E19-E14 S21-S14
      8  333  914 (S26E44 L=218)   673    160        1 ( 259/ 893- 411/ 949)  19.7 E51-E36 S28-S22
      9 1160  863 (S09W12 L=274)   373    164        0 (1093/ 808-1220/ 893)   0.0 W07-W16 S11-S08
     10 1781  776 (N03W57 L=319)   609    124        1 (1731/ 737-1819/ 828)  45.0 W53-W59 N00-N06
     11  646  436 (S51E10 L=252)  1759    356        2 ( 514/ 357- 798/ 525)  -7.6 E24-W06 S53-S47
  

    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