Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1594 1671 (N40W65 L=165)   476     32        0 (1571/1648-1617/1688)   5.2 W59-W70 N39-N41
      2 1343 1641 (N42W35 L=135)   417     62        2 (1316/1606-1368/1664) -78.7 W32-W37 N39-N44
      3  717 1457 (N35E19 L= 81)   254     49        1 ( 699/1430- 736/1486) -63.4 E20-E17 N33-N37
      4  365 1396 (N32E50 L= 50)   851    120        3 ( 321/1360- 415/1425)  18.4 E55-E45 N29-N34
      5  559 1382 (N31E31 L= 69)   337     91        0 ( 520/1362- 606/1402)   8.1 E34-E27 N29-N31
      6  255 1133 (N15E57 L= 44)   575     59        2 ( 229/1112- 287/1161)  -8.1 E59-E52 N14-N17
      7  420  634 (S16E46 L= 54)   419    135        2 ( 363/ 599- 468/ 666) -29.7 E50-E42 S18-S13
      8 1496  505 (S34W31 L=132)   282     83        0 (1455/ 490-1536/ 530)   8.1 W27-W34 S34-S32
      9 1381  414 (S39W23 L=123)   543    147        0 (1317/ 387-1434/ 441)   0.0 W17-W27 S41-S37
     10 1427  306 (S51W32 L=132)   419     39        1 (1401/ 290-1451/ 324)   0.0 W28-W34 S52-S48
  

    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