Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  970 1802 (N52E24 L=164)   695    127        1 ( 903/1787-1034/1816)  24.4 E29-E18 N49-N54
      2  589 1603 (N29E43 L=145)  2446    352        1 ( 458/1511- 731/1694) -11.6 E62-E28 N24-N33
      3  706 1336 (N14E25 L=163)   442    155        2 ( 658/1267- 735/1387) -36.9 E29-E21 N09-N16
      4  204 1231 (N00E64 L=124)   755     75        1 ( 172/1196- 242/1256)  45.0 E69-E59 S01-N02
      5 1327 1225 (N16W16 L=204)   319    102        0 (1272/1210-1383/1240)  23.2 W12-W19 N14-N17
      6 1763 1018 (N10W51 L=239)  1824    312        2 (1721/ 909-1807/1098) -45.0 W46-W56 N04-N14
      7 1621  959 (N05W39 L=227)  1002    204        0 (1601/ 853-1645/1053) -84.8 W37-W41 S01-N10
      8 1840  967 (N08W60 L=248)   677     83        0 (1801/ 930-1871/ 994) -18.4 W55-W64 N06-N10
      9  784  724 (S22E11 L=177)  1240    417        0 ( 587/ 676- 934/ 794)  27.6 E24-E01 S27-S15
     10 1736  623 (S13W59 L=247)  1792    265        0 (1667/ 522-1801/ 733)  77.9 W55-W62 S20-S05
     11 1791  600 (S14W68 L=256)   529    151        2 (1765/ 547-1815/ 664) -69.4 W65-W71 S17-S09
     12  991  258 (S53W15 L=203)   818    104        2 ( 948/ 219-1032/ 302) -35.0 W11-W21 S56-S49
  

    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