Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1425 1774 (N55W53 L=164)  1100    105        3 (1393/1729-1468/1813)  14.5 W43-W70 N50-N57
      2 1435 1562 (N38W36 L=146)   379     94        1 (1400/1538-1463/1587)   8.1 W31-W38 N35-N39
      3  696 1440 (N31E23 L= 88)  2301    523        3 ( 509/1297- 887/1563) -27.2 E43-E08 N21-N39
      4 1726 1479 (N29W62 L=173)   542     52        0 (1705/1451-1741/1504)  53.1 W59-W66 N27-N31
      5 1243 1262 (N18W15 L=126)   467    204        0 (1211/1184-1299/1330)  53.1 W12-W19 N13-N22
      6  320 1154 (N12E50 L= 61)   296     46        0 ( 303/1131- 347/1173)  36.9 E51-E47 N10-N13
      7  946  998 (N03E05 L=106)   516    118        3 ( 897/ 954-1000/1050)  33.7 E07-E01 N00-N05
      8 1645  934 (S04W42 L=152)   284     62        2 (1613/ 913-1673/ 966) -36.9 W39-W44 S05-S01
      9  819  694 (S16E14 L= 97)  1739    508        5 ( 698/ 495- 898/ 903)  56.8 E25-E08 S29-S03
     10 1558  529 (S30W40 L=151)  4808    767        7 (1294/ 365-1768/ 750)  35.5 W20-W55 S42-S16
     11 1106  451 (S34W05 L=115)  1460    351        0 ( 974/ 354-1230/ 542)  36.4 E06-W13 S41-S27
  

    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