Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  645 1723 (N58E13 L=246)  3411    528        5 ( 440/1617- 839/1829)   6.8 E39-W13 N54-N62
      2 1544 1663 (N27W57 L=316)   854     97        5 (1518/1624-1570/1708)  32.0 W53-W63 N24-N29
      3 1613 1598 (N22W58 L=317)  2437    227       10 (1523/1535-1695/1676) -19.7 W48-W65 N17-N26
      4 1185 1536 (N29W23 L=282)   541    204        0 (1116/1459-1250/1603) -76.0 W21-W24 N22-N34
      5 1527 1514 (N19W45 L=304)   914    275        4 (1458/1399-1586/1603)  85.9 W42-W48 N11-N25
      6 1673 1518 (N16W58 L=317)   477     93        2 (1632/1491-1709/1553) -51.3 W55-W61 N13-N18
      7  390 1433 (N42E39 L=220)   725    197        2 ( 309/1407- 476/1454)  -3.2 E48-E30 N39-N43
      8  274 1396 (N41E53 L=206)   654     60        1 ( 245/1367- 301/1428) -45.0 E57-E48 N39-N43
      9 1673 1305 (N04W48 L=307)   343     51        0 (1643/1284-1705/1327) -45.0 W45-W49 N02-N05
     10  572  776 (S02E32 L=227)   446     88        3 ( 554/ 722- 592/ 823)  68.2 E33-E30 S05-N00
     11 1617  751 (S27W33 L=292)   304     45        1 (1592/ 725-1638/ 772) -45.0 W30-W34 S28-S24
     12 1341  540 (S33W09 L=268)  1841    534        1 (1128/ 436-1514/ 694)   7.9 E06-W23 S37-S27
     13 1538  471 (S44W25 L=284)   609    150        0 (1498/ 415-1590/ 528)  36.9 W20-W29 S48-S40
     14  421  422 (S20E60 L=199)   784    109        0 ( 380/ 389- 476/ 452) -45.0 E64-E56 S23-S17
     15  369  399 (S21E69 L=190)  1087     64        1 ( 336/ 379- 400/ 422)   3.6 E79-E63 S22-S19
     16 1466  406 (S48W18 L=277)   366     56        2 (1437/ 385-1497/ 420) -31.0 W15-W20 S49-S45
     17 1073  153 (S58E31 L=228)   484     62        1 (1040/ 138-1105/ 170) -45.0 E35-E26 S60-S56
  

    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