Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  178 1265 (N41E68 L= 53)  1008    116        1 ( 164/1204- 192/1320) -69.4 E72-E64 N36-N44
      2  323 1056 (N25E46 L= 74)   978    234        0 ( 236/ 976- 405/1117) -38.7 E55-E40 N18-N30
      3 1870 1030 (S19W59 L=179)   436    104        1 (1852/ 989-1885/1082)  38.7 W56-W62 S20-S16
      4  140 1024 (N26E71 L= 50)  1462    135        3 ( 111/ 969- 172/1079)  31.0 E77-E64 N22-N29
      5  171  931 (N20E66 L= 54)  2288    265        7 ( 124/ 815- 218/1026)  72.9 E75-E58 N12-N25
      6  239  926 (N19E57 L= 63)   724    157        4 ( 220/ 845- 266/ 996)  63.4 E59-E55 N14-N22
      7 1880  936 (S25W59 L=179)   861    122        3 (1849/ 897-1923/ 978) -33.7 W55-W67 S29-S21
      8 1446  867 (S14W20 L=140)   888    256        0 (1341/ 822-1541/ 933)   3.8 W12-W27 S15-S12
      9 1781  691 (S37W46 L=166)   420     41        0 (1757/ 672-1802/ 714)   0.0 W42-W48 S38-S35
     10  839  381 (S26E32 L= 89)  3329    876        1 ( 667/ 244-1113/ 607)  11.8 E49-E06 S34-S18
     11 1177  492 (S29E06 L=114)   615    268        0 (1093/ 426-1214/ 607)  31.6 E13-W00 S31-S22
     12 1654  557 (S41W32 L=152)   601    108        0 (1609/ 526-1706/ 593)   5.2 W26-W37 S42-S39
  

    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