Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1142 1736 (N44W04 L= 63)   283     43        0 (1116/1717-1166/1752)  36.9 W01-W05 N41-N45
      2 1266 1674 (N40W14 L= 73)   276     45        0 (1242/1645-1288/1694) -63.4 W11-W15 N37-N41
      3 1339 1666 (N40W19 L= 79)   436     58        0 (1312/1640-1373/1688)  38.7 W17-W22 N37-N41
      4  896 1617 (N32E13 L= 46)   849    166        0 ( 840/1543- 946/1671) -45.0 E18-E09 N27-N36
      5 1436 1623 (N37W27 L= 86)   565     85        1 (1404/1576-1466/1649) -80.5 W24-W29 N33-N39
      6 1655 1507 (N31W43 L=103)  1779    315        3 (1550/1436-1776/1565) -11.9 W34-W54 N27-N34
      7  775 1308 (N11E17 L= 42)   261     46        0 ( 739/1288- 799/1328) -14.0 E19-E15 N09-N12
      8  211 1260 (N07E61 L=359)   493     90        0 ( 185/1230- 236/1288) -38.7 E64-E57 N04-N08
      9  831  689 (S26E11 L= 49)   250     74        0 ( 810/ 664- 853/ 733)  63.4 E11-E09 S27-S23
     10  398  577 (S35E45 L= 14)   816    153        0 ( 336/ 549- 473/ 613) -15.3 E50-E38 S37-S32
  

    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