Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  464 1521 (N18E45 L=351)  2296    420       11 ( 334/1457- 617/1588)  -6.1 E59-E31 N14-N22
      2 1802 1465 (N35W57 L= 93)   535     49        2 (1775/1444-1833/1488)   6.3 W52-W62 N33-N37
      3  416 1431 (N12E46 L=351)  1047    142        3 ( 350/1378- 472/1485)  40.6 E50-E42 N08-N15
      4  346 1457 (N13E53 L=343)   529     65        0 ( 306/1435- 385/1476)   8.1 E57-E49 N11-N14
      5  235 1368 (N07E61 L=335)  3089    300       11 ( 174/1287- 318/1448) -84.3 E70-E53 N02-N12
      6 1090 1094 (S01W03 L= 39)   220     93        0 (1079/1048-1101/1140) -76.0 W02-W03 S03-N01
      7  384 1047 (S11E42 L=355)   777    195        2 ( 320/ 976- 431/1132)  56.3 E46-E38 S15-S06
      8  265 1066 (S11E52 L=344)   568    102        2 ( 240/1025- 283/1114)  63.4 E54-E50 S13-S07
      9 1941 1028 (N10W70 L=106)  1087    133        2 (1917/ 962-1964/1086)  63.4 W66-W74 N05-N13
     10  781  711 (S27E11 L= 25)  2270    661        3 ( 530/ 622-1034/ 916)  35.7 E30-W02 S35-S12
     11  291  825 (S25E50 L=346)  2444    419        3 ( 136/ 762- 475/ 870)   4.6 E70-E33 S29-S21
     12  687  848 (S20E19 L= 18)   306    119        0 ( 620/ 833- 739/ 863)   7.1 E23-E15 S21-S18
     13 1179  726 (S21W14 L= 50)   651    180        1 (1160/ 638-1197/ 811) -73.3 W11-W15 S26-S16
     14 1638  622 (S20W49 L= 85)   840    232        2 (1592/ 560-1719/ 693) -41.2 W42-W54 S23-S16
     15 1401  642 (S23W30 L= 66)   360     99        0 (1355/ 616-1436/ 682)  45.0 W27-W31 S24-S19
     16 1216  374 (S43W27 L= 64)  2685    400        8 (1024/ 325-1404/ 432)   8.6 W11-W44 S46-S39
  

    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