Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1363 1774 (N53W48 L= 53)   956     84        1 (1321/1742-1405/1801)   6.3 W39-W57 N50-N55
      2 1250 1763 (N54W33 L= 38)   414     36        0 (1229/1736-1270/1786)  20.6 W29-W37 N51-N56
      3 1596 1666 (N41W64 L= 69)   664    102        1 (1565/1633-1647/1694) -38.7 W60-W68 N37-N43
      4  357 1538 (N41E60 L=304)   607     94        0 ( 308/1525- 402/1553)   3.8 E69-E54 N39-N41
      5  836 1452 (N33E11 L=354)  1129    251        2 ( 741/1392- 928/1540) -35.5 E18-E04 N28-N40
      6  661 1480 (N36E25 L=340)   377    125        0 ( 637/1425- 693/1535)  74.1 E26-E22 N32-N39
      7  868 1222 (N18E09 L=356)   364    144        0 ( 824/1176- 917/1272)  39.8 E11-E05 N15-N20
      8 1438 1116 (N08W27 L= 32)   225     57        1 (1420/1087-1455/1138)  56.3 W25-W28 N05-N08
      9 1722  846 (S12W48 L= 53)   585    106        1 (1685/ 801-1759/ 889)  45.0 W44-W51 S14-S09
     10  388  514 (S26E55 L=310)  1367    258        3 ( 313/ 434- 450/ 572) -52.1 E63-E49 S31-S21
     11  886  505 (S28E13 L=352)   299     98        0 ( 841/ 493- 943/ 517)   8.1 E16-E09 S28-S27
     12 1219  387 (S40W11 L= 15)   416     78        0 (1183/ 365-1260/ 416)  33.7 W07-W13 S41-S37
  

    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