Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  360 1575 (N26E61 L=148)  1116    176        3 ( 327/1529- 394/1635) -54.5 E63-E57 N23-N30
      2  543 1389 (N17E37 L=171)   760    236        2 ( 522/1291- 560/1511) -70.3 E40-E34 N10-N24
      3 1509 1414 (N31W30 L=238)   698    196        0 (1427/1359-1587/1452) -14.0 W23-W35 N28-N33
      4 1297 1349 (N24W14 L=222)   207     65        2 (1281/1316-1309/1388) -76.0 W12-W14 N22-N26
      5 1129 1296 (N19W03 L=212)   197     51        0 (1104/1274-1155/1318) -26.6 W01-W05 N17-N19
      6 1842  994 (N08W61 L=269)  5275    819       12 (1785/ 756-1906/1180) -80.9 W53-W70 S05-N20
      7  194  892 (S18E61 L=148)  2874    533        8 ( 151/ 799- 252/ 985)  15.9 E66-E54 S24-S11
      8  346  934 (S13E45 L=163)   343     59        3 ( 325/ 907- 367/ 969)  90.0 E47-E43 S15-S11
      9  986  863 (S10E00 L=208)   187     74        0 ( 969/ 825-1000/ 901)  90.0 W00-W00 S12-S07
     10  298  778 (S24E50 L=159)  1144    301        3 ( 240/ 702- 385/ 828)  37.6 E55-E41 S29-S19
     11 1611  698 (S12W44 L=253)  1467    375        4 (1491/ 606-1716/ 785)  47.5 W37-W53 S19-S06
     12  842  724 (S20E08 L=201)   376    129        1 ( 788/ 699- 913/ 760)  -7.1 E11-E03 S21-S18
     13 1423  737 (S12W29 L=238)   411    125        2 (1359/ 724-1478/ 752)  14.0 W25-W33 S13-S10
     14 1259  636 (S20W20 L=229)   320     96        1 (1215/ 603-1305/ 677)  51.3 W18-W22 S22-S17
     15 1629  633 (S16W48 L=256)   602     78        2 (1600/ 610-1664/ 656)  45.0 W45-W51 S17-S14
     16 1595  527 (S23W49 L=257)   540    157        2 (1557/ 472-1634/ 569) -45.0 W45-W52 S26-S20
     17  604  533 (S37E25 L=184)   412    105        0 ( 562/ 517- 641/ 562)  29.7 E28-E21 S38-S33
  

    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