Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  467 1330 (N34E30 L=178)  1036    176        2 ( 393/1287- 531/1370) -28.6 E36-E25 N29-N36
      2  115  919 (N15E70 L=138)   759    193        1 (  90/ 883- 145/ 948)  16.7 E75-E64 N13-N16
      3  564  763 (S02E33 L=175)  2116    433        6 ( 462/ 650- 707/ 838)  15.5 E43-E24 S06-N02
      4 1072  775 (S12E03 L=205)   204     83        0 (1023/ 763-1112/ 789) -31.0 E05-W00 S13-S10
      5  708  581 (S15E30 L=178)  2111    459        6 ( 598/ 441- 777/ 732)  33.7 E42-E21 S21-S07
      6 1059  678 (S17E06 L=202)   197     56        0 (1028/ 662-1083/ 701)   0.0 E07-E04 S18-S16
      7  486  580 (S11E45 L=163)   318     72        2 ( 465/ 544- 507/ 609)  53.1 E46-E42 S12-S08
      8  858  488 (S24E24 L=184)   555    144        2 ( 803/ 450- 928/ 519)  11.3 E28-E18 S25-S22
      9  497  393 (S22E53 L=155)   400     58        2 ( 467/ 377- 532/ 410) -36.9 E55-E50 S23-S20
  

    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