Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1059 1626 (N39W23 L=161)   898    215        1 (1014/1541-1111/1694)  33.7 W16-W29 N33-N41
      2 1473 1597 (N24W47 L=185)  1054    152        2 (1406/1562-1544/1623) -45.0 W41-W51 N19-N26
      3  519 1244 (N31E26 L=112)   876    349        1 ( 382/1205- 677/1278) -22.2 E36-E14 N25-N36
      4 1524 1159 (S01W32 L=170)  1007    244        3 (1425/1121-1619/1195) -22.6 W25-W38 S05-N01
      5 1689  982 (S16W39 L=177)   610    122        2 (1639/ 951-1727/1020)   0.0 W35-W42 S17-S14
      6 1417  831 (S16W17 L=155)   613    196        0 (1394/ 752-1452/ 934)  74.7 W15-W20 S21-S10
      7 1093  684 (S15E05 L=133)   706    159        0 (1002/ 644-1169/ 715) -45.0 E09-E01 S19-S11
      8 1184  562 (S25E04 L=135)   217     75        0 (1170/ 525-1199/ 596)  68.2 E04-E02 S27-S22
      9  783  400 (S24E34 L=104)   448    108        3 ( 755/ 359- 814/ 442)  80.5 E35-E31 S27-S21
     10 1034  218 (S46E30 L=108)   990    142        2 (1008/ 172-1060/ 261)  35.0 E35-E25 S49-S41
  

    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