Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  956 1572 (N24E21 L= 23)   451    108        1 ( 918/1540-1012/1600)  -8.1 E23-E16 N22-N25
      2 1551 1524 (N41W20 L= 64)   463     59        2 (1521/1501-1582/1555) -21.8 W17-W23 N39-N42
      3  403 1438 (N02E52 L=352)   346     44        0 ( 378/1415- 431/1457)  33.7 E53-E49 N00-N03
      4 1763 1232 (N29W43 L= 87)   872    230        1 (1726/1127-1807/1321) -42.0 W39-W49 N24-N33
      5 1542 1011 (N08W30 L= 74)   319    102        0 (1513/ 985-1569/1036)  76.0 W28-W31 N06-N10
      6  597  989 (S19E24 L= 19)   455    176        1 ( 523/ 956- 678/1031)  41.6 E28-E19 S22-S14
      7 1860  964 (N17W59 L=103)   909    145        2 (1815/ 912-1890/1021)  76.0 W54-W63 N12-N20
      8  430  879 (S29E35 L=  9)   282     88        0 ( 397/ 848- 468/ 898)  26.6 E37-E31 S31-S27
      9  237  816 (S38E55 L=348)   472     60        1 ( 203/ 802- 279/ 836)   5.2 E60-E49 S39-S36
     10 1414  699 (S12W32 L= 76)   254     43        0 (1384/ 685-1444/ 714)  26.6 W30-W34 S13-S10
     11 1486  300 (S31W65 L=108)   946    104        2 (1453/ 267-1528/ 348)  -6.7 W57-W74 S32-S28
  

    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