Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  815 1772 (N52E10 L= 68)   826    133        1 ( 753/1728- 888/1813)  24.8 E15-E02 N48-N55
      2  754 1638 (N41E15 L= 64)   822    207        1 ( 662/1613- 849/1673) -17.4 E22-E06 N38-N44
      3  518 1483 (N32E32 L= 46)   749    154        1 ( 473/1449- 576/1512) -31.0 E36-E26 N28-N34
      4 1173 1317 (N15W12 L= 90)   205     47        0 (1142/1305-1204/1329)   0.0 W09-W13 N13-N15
      5  268 1284 (N21E51 L= 27)   421    136        1 ( 226/1264- 308/1308) -20.6 E55-E47 N19-N23
      6  725 1215 (N12E17 L= 62)   207     67        0 ( 688/1194- 767/1237)  11.3 E19-E13 N11-N13
      7  146 1016 (N06E65 L= 13)   567    117        0 ( 123/ 978- 164/1057)  29.7 E69-E62 N04-N08
      8  791  801 (S13E16 L= 62)   359     63        0 ( 772/ 764- 809/ 842) -90.0 E17-E14 S14-S10
  

    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