Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1729 1547 (N13W68 L=131)  1155     71        1 (1702/1510-1754/1585)  68.2 W63-W72 N10-N15
      2 1559 1540 (N18W51 L=114)   418     89        1 (1524/1509-1588/1571) -78.7 W49-W52 N15-N20
      3 1246 1467 (N24W26 L= 89)  2421    768        7 ( 962/1419-1604/1526) -29.4 W10-W50 N12-N34
      4  288 1474 (N52E56 L=  7)  2147    192        5 ( 195/1441- 371/1511)   0.0 E83-E42 N48-N53
      5  878 1401 (N30W02 L= 65)   255     58        0 ( 854/1383- 910/1421) -53.1 W00-W03 N28-N32
      6  501 1188 (N27E29 L= 34)   296     73        1 ( 461/1176- 540/1206)  -9.5 E31-E25 N26-N28
      7  427  998 (N18E38 L= 25)  1451    392        4 ( 396/ 850- 491/1189)  84.6 E41-E33 N08-N29
      8 1594 1103 (S06W35 L= 98)   245     42        1 (1571/1074-1613/1124)  71.6 W34-W37 S07-S04
      9 1818  836 (S29W50 L=113)  1409    335        3 (1769/ 683-1862/ 989)  59.5 W44-W57 S38-S20
     10 1258  677 (S21W04 L= 67)   565    230        0 (1149/ 665-1364/ 694) -21.0 E02-W11 S23-S17
     11  458  494 (S12E52 L= 10)   555     82        0 ( 430/ 460- 496/ 521)  12.5 E56-E47 S13-S10
     12  668  356 (S25E45 L= 18)   673    133        0 ( 611/ 318- 696/ 414)  18.4 E51-E39 S27-S21
  

    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