Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1595 1644 (N22W63 L=209)  1265    126        1 (1549/1598-1650/1686)   0.0 W54-W74 N19-N22
      2  358 1538 (N57E46 L=100)   488     48        0 ( 331/1518- 395/1554)  -5.2 E51-E40 N55-N58
      3  775 1453 (N37E03 L=144)   371    125        1 ( 703/1436- 829/1470) -26.6 E08-W00 N34-N39
      4 1820 1399 (N02W68 L=214)  1356    145        3 (1783/1335-1852/1455) -81.9 W63-W71 S02-N05
      5 1779 1033 (S16W48 L=195)   345     59        2 (1753/1011-1801/1067)  26.6 W45-W50 S17-S13
      6  154  866 (N17E70 L= 77)  1410    126        5 ( 124/ 810- 183/ 937) -74.1 E75-E65 N13-N20
      7 1385  825 (S15W15 L=161)   989    419        2 (1245/ 716-1461/ 934)  21.3 W04-W22 S18-S10
      8  994  735 (S09E09 L=137)   661    260        2 ( 918/ 625-1050/ 839) -90.0 E11-E07 S16-S01
      9 1418  698 (S24W14 L=160)   372    202        0 (1373/ 659-1459/ 739)   8.1 W10-W17 S26-S21
     10 1446  616 (S30W14 L=161)   360     77        0 (1405/ 597-1483/ 634) -26.6 W11-W17 S31-S27
     11  953  578 (S17E17 L=130)   220    102        0 ( 931/ 546- 967/ 629)  53.1 E18-E15 S19-S14
     12 1098  411 (S32E14 L=132)  1003    202        0 (1047/ 323-1132/ 482)  69.0 E19-E11 S38-S25
  

    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