Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1121 1891 (N62E28 L= 35)  1192    109        1 (1082/1842-1167/1926) -10.0 E37-E19 N55-N66
      2 1182 1751 (N49E09 L= 54)  2542    295        9 (1093/1658-1271/1827) -21.8 E19-W01 N42-N56
      3  399 1541 (N15E57 L=  6)  1858    243        5 ( 366/1442- 448/1648) -45.0 E67-E52 N08-N22
      4  488 1520 (N16E48 L= 15)   341     68        1 ( 452/1496- 511/1546)  53.1 E49-E45 N13-N17
      5  625 1202 (S00E27 L= 36)  1312    295        2 ( 542/1086- 710/1312)  71.6 E30-E24 S08-N07
      6 1384 1262 (N19W17 L= 80)   228     99        0 (1336/1239-1425/1294)  -9.5 W13-W19 N18-N20
      7  952 1094 (N00E06 L= 57)   239    102        0 ( 936/1041- 963/1146) -80.5 E07-E04 S03-N03
      8 1205 1025 (N02W10 L= 73)   213     85        0 (1190/ 983-1226/1066)  90.0 W09-W10 S01-N04
      9  846  990 (S08E10 L= 53)   320    112        0 ( 800/ 946- 895/1023)  56.3 E11-E07 S11-S05
     10 1796  772 (N01W59 L=122)  3633    325        4 (1732/ 637-1853/ 919) -31.8 W49-W71 S06-N07
     11 1458  850 (S03W30 L= 93)   571     97        1 (1418/ 810-1502/ 887)  59.0 W27-W32 S05-N00
  

    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