Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1144 1597 (N41W12 L=293)   320     75        2 (1113/1570-1178/1626) -45.0 W10-W14 N39-N43
      2 1013 1437 (N30W01 L=281)   300     49        1 ( 987/1408-1039/1462)  76.0 E01-W02 N28-N32
      3 1756 1280 (N16W56 L=337)   558     92        1 (1711/1243-1795/1308) -29.7 W51-W59 N13-N17
      4 1498 1198 (N13W32 L=312)   684    166        2 (1433/1151-1557/1249) -45.0 W27-W35 N09-N16
      5  235  869 (S05E59 L=222)  1215    219        5 ( 163/ 821- 311/ 955) -23.6 E67-E51 S07-N00
      6  663  853 (S05E23 L=257)   427    147        1 ( 593/ 812- 712/ 878) -20.6 E28-E20 S08-S03
      7  703  603 (S22E23 L=257)   407    154        0 ( 669/ 552- 732/ 657) -90.0 E26-E21 S25-S17
      8  308  600 (S22E59 L=222)   681    130        1 ( 275/ 572- 356/ 642) -29.7 E62-E53 S24-S19
      9  583  602 (S22E32 L=248)   591    198        2 ( 534/ 567- 638/ 641)  14.0 E35-E27 S24-S18
     10  875  462 (S32E13 L=267)   270     61        0 ( 853/ 443- 899/ 483)  36.9 E15-E11 S33-S30
     11 1301  421 (S37W19 L=300)   756    138        0 (1271/ 358-1341/ 475)  71.6 W17-W21 S42-S33
  

    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