Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1527 1755 (N33W66 L=  3)   905    152        2 (1458/1727-1602/1774) -32.5 W59-W73 N28-N35
      2  517 1678 (N59E28 L=269)   638     94        2 ( 484/1637- 554/1711)  12.5 E33-E23 N55-N61
      3 1203 1505 (N26W23 L=321)   264    128        0 (1155/1475-1262/1532)  -7.1 W19-W27 N25-N27
      4 1165 1303 (N16W15 L=312)   213     88        1 (1135/1282-1199/1326) -63.4 W12-W16 N13-N17
      5  990 1171 (N12W01 L=299)   552    197        3 ( 935/1101-1016/1261)  53.1 E03-W03 N08-N17
      6 1515  950 (S12W27 L=325)   230     74        2 (1501/ 910-1535/ 980)  53.1 W26-W29 S14-S10
      7 1842  645 (S42W61 L=358)   527     55        1 (1807/ 632-1869/ 659)  -8.1 W54-W68 S43-S39
      8  892  373 (S33E27 L=271)  1378    469        0 ( 736/ 273-1054/ 531)  12.7 E41-E10 S37-S26
      9  842  441 (S27E27 L=271)   270    113        0 ( 815/ 414- 870/ 477)  31.0 E28-E23 S28-S25
  

    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