Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1048 1538 (N28E02 L=358)   344     80        1 (1024/1506-1085/1562) -33.7 E03-W00 N25-N29
      2 1728 1445 (N27W49 L= 49)   901    105        1 (1683/1391-1760/1480) -56.3 W45-W51 N23-N29
      3 1643 1405 (N23W40 L= 40)   303     89        0 (1604/1376-1684/1425) -18.4 W37-W43 N21-N24
      4 1416 1251 (N12W23 L= 23)   410    170        2 (1338/1206-1490/1304) -24.4 W17-W28 N09-N14
      5 1703 1120 (N06W44 L= 44)   911    268        2 (1631/1065-1779/1214) -30.3 W38-W51 N03-N11
      6 1476 1104 (N03W27 L= 27)   282    120        1 (1422/1075-1526/1146)  35.5 W23-W30 N01-N06
      7  216 1005 (S08E57 L=303)   508     79        3 ( 192/ 967- 239/1042)  68.2 E59-E54 S10-S05
      8  831  974 (S08E11 L=349)   284     74        1 ( 815/ 933- 849/1021)  68.2 E12-E10 S10-S05
      9  462  885 (S15E36 L=324)   308    113        1 ( 441/ 841- 480/ 947) -90.0 E37-E34 S17-S11
     10 1222  907 (S10W13 L= 13)   212    115        0 (1167/ 890-1287/ 919)  15.9 W09-W16 S11-S09
     11  677  867 (S16E21 L=339)   430    109        2 ( 635/ 844- 740/ 894)  15.9 E23-E16 S16-S13
     12  199  796 (S21E60 L=300)   972    229        3 ( 186/ 697- 222/ 890)  71.6 E63-E56 S26-S14
     13 1744  761 (S15W52 L= 52)   499    111        0 (1716/ 702-1774/ 823)  76.0 W50-W53 S18-S10
     14 1419  783 (S16W27 L= 27)   266     97        0 (1388/ 750-1460/ 806)  21.8 W24-W29 S18-S14
     15 1108  672 (S25W08 L=  8)   414    114        2 (1055/ 624-1154/ 715)  45.0 W04-W10 S28-S22
     16 1302  561 (S31W23 L= 23)  2652    636        4 ( 998/ 431-1592/ 692)  30.1 W02-W40 S42-S20
     17  503  555 (S36E37 L=323)   277    154        0 ( 471/ 511- 534/ 603)  45.0 E40-E33 S39-S32
  

    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