Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1105 1853 (N51E22 L=195)   394     46        0 (1078/1830-1134/1875) -15.9 E25-E18 N48-N52
      2 1258 1444 (N23W04 L=221)   246    103        0 (1218/1414-1304/1477)  -9.5 W01-W07 N21-N24
      3  571 1415 (N06E36 L=180)  1569    511        2 ( 358/1364- 762/1461)  14.0 E51-E23 N01-N09
      4 1294 1301 (N15W09 L=226)   977    244        1 (1226/1208-1361/1395) -45.0 W04-W13 N09-N19
      5 1648 1205 (N19W35 L=252)   394     76        0 (1626/1162-1673/1245) -53.1 W32-W37 N16-N20
      6  180 1128 (S16E64 L=153)   826    145        1 ( 159/1047- 204/1191)  77.5 E66-E61 S20-S11
      7 1722 1123 (N16W42 L=259)   564     81        2 (1672/1103-1770/1146)  12.5 W37-W46 N14-N17
      8  928  803 (S22E01 L=216)   239    136        0 ( 875/ 788-1005/ 826)  39.8 E03-W03 S23-S18
      9  229  735 (S40E60 L=156)   412    154        1 ( 213/ 666- 246/ 801)  52.1 E63-E55 S44-S35
     10  470  663 (S40E32 L=185)   888    119        2 ( 416/ 612- 524/ 725)  52.1 E37-E27 S44-S35
     11 1748  546 (S14W67 L=283)   860    100        2 (1716/ 504-1770/ 580) -63.4 W61-W71 S16-S12
  

    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