Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1459 1725 (N32W55 L=246)  1133    203        3 (1390/1699-1530/1745) -24.4 W49-W61 N29-N34
      2 1057 1446 (N29W15 L=206)   312     89        0 (1021/1415-1103/1469)   8.1 W11-W18 N28-N30
      3 1512 1428 (N14W42 L=233)  2852    690        7 (1249/1373-1735/1461) -23.2 W24-W59 N05-N21
      4  748 1280 (N28E10 L=181)   506    124        2 ( 711/1239- 800/1339) -81.9 E11-E06 N25-N32
      5 1870 1203 (S10W63 L=254)  1207    211        5 (1839/1117-1895/1265)  83.7 W60-W67 S14-S05
      6 1168 1160 (N10W12 L=203)   357    126        1 (1155/1091-1189/1223)  60.3 W09-W14 N06-N13
      7 1794 1115 (S11W52 L=243)  1231    259        4 (1746/1072-1863/1162) -41.2 W47-W58 S14-S07
      8 1709 1093 (S09W44 L=235)   266     41        1 (1695/1065-1724/1117) -90.0 W42-W45 S11-S07
      9 1207 1039 (N03W10 L=201)   428    109        2 (1155/1017-1260/1083) -51.3 W07-W12 N00-N05
     10 1812 1018 (S17W51 L=242)   456    113        1 (1788/ 984-1845/1065) -80.5 W48-W53 S20-S14
     11  348  959 (N19E45 L=146)   385     82        1 ( 335/ 914- 360/1001)  78.7 E46-E43 N16-N21
     12 1567  797 (S22W26 L=217)  1828    387        2 (1480/ 648-1605/ 954)  87.3 W22-W28 S31-S09
     13 1379  672 (S23W11 L=202)   333    123        0 (1354/ 610-1420/ 723)  39.8 W08-W14 S26-S21
     14  511  547 (S08E45 L=146)   427    135        0 ( 461/ 519- 572/ 583) -39.8 E47-E41 S10-S05
     15  293  535 (S05E67 L=124)  1222    145        2 ( 259/ 496- 333/ 581) -90.0 E71-E62 S08-S02
     16  718  467 (S18E35 L=156)   503    171        1 ( 677/ 405- 749/ 535)  31.0 E40-E29 S20-S14
     17  789  273 (S33E41 L=150)  1556    383        3 ( 707/ 200- 927/ 346)   9.5 E52-E28 S36-S30
  

    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