Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1701 (N52W06 L= 47)   439     81        1 ( 925/1675- 996/1721)  15.9 W02-W09 N50-N53
      2 1771 1486 (N22W71 L=112)   720    110        0 (1715/1472-1820/1502) -17.4 W63-W82 N19-N24
      3  171 1224 (N24E71 L=330)   823    106        1 ( 154/1172- 190/1284) -41.2 E75-E65 N21-N28
      4  429 1035 (N12E40 L=  1)   354     87        1 ( 413/ 989- 443/1083) -78.7 E40-E38 N09-N14
      5  231 1014 (N12E59 L=342)   377     64        1 ( 216/ 978- 244/1047) -76.0 E60-E57 N09-N13
      6 1883  898 (S15W66 L=106)   995    109        2 (1841/ 855-1908/ 941)  11.3 W59-W70 S18-S12
      7 1242  784 (S11W11 L= 51)   307    104        0 (1217/ 741-1274/ 824)  45.0 W08-W13 S13-S09
      8 1312  684 (S19W14 L= 55)   410    148        0 (1256/ 634-1376/ 717)  26.6 W10-W18 S21-S17
  

    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