Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  534 1593 (N14E48 L=342)   704    141        2 ( 486/1546- 598/1650)  63.4 E50-E46 N10-N18
      2  649 1573 (N15E39 L=351)  1155    270        4 ( 588/1487- 703/1635) -56.3 E44-E35 N10-N19
      3  515 1456 (N05E43 L=346)   930    287        3 ( 454/1412- 564/1529) -29.1 E49-E39 N01-N08
      4 1619 1479 (N40W26 L= 56)  1003    179        3 (1551/1428-1691/1518)   8.7 W20-W33 N37-N43
      5 1779 1372 (N39W43 L= 73)   383    110        1 (1747/1352-1809/1393)  35.5 W39-W46 N36-N41
      6 1651 1012 (N12W37 L= 67)   334     84        1 (1622/ 973-1681/1043) -36.9 W34-W39 N09-N13
      7  734  979 (S16E15 L= 15)   300     45        2 ( 714/ 953- 757/1010) -33.7 E16-E13 S17-S14
      8  849  936 (S15E07 L= 22)   212     53        1 ( 819/ 915- 880/ 963)   0.0 E09-E05 S16-S13
      9 1371  385 (S30W45 L= 75)   433    130        0 (1320/ 357-1414/ 416)   6.3 W39-W48 S30-S29
  

    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