Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  788 1445 (N34E13 L= 51)   635    149        2 ( 724/1408- 866/1466)   5.7 E17-E07 N31-N35
      2  469 1323 (N28E38 L= 26)   440    162        0 ( 450/1251- 492/1394) -60.9 E41-E35 N23-N32
      3 1238 1348 (N23W18 L= 82)   210     81        1 (1202/1325-1273/1365)   9.5 W14-W20 N22-N24
      4 1005  686 (S15E04 L= 60)   217    123        1 ( 973/ 657-1045/ 725)  38.7 E06-E01 S17-S13
      5  617  599 (S18E32 L= 33)   360     65        3 ( 587/ 570- 642/ 628)  45.0 E33-E29 S20-S16
      6  385  553 (S21E53 L= 11)   629     95        2 ( 368/ 508- 411/ 594)  50.2 E56-E50 S23-S17
      7  607  419 (S31E39 L= 25)   372     64        1 ( 584/ 395- 637/ 441) -76.0 E41-E35 S33-S29
  

    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