Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1659 1362 (N18W44 L=253)   728    189        1 (1576/1311-1728/1418) -32.5 W38-W49 N14-N21
      2  500 1295 (N14E34 L=174)   651    206        0 ( 413/1264- 596/1316)  -4.1 E40-E26 N12-N15
      3  195 1192 (N09E61 L=147)  3831    286       11 ( 148/1080- 249/1305) -51.8 E69-E54 N02-N16
      4  696  834 (S14E21 L=187)   940    286        0 ( 645/ 712- 728/ 968)  69.4 E25-E18 S21-S05
      5 1277  795 (S17W16 L=224)   231     81        1 (1237/ 780-1306/ 826)  33.7 W13-W17 S17-S15
      6 1681  451 (S39W60 L=269)   418     36        1 (1658/ 433-1704/ 469)   9.5 W56-W64 S40-S37
      7  694  456 (S39E27 L=181)   284     79        0 ( 653/ 445- 741/ 467)   0.0 E30-E22 S39-S37
  

    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