Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  967 1804 (N52E25 L=164)   716    127        1 ( 899/1788-1029/1820)  24.4 E30-E19 N49-N54
      2  583 1604 (N29E43 L=145)  2358    339        3 ( 458/1515- 721/1695) -12.0 E62-E29 N25-N33
      3  698 1338 (N14E25 L=163)   459    172        1 ( 644/1272- 731/1392) -33.7 E30-E21 N09-N16
      4  202 1230 (N00E64 L=124)   624     74        1 ( 175/1197- 237/1256)  36.9 E68-E60 S01-N02
      5 1313 1226 (N16W15 L=203)   301    105        0 (1259/1211-1370/1239)  23.2 W11-W18 N14-N17
      6 1766 1001 (N09W51 L=240)  2022    377        3 (1716/ 866-1820/1079) -65.2 W46-W58 N01-N14
      7 1614  965 (N05W39 L=227)   944    196        2 (1595/ 853-1640/1059) -84.8 W36-W40 S01-N10
      8 1844  965 (N08W60 L=249)   592     65        2 (1820/ 932-1870/ 997) -38.7 W57-W64 N06-N10
      9  550  887 (S15E28 L=161)   314    106        2 ( 499/ 861- 600/ 914)  29.7 E31-E24 S17-S12
     10  771  725 (S22E11 L=177)  1200    419        0 ( 579/ 678- 926/ 798)  28.6 E24-E02 S27-S15
     11 1745  648 (S12W59 L=247)  1918    306        0 (1662/ 535-1814/ 783)  76.8 W55-W61 S19-S02
     12  987  258 (S53W15 L=204)   844    120        1 ( 948/ 217-1028/ 304) -32.5 W10-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