Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  860 1790 (N56W07 L=304)   395     59        0 ( 825/1769- 896/1810)   7.1 W03-W11 N54-N57
      2 1433 1714 (N35W49 L=346)   745     99        0 (1385/1684-1484/1745)   4.4 W43-W56 N33-N36
      3 1616 1675 (N28W64 L=  1)  1533    127        1 (1552/1645-1685/1692)  -9.5 W57-W77 N26-N30
      4 1595 1520 (N19W50 L=347)   374     46        0 (1563/1497-1623/1538) -36.9 W47-W51 N17-N20
      5  778 1352 (N25E09 L=288)  2476    607        3 ( 491/1275-1034/1416)  -7.7 E29-W08 N22-N29
      6  228 1338 (N35E58 L=239)   381     44        0 ( 211/1315- 256/1354)  11.3 E60-E54 N33-N36
      7 1825 1133 (S07W57 L=354)   493     31        0 (1801/1109-1844/1160)  21.8 W54-W59 S08-S05
      8  381 1004 (N11E41 L=256)   246     34        0 ( 357/ 986- 410/1022) -26.6 E43-E39 N10-N12
      9 1486  857 (S17W25 L=322)   606    177        0 (1444/ 780-1508/ 943)  60.9 W21-W27 S20-S11
     10 1573  741 (S26W31 L=327)   259     89        0 (1553/ 698-1596/ 784)  45.0 W28-W32 S28-S23
     11 1315  732 (S21W12 L=309)   234     76        0 (1291/ 698-1338/ 769)  53.1 W10-W14 S22-S18
     12  521  402 (S27E51 L=246)  1706    328        5 ( 423/ 342- 612/ 508)  16.3 E61-E37 S30-S21
  

    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