Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  524 1649 (N43E50 L= 42)   509     39        0 ( 501/1625- 549/1676) -16.7 E55-E45 N41-N44
      2 1405 1571 (N40W30 L=122)   410     64        0 (1370/1547-1433/1593) -36.9 W27-W32 N38-N41
      3  292 1517 (N31E71 L= 21)  1051    127        1 ( 276/1460- 310/1578) -22.8 E81-E62 N27-N35
      4 1674 1450 (N31W53 L=144)   645    105        2 (1629/1409-1712/1488) -39.8 W48-W56 N28-N33
      5 1205 1368 (N25W12 L=103)   333     98        0 (1170/1330-1249/1400)  38.7 W09-W14 N22-N27
      6 1444 1276 (N19W28 L=120)   527    190        1 (1368/1217-1503/1332) -41.6 W22-W31 N15-N23
      7 1844 1232 (N16W65 L=157)  1845    210        3 (1805/1153-1874/1314)  59.0 W59-W71 N11-N21
      8 1747  855 (S07W52 L=144)   505     33        1 (1724/ 830-1772/ 878)  14.0 W49-W54 S08-S05
      9 1043  837 (S08W02 L= 93)   234    113        0 ( 998/ 818-1085/ 869)  21.8 E01-W04 S09-S06
     10  894  746 (S15E08 L= 84)  1145    270        4 ( 851/ 645- 936/ 867)  82.4 E10-E04 S21-S06
     11  444  786 (S14E39 L= 53)   282     66        1 ( 422/ 767- 481/ 807)  14.0 E41-E36 S15-S12
     12  177  745 (S19E72 L= 20)   705     34        1 ( 151/ 725- 203/ 763)   7.6 E81-E66 S20-S17
     13  273  686 (S22E59 L= 33)   450     42        0 ( 249/ 662- 293/ 713) -53.1 E61-E56 S23-S19
     14 1195  666 (S19W12 L=104)   319     66        0 (1165/ 635-1219/ 696) -45.0 W09-W13 S21-S17
     15 1460  528 (S29W34 L=126)  1566    260        2 (1417/ 431-1505/ 648) -59.0 W30-W41 S35-S20
     16  763  599 (S25E17 L= 75)   647    167        3 ( 689/ 568- 831/ 626)   5.7 E22-E12 S27-S22
     17 1671  494 (S31W57 L=149)   614     93        1 (1626/ 482-1713/ 512)  -8.7 W50-W63 S32-S29
  

    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