Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  830 1634 (N47E13 L=144)   335     46        1 ( 803/1614- 858/1657)   0.0 E14-E09 N45-N48
      2  903 1600 (N44E06 L=150)   301     50        1 ( 882/1576- 927/1626)  63.4 E07-E03 N41-N45
      3  696 1543 (N40E23 L=133)   411     79        0 ( 683/1500- 711/1585) -80.5 E24-E21 N37-N43
      4 1095 1545 (N38W09 L=166)   297     50        1 (1066/1527-1126/1560) -14.0 W07-W11 N36-N39
      5  506 1409 (N31E37 L=119)   457     60        2 ( 469/1384- 538/1431)   0.0 E40-E34 N29-N32
      6  223 1306 (N24E66 L= 90)   515    115        0 ( 212/1258- 235/1370) -66.8 E68-E63 N21-N28
      7 1489  339 (S48W40 L=196)   948    113        1 (1459/ 305-1523/ 387) -34.7 W34-W47 S51-S42
      8 1365  331 (S47W26 L=182)   483     56        1 (1337/ 314-1398/ 350) -23.2 W22-W29 S48-S44
  

    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