Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1091 1540 (N31E08 L=295)   276     54        0 (1071/1509-1114/1566)  90.0 E09-E07 N28-N32
      2  267 1348 (N02E62 L=242)   389     85        2 ( 250/1310- 285/1386) -53.1 E63-E58 N00-N04
      3 1566 1325 (N28W29 L=333)   830    137        2 (1500/1277-1613/1363)  12.5 W24-W33 N25-N30
      4  204  953 (S23E60 L=243)   843    166        0 ( 149/ 894- 262/1019)  32.0 E69-E53 S27-S17
      5  464  917 (S20E34 L=270)   383     84        1 ( 425/ 890- 506/ 948)  39.8 E36-E30 S22-S17
      6 1784  823 (N03W57 L=  0)   825    103        2 (1745/ 773-1823/ 861)  45.0 W53-W60 N00-N05
      7  881  765 (S20E03 L=300)   451    141        2 ( 811/ 744- 951/ 788)  20.6 E07-W01 S22-S18
      8  257  733 (S36E55 L=248)   922    102        2 ( 210/ 709- 304/ 770)   0.0 E62-E48 S37-S34
      9 1355  572 (S21W32 L=336)   237     72        2 (1320/ 560-1395/ 587)  31.0 W29-W34 S22-S19
  

    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