Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1101 1577 (N40W05 L=341)   394    111        3 (1081/1538-1122/1613) -80.5 W03-W07 N36-N42
      2 1652 1431 (N29W48 L= 23)   569     97        0 (1622/1401-1685/1455)  15.9 W44-W51 N26-N30
      3 1788 1419 (N27W61 L= 37)   858    141        1 (1752/1395-1829/1435) -12.5 W57-W66 N25-N28
      4 1168 1315 (N22W09 L=345)   601    111        2 (1136/1261-1203/1379) -60.3 W07-W11 N18-N25
      5 1430 1244 (N17W27 L=  2)   267     78        0 (1390/1205-1457/1287) -51.3 W24-W28 N14-N19
      6 1847 1175 (N12W63 L= 39)  1300    190        0 (1800/1116-1887/1237) -38.7 W57-W68 N07-N15
      7  190 1152 (N08E65 L=271)  1842    190        2 ( 166/1059- 217/1235) -84.8 E68-E61 N02-N13
      8 1715 1124 (N09W48 L= 24)   281     30        0 (1698/1103-1730/1147) -71.6 W46-W49 N07-N10
      9 1522 1050 (N05W32 L=  8)   256     58        1 (1500/1016-1543/1086) -68.2 W30-W33 N02-N07
     10 1633  982 (N01W41 L= 17)  1018    184        2 (1552/ 938-1701/1038)  20.0 W34-W46 S01-N03
     11  198  893 (S08E63 L=273)   956    330        3 ( 170/ 785- 222/ 983) -85.2 E66-E59 S14-S01
     12 1446  459 (S34W33 L=  9)   393     69        0 (1410/ 436-1473/ 489)  36.9 W30-W35 S35-S31
  

    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