Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1371 1802 (N49W61 L=329)   953    107        0 (1333/1766-1415/1832)   2.7 W51-W72 N46-N50
      2 1505 1735 (N40W65 L=333)  1527    128        3 (1463/1690-1551/1786)  53.1 W57-W71 N36-N44
      3 1493 1637 (N34W52 L=320)   777    167        0 (1436/1561-1527/1686)  77.5 W48-W55 N29-N38
      4 1449 1606 (N34W46 L=314)   396     84        0 (1418/1569-1483/1642) -90.0 W45-W47 N30-N36
      5  593 1488 (N43E26 L=242)   299     72        1 ( 562/1475- 634/1503)  -8.1 E28-E21 N41-N44
      6 1747 1369 (N12W60 L=328)   801    117        0 (1705/1327-1776/1412) -63.4 W57-W62 N09-N15
      7 1118 1224 (N17W09 L=277)  1121    376        1 ( 940/1164-1272/1278) -33.7 E01-W17 N10-N22
      8 1603  991 (S06W36 L=304)   617    148        0 (1570/ 931-1636/1057)  45.0 W33-W39 S09-S02
      9 1912  993 (S14W70 L=338)  1065    104        1 (1889/ 951-1936/1035)  23.2 W65-W75 S15-S11
     10  312  589 (S12E61 L=207)   407     49        1 ( 290/ 563- 331/ 611) -56.3 E63-E58 S13-S10
     11 1589  489 (S40W35 L=303)   444     58        2 (1568/ 457-1612/ 523)  90.0 W32-W37 S42-S36
     12 1182  229 (S54E05 L=262)   377     40        0 (1163/ 212-1209/ 250) -59.0 E07-E02 S55-S50
  

    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