Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  576 1754 (N26E59 L=230)  4139    365       11 ( 455/1672- 699/1857)  83.3 E70-E51 N19-N36
      2  939 1754 (N36E32 L=257)  2176    412        3 ( 806/1618-1019/1818)  -8.1 E41-E19 N28-N41
      3  613 1654 (N20E47 L=242)   502    131        1 ( 591/1599- 633/1713) -37.9 E52-E43 N16-N23
      4  427 1611 (N13E60 L=229)  1200    162        4 ( 389/1546- 473/1667) -33.7 E69-E54 N08-N16
      5  183 1260 (S12E69 L=221)  4883    388       12 ( 103/1122- 266/1375)  57.1 E84-E58 S21-S04
      6  201 1181 (S17E63 L=226)  2109    168        8 ( 162/1092- 241/1243)  63.4 E69-E58 S22-S12
      7 1669 1070 (N16W38 L=327)   768    289        0 (1628/ 964-1711/1203) -47.7 W32-W43 N11-N22
      8  253  833 (S37E53 L=236)  1170    167        3 ( 194/ 781- 303/ 907)  33.7 E63-E46 S41-S31
      9  454  814 (S33E32 L=258)   343     57        1 ( 428/ 794- 484/ 847) -33.7 E33-E28 S34-S31
     10 1806  787 (N06W59 L=349)  1171    149        3 (1753/ 725-1855/ 836)  63.4 W54-W63 N01-N09
     11 1213  793 (S13W17 L=307)   575    253        1 (1117/ 748-1311/ 834)  22.6 W12-W24 S17-S10
     12 1867  756 (N07W69 L=359)   753     76        4 (1841/ 716-1892/ 787) -76.0 W65-W72 N04-N08
     13 1153  449 (S34W29 L=318)   366    107        1 (1099/ 436-1210/ 467)  35.5 W25-W32 S36-S31
  

    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