Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  971 1827 (N64E06 L=  5)   529     66        2 ( 938/1805-1002/1847)  24.0 E10-E01 N61-N66
      2  277 1457 (N30E67 L=304)   520     62        0 ( 252/1424- 293/1477) -18.4 E72-E62 N27-N31
      3  864 1372 (N26E11 L=  0)   222     70        0 ( 833/1351- 905/1388)  21.8 E12-E07 N24-N26
      4 1214 1331 (N23W13 L= 24)   328    134        0 (1160/1284-1259/1375) -45.0 W09-W15 N19-N25
      5 1758  856 (S09W53 L= 64)   798    144        2 (1719/ 813-1820/ 884)  24.0 W49-W59 S11-S06
      6  529  820 (S09E33 L=338)   287     42        1 ( 505/ 792- 554/ 837)  14.0 E34-E30 S10-S08
      7 1629  763 (S14W42 L= 53)   685     87        1 (1591/ 727-1662/ 811)  71.6 W38-W44 S16-S10
      8 1730  749 (S15W52 L= 63)  1326    199        3 (1668/ 692-1800/ 809)  26.6 W46-W58 S19-S11
      9 1229  679 (S18W13 L= 24)   453    111        1 (1183/ 647-1279/ 714)  33.7 W10-W16 S20-S15
     10 1682  610 (S25W51 L= 62)  1395    266        2 (1621/ 560-1747/ 683)  33.7 W44-W56 S28-S19
     11  690  480 (S32E26 L=345)  2664    370        8 ( 592/ 398- 804/ 605) -32.7 E32-E17 S38-S22
  

    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