Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1112 1580 (N40W06 L=341)   396     69        1 (1082/1544-1130/1612)  68.2 W03-W07 N37-N42
      2 1729 1426 (N28W55 L= 30)  1730    374        3 (1636/1396-1836/1454)   0.0 W45-W67 N26-N30
      3  456 1392 (N25E42 L=293)   288     98        1 ( 424/1369- 493/1420) -23.2 E45-E38 N23-N26
      4 1190 1306 (N21W11 L=345)   756    168        2 (1157/1230-1231/1394) -73.3 W08-W13 N16-N26
      5  198 1150 (N08E63 L=271)  1846    203        4 ( 172/1056- 229/1239) -74.7 E66-E59 N02-N13
      6 1854 1181 (N12W64 L= 39)  1172    162        0 (1804/1121-1885/1236) -37.9 W58-W67 N08-N15
      7 1541 1045 (N05W34 L=  8)   315     71        0 (1517/1007-1565/1083) -53.1 W31-W35 N02-N06
      8 1649  985 (N01W42 L= 17)  1092    191        2 (1571/ 946-1716/1042)  20.0 W35-W47 S01-N04
      9 1470  462 (S34W35 L= 10)   345     42        0 (1447/ 441-1495/ 481)  56.3 W33-W37 S35-S32
  

    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