Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1166 1681 (N41W33 L=324)   409     50        0 (1135/1659-1200/1699) -21.8 W30-W35 N39-N43
      2 1028 1524 (N36W15 L=306)  1125    290        2 ( 916/1479-1168/1548) -12.5 W06-W24 N32-N39
      3 1834 1072 (S13W55 L=346)  1257    180        2 (1800/ 983-1851/1167)  63.4 W50-W59 S17-S07
      4 1338  973 (S03W16 L=307)   591    181        1 (1307/ 893-1372/1059)  45.0 W13-W20 S07-N00
      5  359 1021 (N21E44 L=247)   321     81        0 ( 327/1001- 397/1054) -33.7 E46-E40 N19-N23
      6  331  944 (N17E47 L=244)   939    247        3 ( 240/ 889- 413/ 984) -10.6 E56-E40 N12-N18
      7  317  856 (N12E50 L=241)   593     84        0 ( 293/ 811- 341/ 897) -78.7 E52-E47 N09-N14
      8 1589  621 (S32W27 L=318)  1994    310        5 (1469/ 521-1675/ 731)  21.3 W16-W34 S36-S26
      9  491  378 (S20E56 L=235)  3361    444        8 ( 388/ 267- 620/ 555) -90.0 E66-E49 S30-S07
     10  337  427 (S15E71 L=220)   908     99        4 ( 303/ 397- 369/ 455) -12.5 E78-E65 S17-S12
     11  895  207 (S42E38 L=253)  2385    302        3 ( 792/ 121-1021/ 288)  12.4 E64-E23 S48-S37
  

    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