Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1078 1773 (N52E12 L=162)   697    143        0 (1008/1751-1145/1796)  28.6 E17-E06 N48-N54
      2  733 1566 (N29E30 L=144)  2005    393        1 ( 574/1465- 897/1643)  -9.5 E45-E15 N24-N32
      3 1613 1564 (N45W42 L=215)   760    115        1 (1548/1541-1655/1593)  12.1 W33-W47 N41-N47
      4  902 1302 (N15E12 L=162)   395    138        0 ( 850/1249- 941/1347) -29.7 E15-E08 N12-N16
      5  298 1280 (N05E55 L=119)   314     45        1 ( 282/1248- 318/1303) -45.0 E56-E52 N03-N06
      6 1887 1031 (N13W66 L=239)  4403    271        7 (1845/ 896-1937/1129) -54.5 W59-W75 N04-N18
      7 1785  939 (N06W54 L=228)   887    165        1 (1768/ 863-1800/1022) -71.6 W51-W57 N01-N10
      8  236  954 (S15E56 L=118)   603    100        1 ( 219/ 902- 255/1002) -90.0 E57-E53 S18-S12
      9  923  853 (S12E04 L=170)   404    125        0 ( 872/ 822- 968/ 895)  45.0 E06-E01 S14-S08
     10  998  670 (S22W04 L=177)  1810    510        0 ( 772/ 625-1168/ 772)  32.5 E10-W12 S27-S13
     11 1122  228 (S53W30 L=203)   843    103        3 (1089/ 189-1157/ 271) -32.5 W24-W36 S56-S48
  

    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