Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  826 1850 (N44E50 L=322)  1112    163        0 ( 771/1803- 872/1885)  -6.7 E60-E41 N40-N46
      2  545 1747 (N28E62 L=310)  1272    133        3 ( 493/1708- 595/1785)  56.3 E66-E58 N24-N31
      3  409 1623 (N17E64 L=308)   915     87        1 ( 360/1599- 456/1642)  -3.6 E73-E57 N16-N18
      4  946 1576 (N27E21 L=352)   637    156        0 ( 917/1514- 989/1638) -73.3 E22-E18 N22-N32
      5  692  927 (S17E17 L=355)   482    145        0 ( 639/ 898- 742/ 973) -14.0 E21-E13 S18-S15
      6  615  904 (S21E22 L=351)   297     87        1 ( 594/ 880- 648/ 944) -45.0 E23-E19 S22-S18
      7  827  696 (S28E03 L=  9)   870    295        0 ( 684/ 666- 952/ 728)  31.0 E12-W03 S32-S22
      8 1468  406 (S27W52 L= 64)   829     95        1 (1436/ 378-1506/ 446)  -7.1 W48-W56 S28-S24
      9 1306  320 (S37W47 L= 59)  3242    377        5 (1144/ 266-1435/ 402)  39.1 W34-W53 S44-S28
  

    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