Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  798 1688 (N40E19 L=170)   329     70        0 ( 762/1666- 827/1706) -31.0 E21-E16 N37-N41
      2  492 1685 (N40E45 L=144)   757    145        2 ( 442/1658- 539/1705)   0.0 E50-E41 N37-N41
      3  226 1222 (N09E57 L=132)   656    133        2 ( 205/1173- 253/1275) -54.5 E59-E53 N05-N12
      4  172 1128 (N03E63 L=126)  2014    352        5 ( 122/1021- 208/1214) -73.3 E70-E58 S02-N08
      5 1872 1049 (N01W61 L=250)  1027    192        0 (1833/ 973-1913/1127) -45.0 W56-W67 S03-N05
      6 1005  856 (S14E01 L=188)   368    119        0 ( 975/ 788-1044/ 917)  66.8 E02-W01 S17-S10
      7  334  735 (S21E48 L=141)   384    159        1 ( 291/ 686- 372/ 777)  32.0 E52-E44 S23-S18
      8  512  559 (S32E37 L=152)   467     87        2 ( 469/ 538- 558/ 579)   0.0 E40-E33 S33-S31
  

    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